Search Documents

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 P. B. Crawford, W. L. Huskey

A study has been made to determine the effects of random vertical fractures existing in a reservoir matrix on the effective matrix properties and producing characteristics of a well. The study indicated that short vertical fractures may distort the isopotentials by 1 percent in the vicinity of the fracture. A correlation was developed which showed the effect of fracture density on the producing capacity of a well. The correlation indicated that fracture shape had little effect on producing capacity, but that total fracture length was closely correlated with the producing capacity. Effective values of reservoir permeability were found to correlate closely with fracture density when measured parallel to the direction of flow. INTRODUCTION In reviewing the petroleum literature, it is found that in most analytical solutions of hydrocarbon reservoir behavior the assumption is made that the rock is uniform throughout although core analysis and well logging surveys indicate that substantially all reservoirs are heterogeneous and some are fractured. Only in recent years have investigators studied the effects of reservoir hetrogeneity. Russell1 discussed the effect of fractures in reservoirs. The literature indicates that fractures greatly affect the exploitation of man reservoirs. This is described by Elkins and Skov2 and Little-field et al. 3 Several types of reservoir heterogeneity have been studied by unsteady-state methods. In 1959, Landrum et al. 4 showed that transient phenomena in hydrocarbon reservoirs could be studied by thermal models. pickering5 studied the flow of heat in stratified linear reservoirs consisting of plates of different metals separated by sheets of rubber which were suddenly exposed to a low-temperature source at one end. Cotman6 studied both laterally and vertically heterogeneous thermal model reservoirs with special reference to cross flow. Another type of reservoir heterogeneity involves the presence of highly permeable, circular lenses around or near the wellbore which was studied by Miesch7 using a thermal and a steady-state model. Many studies have been made of the effects of induced hydraulic fractures on well performance. 8-10 Work regarding vertical fractures that exist initially in the reservoir was reported by Givens, et al.11 Since there appears to be little published work on the importance and influence on reservoir performance of fractures that exist out in the rock matrix, this study was conducted to determine the effects of vertical fractures on pressure distribution, producing capacity and the effective permeability of the reservoir rock. The analogy between the flow of electrical current through a conducting body and the flow of fluid in a porous medium provides a convenient method of studying various aspects of petroleum reservoirs. For electrical radial steady-state flow,12 For steady-state fluid flow, With the above analogy, it is only necessary to establish the ratio of the values of the electrical units in terms of the corresponding flow units to make the analogy quantitative. The ratios are shown below.

By H. N. Hall

Reservoir and producing characteristics can govern the decision to use either a one-, two- or three-dimensional model for making predictions for gravity-drainage reservoirs. Examples of conditions requiring one-, two- and three-dimensional calculations are given. In 1961 the author presented a method for predicting one-dimensional gravity-drainage performance. This work has been exrendcd to obtain a three-dimensional model, which is described and for which a sample problem is presented. Introduction Many papers on gravity drainage or gas-cap drive'-3 have been presented for conditions where the reservoir is treated as a single areal segment so that, in effect, only one-dimen-sional flow is considered. Other authors'-Q ave given evidence showing that the one-dimensional concept would be unsatisfactory to predict reservoir performance for certain types of reservoirs. Experience has indicated that specific reservoir and producing conditions govern whether one-, two- or three-dimensional concepts must be used as a basis for reliable reservoir performance predictioils. Previously the author presented a paper outlining a method for predicting the performance of gravity-drainage reservoirs.' That paper treated the reservoir as a single areal segment so that it was applicable for only one-dimensional flow. The work has been extended to obtain a three-dimensional model that can be used to predict performance of gravity-drainage reservoirs. Although capable of being three-dimensional, it can be used equally well for one- and two-dimensional predictions. This paper discusses those reservoir and producing conditions governing the choice of model, and also describes a three-dimensional model suitable for predicting reservoir performance in the more complex situations. Discussion The Choice of Model To predict reservoir performance, it is necessary to represent the reservoir by a physical model. A mathematical model, based on the assumed physical model, then is used to predict the movement of oil and gas within the reservoir. The word "model" is used here to connote the combination of both a physical model and a mathematical model based on the physical model. Certain types of reservoir situations may be predicted suitably with a one-dimensional model, whereas two- and three-dimensional models may be required for other situations. This section illustrates reservoir conditions associated with the use of a specific type of model. Many types of reservoir structures are such that gravity drainage possibly could be an important factor in oil production. Production practices, as well as the nature of the structure. can influence the type of model that should be used. obviously, it is impossible to discuss all combinations of reservoir structure and production practices that might be encountered in gravity-drainage reservoirs. The following are typical examples of conditions where a choice must be made among a one-, two- or three-dimensional model to predict reservoir performance satisfactorily. Conditions Requiring a One-Dimensional Model Fig. la shows contours for an asymmetrical anticline in which gravity drainage might be expected to be important. Fig. Ib depicts a cross-section along Line A-A in Fig. la. Existence of a gas cap is a good clue that gravity drainage will be an important factor in producing oil from this type of reservoir. However, it is not essential since a gas cap will form as oil is produced from the reservoir. If the proper combination of reservoir permeability and withdrawal rates is present, the reservoir will behave like a

Jan 1, 1969

By W. A. Klikoff, I. Fatt

The concept of fractional wet wattability is examined. Fractional water wettability of a reservoir rock is defined as the fraction of the internal surface urea that is in contact with water. Capillary pressure and relalive permeability of unconsolidated sand are shown to be functions of fractional wettability. INTRODUCTION The petroleum industry has long recognized that wettability of reservoir rock has an important effect on multiphase flow of oil, water and gas through reservoirs. As early as 1928 the American Petroleum Institute sponsored a study of wettability as part of API Project 27 at the U. of Michigan.' Despite 30 years of research, there is still little exact knowledge of the wettability of reservoir rocks. There are two parts to the wettability problem. After agreeing to a uniform nomenclature in regard to wettability,' the first question to be answered is, "What is the in situ wettability of a given reservoir rock?" If this can he answered the next question is, "What part does wettability play in determining the characteristics of multiphase fluid flow through the rock?" This paper represents an oblique attack on the problem of wettability. No attempt is made here to answer the basic question of wettability in situ. instead the consequences of the concept of fractional wettability are examined. Multiphase flow in sandpacks is shown to he highly influenced by fractional wettability. Jennings' has given 3 definition of wettability and the other terms used in discussing wettability. These terms must be applied to the physical situation existing In reservoir rock. A survey of the pertinent literature from 1928 to 1956 indicates that the concept of a contact angle was applied to reservoir rock in the same way it would be applied to a flat, homogeneous surface. Attempts were made to state quantitatively the wettability of a reservoir rock in terms of a contact angle which was presumably constant at all points on the very rough and heterogeneous interior surface of a porous rock. Calhoun, et al,3-6 prepared synthetic consolidated and un-consolidated porous media in which they claimed there was a known uniform contact angle. They then showed the effects on the capillary pressure and relative permeability characteristics of varying this angle. The API Project 47 at the U. of Texas' and others' have made extensive studies of an indirect approach to the contact angle through the use of heat of wetting data. Even if successful. however, this approach also states the wettability of porous rock in terms of a contact angle which is uniform over the entirc surface. If the angle varies from one part to another on the internal surface, there is no way of determining from the measurements the area distribution of contact angles. In 1956 Brown and Fatt5 suggested that the concept of a contact angle, as applied to reservoir rock, be abandoned. This suggestion was made because it is known that the internal surface of most reservoir rocks is composed of many different minerals, cach with a different surface chemistry and a different capacity to adsorh surface active materials from reservoir fluids. Furthermore, the operation of a contact angle in determining the form of a fluid-fluid interface is difficult to picture in the very complex geometry of a pore. Brown and Fatt proposed that the wettability of reservoir rock be stated in terms of the fractional internal surface area that is in contact with water or oil. All surfaces on which there is water are called water-wet; surfaces on which there is oil are called oil-wet. The fractional water wettability is then stated as a number which represents the fraction of the internal surface that is in contact with water. A symmetrical statement can be made for the fractional oil wettability. The concept of a fractional wettability as previously stated has in its favor the recognition of the heterogeneous mineral composition of most reservoir rocks. Another point in its favor is that fractional wettability can be measured quantitatively with relative ease. Hol-brook and Bernard10 se a simple dye adsorption test to obtain fractional wettability of reservoir rocks. Amott11 uses a combination of imbibition and displacement to arrive at a wettability index of reservoir rocks which seems related to fractional wettability in the range 0.25 to 0.75 fractional water wettability. Jennings" has shown the changes in relative permeability that take place when a porous material is changed from unity to zero fractional water wettability. He also shows that reservoir rocks in their natural state, but at room temperature and atmospheric pressure, have relative permeability characteristics which would indi-

By H. J. Ramey

Conventional calculation of total system isothermal compressibility for a system containing a free gas phase involves, among other things, evaluation of the change of oil and gas formation volume factors and the gas in solution with pressure. Preferably, this information should be obtained from laboratory measurements made with particular oils and gases. Often, experimental measurements are not available. In this case, it is necessary to obtain pressure-volume-temperature relationships from general correlations such as those of Standing' for California oils. In order to speed estimates of compressibility, generalized plots have been prepared of the change of both oil formation volume factors and gas in solution, with pressure from Standing's correlations. A generalized plot for estimating the change in the two-phase (oil and gas) formation volume factor with pressure is also presented. Usually, the effect of gas dissolved in reservoir water upon the total system compressibility is neglected for gas saturated systems, due to the low solubility of gas in water. Results of this study indicate that the increase in total system compressibility caused by solution of gas in water is often as large as the compressibility of wafer, and can be magnitudes larger for low pressure systems. Generalized results for estimating the change of gas in solution in water with pressure are presented in tabular and graphical form. INTRODUCTION During the past decade, pressure build-up and drawdown techniques have gained an important place in reservoir engineering. Build-up and drawdown analyses are only two special applications of the broad field of transient fluid-flow theory. All solutions of transient fluid-flow problems contain a parameter called the total system isothermal compressibility. This property of fluids and porous rock is a measure of the change in volume of the fluid content of porous rock with a change in pressure, and it may vary considerably with pressure. Evaluation of total system isothermal compressibility is not difficult, but it is tedious and time-consuming. Often compressibilities are estimated roughly, or transient flow methods are neglected completely. The benefits of using accurate system compressibility in properly-executed build-up or drawdown analyses are: 1. Better planning of pressure build-ups may be achieved to avoid unnecessary loss of revenue due to excessively long shut-in periods, or to shut-in periods too short to yield useable data. 2. Better and more reliable estimates of static formation pressures for reserves estimates and rate performance estimates. 3. Reliable information for evaluation of well completion effectiveness, and planning and interpretation of well stimulation efforts. The purpose of this paper is to clarify the nature of the total system isothermal compressibility, and to present useful methods for estimation of compressibility, particularly for systems containing a gas phase. DEVELOPMENT Numerous publications have presented solutions to transient single-phase flow of slightly compressible fluids, stressing pressure build-up applications. In transient flow, a compressibility* term arises to permit volume content of fluids in porous rock to change as pressure changes. The basic nature of the compressibility term is usually taken for granted. Problems arise in practical applications of transient fluid theory because most published works consider only one flowing fluid—in an ideal porous system containing only one fluid. In 1956, Perrine2 presented an intuitive extension of single-phase flow pressure build-up methods to multiphase flow conditions. Later, Martin- established conditions under which Perrine's multiphase build-up method had a theoretical foundation. Perrine has shown that improper use of single-phase build-up analysis in certain multiphase flow situations can lead to gross errors in estimated static formation pressure, permeability and well condition. It is likely that, much pressure build-up data for oil wells should be analyzed on the basis of multiphase flow. For both single-phase and multiphase build-up analysis, the isothermal compressibility term in dimensionless time groups often should be interpreted as the total system compressibility. All real reservoirs contain one or more compressible fluid phases. In addition, rock compressibility can contribute in an important way to the total system compressibility. The proper total system compressibility expression may contain terms for compressibility of oil, gas, water, reservoir rock and terms for the change of solubility of gas in liquid phases.

Jan 1, 1965

By W. L. Martin, J. N. Dew, H. B. Steves, M. L. Powers

This paper presents and discusses the results obtained during a pilot test in the Loco field in southern Okla homa. The test was conducted in a 2%-acre pattern that was part of a 20-acre conventional waterflood pilot area. The conventional flood was well past its economic limit when the hot waterflood was initiated to obtain technical information and operating experience. Temperature data from nine observation wells showed that about 75 percent of the pattern was affected by heat, and that heat losses were severe in the pilot pattern. About 60 percent of the injected heat was lost to overburden and underburden zones during hot-water injection. Wellbore heat losses were held at tolerable levels at the shallow depth of the test by providing a low-pressure air annulus between the injection tubing and well casing. Hot water provided water injectivity increases of 200 to 400 percent. The hot water channeled across the lower portion of the oil sand in three directions through zones of relatively high water saturation. There was no conclusive evidence that natural fractures or pressure partings aflected the flow of fluids in the pilot pattern; however response undoubtedly was aflected by localized pressure gradients and by injection-production rate ratios. The results showed that hot waterflooding increased oil recovery in a reservoir containing 600 cp oil. The total tertiary oil production from the pilot pattern area was 3,896 bbl or about 156 bbl/acre-ft swept by heat. The corre.sponding WOR was about 34:I. Description of the Hot Waterflood Process To hot waterflood an oil reservoir, water is injected that has been heated to a temperature substantially higher than the original reservoir temperature, but lower than the vaporization temperature of water at the prevailing pressures. In the reservoir the hot water flows continuously into cooler sand and rapidly loses heat to the sand until it has been cooled to the original reservoir temperature. Thus, a heated zone and a region or "bank" of cooled water begins to accumulate immediately after hot water injection is started. This bank of cooled water continues to grow ahead of the heated zone, which also grows, but at a slower rate. This occurs because heat transfer is almost instantaneous, and the ratio of heat capacities of water to rock is such that two or three unit PV of hot water must be injected to heat a given unit bulk volume of the reservoir. The primary displacement mechanism is the same for both hot and conventional cold waterfloods; i.e., "piston" displacement occurs at the original reservoir temperature. The incremental benefits of hot waterflooding usually occur long after the breakthrough of cold water at producing wells, and the increased oil recovery necessarily is accomplished with high WOR's (water-oil ratios). Heat decreases the viscosity and density of oil and water. These effects result in more rapid and increased recovery of secondary or tertiary oil. If the cost of the required heat is low enough, the ultimate oil recovery of a hot waterflood should be increased substantially over what would be expected at the economic limit of a conventional cold waterflood. The economic benefits of any hot-fluid injection project depend primarily upon the cost of the heat required to produce more oil at an increased rate. This cost depends in part upon the amount of heat lost to surrounding formations. Heat loss depends upon reservoir thickness, water injection rate and temperature, the depth of the formation, well spacing, and the characteristics of the reservoir rock and surrounding formations. In general, percentage heat losses decrease as injection rate and reservoir thickness increase. Although it is an old idea, hot-water injection has not received widespread field application in the oil industry as a drive process. Much of the original field work with hot water was done in Pennsylvania fields where water permeabilities and injection rates are low. In these cases, hot-water injection was used primarily as a means of increasing in-jectivities — not as a recovery process. Ramey' recently has published an excellent review of the development of hot fluid injection processes. Recent publications"' and our own experiences now indicate that hot water or steam injection also can effectively increase the oil recovery from reservoirs containing viscous crudes. However, the economics of these methods as displacement processes have not been established. Several theoretical predictive techniques have been published (for a review see Ramey'), but simplifying assumptions make

Jan 1, 1969

By H. H. Rachford, J. Douglas, D. W. Peaceman

A numerical solution of equations describing two-phase flow in porous media shows promise in providing a technique for predicting the displacermet from satlds of oil by water or gas. The description includes the influence of relative permeability, fluid viscosities and densities, graviry, and capillary pressure, and, tholrgh tested ortly for two-dimensional cases, should be eqlmlly applicable in three-dimensional geometry. Two memorical methods ore presented: the first method is quite general in its applicability; the second method can be used only with certain types of boundary conditions but requires less computing. Comparisons of computed results with data from Inhoratory models are presented. These data were taken on a water flood of a strarified model and on water floods of a five-spot model for favorable and unfavorable mobility ratios. On the stratified model, excellent agreement with recovery at breakthrough was obtained; agreement with recovery after breakthrough was poor. In the five-spot model, good agreement was obtainecl with recovery at and after breakthrough. The purpose of this paper is twofold: first, to present a reservoir engineering method requiring knowledge only of rock geometry and the normally measured rock and fluid properties for calculating the multi-dimensional flow of water displacing oil from porous, water-wet rock containing connate water; and, second, to investigate the validity of the method by comparing results of calculations with previously observed displacements in laboratory models. The availability of such a verified technique of reservoir analysis would afford major advantages. First, it would demonstrate that the fundamental macroscopic concepts of two-phase fluid mechanics, i.e., relative permeability and capillary pressure, yield an adequate description of the physical process. Also, since most reservoir rocks appear to be water-wet and contain connate water along with oil, it would provide a method of technical and economic value for calculating the course of oil displacement by water directly from measured reservoir and fluid properties. It is important to emphasize that the calculation is subject only to the limitations of detail and precision of reservoir information, but not to limitations introduced by simplifying assumptions. Further, it would providc almost the only means of examining the influence of factors such as the size and extent of reservoir inhomogeneities and uncertainties in the basic reservoir data. Knowledge of this influence is essential in stating quantitatively the detail required to define a reservoir and in establishing the relation between the uncertainty in the definition of reservoir properties and the reliability of predicted performance. Moreover, in the process of solving a particular problem much detailed information would become available about the displacement process. For example, at each point in the reservoir for all stages of the displace ment. the calculations would yield not only the water and oil saturations but also the direction and magnitude of fluid velocities and the local fluid mobilities. Such dctail is potentially of great value in giving insight into the mechanics of particular displacements. The development of such a method has long been the goal of research in the application of numerical an-alysis to petroleum reservoir enginecrinz. The most complete treatment of the displacement process published to date is that of Douglas. Blair, and Wagner,' which, however, was limited to flow in a single dimension. The method presented in the present work is hascd on the numerical solution of a finite difference analogue of the multi-dimensional differential system describing the displacement process. Although current work has considered only displacement of oil by water from water-wet sands, the differential system for other immiscible displacements, such as gas displacing oil with which it is in phase equilibrium, is quite similar. It should therefore be expected that the technique described will be applicable to displacement by gas as well. The description of the simultaneous flow of fluids through porous media in terms of relative perrncability and capillary pressure has been adequately discussed in the literature and standard textbooks. Sec, for exampje, Muskat. Chapter VII.'

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 M. W. Clegg

The application of the Laplace transformation to problems in the flow of compressible fluids in porous media has provided a large number of exact solutions. For plane radial flow, however, these solutions are either complex integrals or infinite series and are of little value to the field engineer. In the case of production at constant well pressure, the available approximate solutions me valid for large times only. In this paper it is shown that an approximate inversion formula for the Laplace transform, developed for the solution of viscoelastic problems, is applicable to radial flow problems and provides simple analytical solutions to constant terminal pressure problems. The method may be used to obtain approximate solutions to many problems, including media with radial permeability discontinuities, multi-layer formations and pressure buildup in wells after shut-in. The results are compared with the few available computer solutions as well as the large time solutions, and it is shown that this approximate method greatly extends the time interval over which a simple analytical solution is acceptable. INTRODUCTION The study of transient problems in the flow of fluids through porous media has benefited greatly from the application of transform methods. The use of the Laplace transformation for solving parabolic equations has been widely discussed in the field of heat conduction and diffusion as well as in the petroleum literature. Removal of the time variable with the Laplace transformation generally reduces the problem to a boundary value problem which may be solved by standard techniques. A much more formidable problem then faces the engineer, however, for frequently the transform does not possess a simple inverse. The result is that the general inversion integral must be used and this leads to either an infinite integral or an infinite series, both of which are difficult to handle from a computational standpoint. Asymptotic approximations for the inverse have been known for some time and these yield approximate inverse functions that are valid for very large or very small times — but frequently the times of interest lie somewhere between these two extremes. Therefore, some acceptable approximation valid over a larger interval of time is desirable. During the past few years a number of methods for achieving this have been developed and some of these are discussed briefly in this paper. The relative merits of the various methods are not evaluated here, but some general conclusions reached by other authors are given. One of these methods has been applied to problems associated with the radial flow of compressible liquids to producing wells. In the case of production at constant well pressure, the method leads to simple analytical solutions for a number of standard problems; e.g., homogeneous formation, permeability discontinuities, pressure buildup. These solutions greatly extend the range of validity of the asymptotic ones (valid for large times only) and should be of value in studying the behavior of wells producing under constant pressure conditions. METHOD FOR THE APPROXIMATE INVERSION OF LAPLACE TRANSFORMS Much of the work that has been done on the approximate inversion of Laplace transforms has originated in viscoelastic problems. Apart from the asymptotic approximations,1 one of the first approximate methods was due to Alfrey.2 If (s) is the Laplace transform of f(t) defined by

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 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 F. J. Fayers, R. L. Perrine

For many years the problem of increasing ultimate recovery of oil from a reservoir has been a subject of interest to the oil industry. At present, a standard secondary recovery technique is to flood the reservoir with water once natural forces of the system are depleted. The permeability characteristics of the rock, together with interfacial forces which exist between the two fluids and the porous material. allow only part of the oil in the rock formation to be recovered by the water-flooding process. Laboratory experiments have shown that when detergents are added to the flooding water, the flow properties of oil in a porous rock arc improved. This improvement results in greater oil recovery. Even though higher recoveries are obtained using detergents, the extra cost associated with the method may make it economically unsatisfactory. To establish economic feasibility requires experimental laboratory measurement of fluid and rock properties together with a calculational procedure to predict field results. Pilot tests in the field may also be required. This paper presents a mathematical procedure for predicting detergent flood behavior. A number of reasonable simplifying assumptions are made to obtain a tractable system of equations. The resulting method may be regarded as a generalization of Buckley-Lev-crett' theory to solve simultaneously for water saturation and detergent concentration profiles. Consequently this type of calculation can be used in much the same manner as Buckley-Leverett theory is used to evaluate a water flood. The mathematical method employed is the so-called method of characteristics for solving hyperbolic partial differential equations. This application of the method is instructive because several other problems in reservoir analysis may be analyzed in a similar manner. FLUID FLOW EQUATIONS The flow behavior of oil and water without detergent through porous reservoir rock is usually expressed by Darcy's law and a continuity equation for each phase. If gravitational forces are neglected, vrg = kr (s)µr v p rg.........(1) vrg = kr (s)µr v p rg.........(1) ? pa va = F ?/?t (sd) where v is the volumetric flux, k represents permeability, S is saturation, and P the pressure. Fluid viscosity is denoted by 11, p is the density, 4 is porosity, and I the time. The subscripts w and o refer to water and oil. The capillary forces in the system lead to a difference in pressure between phases. This may be expressed by the relation: Pm -P0 =Pc(S).........(5) The capillary pressure? P,, is determined experimentally as a function of saturation. DETERGENT BALANCE EQUATION Addition of detergent to the water phase affects the flow properties of the system. Capillary pressure and relative permeability become functions of detergent concentration, c, as well as fluid saturation. That is: P--P (S,c) and k,,,.- - -- k,,,. At low oil saturation values, the presence of detergent reduces capillary pressure and reduces the relative permeability ratio. Detergents are surface-active and adsorb on rock surfaces. The quantity adsorbed, a, will depend on the concentration c of detergent solution and may also depend on the saturation; thus, a = a(S,c). In this analysis we have assumed that the detergents used are water soluble but not oil soluble, and their movement is therefore restricted to the water phase. The dispersion of detergent in water is neglected by comparison with transfer by flow within the system. Detergent movement is determined by the water flux and adsorption: ? cv m = -F ?/?t (cs + a) We will assume that the quantity of detergent in solution is too small to significantly change water density or viscosity. When the permeability and capillary pres-

By C. Chu

A theoretical investigation has been made of the forward combustion process using a one-dimensional linear mathematical model, taking into consideration the effect of the vaporization-condensation which occurs on the leading edge of the heat wave. This work involves the solution of five coupled partial differential equations. Besides the vaporization-condensation phenomenon, these equations account for conduction, convection, combustion, heat loss, diffusion and bulk fluid flow. For the one-dimensional linear model studied, the vaporization-condensation phenomenon does not induce appreciable change in the temperature at the combustion front; and its primary effect is to create a steam plateau and to increase the length of the heated zone ahead of the combustion front. This effect becomes more pronounced at lower pressures, higher porosities or reduced gas saturations. The peak temperature and the temperature profile on the leading edge of the heat wave stabilize after a certain period. The length of the steam bank remains practically constant, although the length of the water bank increases as the heat wave advances. INTRODUCTION The existence of the vaporization-condensation phenomenon in the heat-wave process and the important role played by the phenomenon have been recognized by several investigators. Kuhn and Kochl stated that steam plateaus were frequently observed on the temperature records of their experiments. The steam plateaus were attributed primarily to the vaporization and subsequent condensation of the interstitial water existing in the oil sand. Szasz2 suggested that both lighter hydrocarbons and water are vaporized on the leading edge of the heat wave, carried forward in the gas stream, and then condensed to create banks of oil and water. Martin et al suggested that the vaporization-condensation phenomenon is one of the main mechanisms of the heat-wave process, along with thermal expansion and viscosity reduction. Wilson et al. reported the existence of a steam plateau several inches in length in their small-scale tube-run experiments. However, this important phenomenon has never been taken into consideration in the numerous theoretical analyses by various authors. 5 The purpose of this work was to study the thermal aspects of a linear heat wave, taking into consideration the vaporization-condensation on the leading edge of the wave, to determine the effect of this phenomenon on the temperature profile of the reservoir, and to investigate how this effect varies when other process variables are changed. THEORY We consider a reservoir of porous medium of cross-sectional area A, extending from x = 0 to x = L. This reservoir contains, aside from the solid matrix itself, a gas phase and a "combined liquid phase" which is a combination of two immiscible liquid phases — namely, an oil phase and a water phase. The oil present in the reservoir is assumed to consist of three fractions, a noncondensable gas, a nondistillable residuum, and a vaporizable oil fraction which may be present in both the gas phase and liquid phase. Before the heat-wave process begins, preheating has taken place and has imparted an initial temperature distribution T,(x) to the reservoir. At the start of the process, a stream of oxygen-bearing gas is introduced through the face at x = 0. This gas supports the combustion of the residual fuel and supplies the heat throughout the process. In order to develop a tractable system of equations, the model properties are postulated as follows: (1) the system is isotropic and homogeneous; (2) there is a uniform fuel content throughout the reservoir; (3) combustion takes place in a moving vertical planar region of infinitesimal thickness; (4) the temperatures of the fluids and the solid matrix are identical; (5) radiation is negligible; (6) the gas and liquid saturations are uniform throughout the reservoir and remain unchanged

Jan 1, 1965

By N. E. Truitt, D. G. Russell

The transient pressure behavior of a well which produces a single compressible fluid through a singte-plane wrticat fracture has been investigated mathematically. The fracture is assumed to possess infinite flow capacity, to be of limited mdial extent, and' to penetrate the producing formation completely in the vertical direction. Previous studies of vertically fractured wells have been concerned primarily with production rate performance or semisteady-state pressure behavior. This study was undertaken to ascertain the influence of vertical fractures on transient pressure tests such as pressure build-ups and flow tests. In a vertically fractured system, flow in the region nearest the fracture is practically linear, whereas farther away from the fracture essentially radial flow prevails. Thus, transient pressure analyses based on radial flow theory are somewhat inaccurate. As fracture penetration increases radially, kh values calcutated from pressure build-up and flow test curves become increasingly larger than true values. Failure to consider the effect of fracture penetration also introduces inaccuracies into the catculation of fracture length from the apparent skin factor and into the determination of average reservoir pressure. If the total length of the fracture is 20 per cent, or greater, of the drainage radius of the well, corrections must be made to pressure build-up and flow test results. Methods for correcting such results are discussed in this paper. For wells with prefracturing pressure build-up or flow test data, it is possible to estimate fracture length by comparison with postfracturing build-up or flow test results. In new wells or wells without prefracturing build-up or flow test data, fracture length must be estimated to correct the values obtained from analysis of pressure tests after fracturing. Fracturing efficiency calculations should be made whenever possible to provide an estimate of fracture length. Tables of the dimensionless pressure drop as a function of time and fracture penetration are included in this paper. Using these values should permit analysis of other types of transient pressure behavior in vertically fractured wells. INTRODUCTION Hydraulic fracturing has been used quite successfully for over a decade as a completion and stimulation technique in oil and gas wells completed in low-permeability reservoirs During this period a considerable amount of theory has evolved on the performance of hydraulically fractured reservoirs and on more efficient means of artificial fracturing. Although theory has been developed, no rigorous investigation has been made of pressure build-up and flow test behavior in such wells. Prats et al.1 first discussed the performance of vertically fractured reservoirs for the case of a compressible fluid. Their work was primarily concerned with production performance at constant flowing pressure. These authors also considered large-time (semisteady-state) constant production rate behavior for vertically fractured wells; however, transient pressure behavior at constant rate was not investigated. McGuire and Sikora10 and Dyes, Kemp, and Caudle2 employed an electrical analog to investigate the influence of artificial vertical fractures on well productivity and pressure build-up. They found that fractures which extend beyond 15 per cent of the drainage radius away from the well alter the position and slope of the straight-line portion of the build-up curve. They concluded that these effects must be considered both in the determination of the effective permeability of the formation and in any calculations of final build-up pressure. Although these authors did not undertake an exhaustive study of the influence of vertical fractures on pressure build-up performance, their limited results were quite interesting from the standpoint of the effects they demonstrated. In a more recent paper, Scott- reported the results of an investigation of the effect of vertical fractures on pressure behavior, which was conducted with a heat flow model. Scott's results appear to be consistent with those reported in Refs. 1 and 2. However, the effects of different fracture lengths on performance were not investigated. Pressure build-ups and transient flow tests are among the most diagnostic tools available to the reservoir engineer or production engineer. Since a very high percentage of present-day well completions incorporate the hydraulic fracturing technique, a definite need exists for information on the effect of fractures on transient pressure performance. For these reasons we have undertaken a rigorous study of pressure build-up and flow test behavior in vertically fractured reservoirs. The objectives of this study were (1) to obtain synthetic pressure build-up and flow test

Jan 1, 1965

By A. G. Spillette, R. L. Nielsen

The purpose of this paper is to further the understanding of reservoir response to hot-water injection by desuribing a two-dimensional, mathematical model of the process. Key assumptions are that no gas phase is present, and that the injected fluid reaches thermal equilibrium instantaneously with the reservoir fluids and sand. The resulting system of three partial differential equations is solved simultaneously through the use of a "leap-frog" application of standard alternating direction implicit methods for the solution of the mass-balance equations and the method of characteristics for solution of the energy-balance equation. The utility of the mathematical model is demonstrated by comparing numerical and analytic temperature distributions for hot-water bank injection and by comparing calculated with observed field behavior. Additional calculations show that hot waterfioding can recover significantly more oil than cold waterflooding, and that a hot-water bank recovers, with less energy input, nearly as much oil us continuous injection. Introduction Consistent field success with hot fluid injection processes requires good reservoir description and a thorough understanding of recovery mechanisms. The latter is fostered by a combination of well designed experimental studies and physically sound mathematical modeling. The purpose of this paper is to further the understanding of reservoir response to hot-water injection by describing a two-dimensional mathematical technique, indicating its validity, and demonstrating its utility for studying the effects of reservoir and operating parameters. Published experimental studies of hot fluid injection processes are few." Even if extensive data were available, two considerations discourage exclusive reliance on laboratory or field work. First, because of the severity of scaling requirements," laboratory results must be interpreted with care. Second, because of the one-shot nature of field experiments — injectivity/productivity tests and pilot floods—results seldom are available over the desired range of operating conditions. These factors emphasize the need for devoting attention to mathematical modeling to complement laboratory and field work. Through the use of mathematical models, scaling uncertainties can be bridged and response can be predicted for unique combinations of reservoir and operating conditions. Numerous mathematical models developed during recent years enable us to calculate temperature distributions, thermal efficiency (or conversely, fraction of heat lost), and oil recovery behavior for hot fluid injection. Ramey" and Spillette' recently have provided reviews of these methods. Most models have concentrated on predicting temperature distributions and thermal efficiency; few have been directed at predicting oil recovery.'."1° Although these techniques are useful for estimating effects of reservoir parameters and operating conditions on process performance, results are limited by the assumptions made and by the methods of coupling independently solved fluid-flow and energy-balance equations. A computer-based method is presented for predicting total reservoir response to hot-water injection that obviates most simplifying assumptions. The model simulates fluid flow and heat transfer in two dimensions within a vertical cross-section spanning the oil sand and adjacent unproductive strata. Known numerical procedures are used to solve the governing partial differential equations. The mathematical model handles the effects of reservoir heterogeneity, gravity, capillarity, relative permeability and temperature-dependent fluid properties. In addition, a wide range of operating conditions can be modeled, including hot-water followed by cold-water injection. Mathematical Description of Hot-Water Injection Key assumptions in the mathematical model for hot-water injection are that (1) no gas phase is present and (2) the injected fluid reaches thermal equilibrium instantaneously with the reservoir fluids and sand. Relative permeability and capillary pressure are assumed independent of temperature. Darcy's law for multiphase fluid flow relates the velocities of the oil and water phases to the phase pressure gradients by the following relations:

Jan 1, 1969

By R. A. Morse, W. O. Keller

The properties of limestone reservoir rocks such as the distribution and degree of continuity of the pore systems, and the relative volumes and permeabilities of the systems making up the complex cause large variations between performance of individual limestone reservoirs and their susceptibility to secondary recovery methods. The effects of these factors on the mechanism of fluid flow cannot be adequately evaluated with presently developed concepts, laboratory data, and geological information. The observance and interpretation of the performance of individual limestone reservoirs provides, at present, the most adequate approach for evaluating the integrated effects upon performance of the many. now immeasurable variations in the properties of limestone reservoirs. INTRODUCTION The development and application of techniques for increasing the efficiency of oil recovery from natural reservoirs is a problem of primary importance to the industry. During the past several years, a great deal of effort has been expended by the technical personnel of the industry toward the improvement of present known methods of oil recovery. and evaluating the factors which control the susceptibility of particular reservoir types to economic application of secondary recovery methods. Providing adequate and accurate laboratory data on the properties of the reservoir rock and its contained fluids. together with good production statistics, are available, methods have been evolved for estimating with reasonable accuracy the performance of an oil reservoir under either continued natural depletion or conditions imposed by introducing additional displacing fluid into the reservoir from an extraneous source through the injection of gas and/or water. The susceptibility of limestone reservoirs to secondary recovery by gas injection is very much dependent upon gas-oil relative permeability relationships for the reservoir in question. In a rock formation in which the porosity is of the intergranular type such as is found in sand stones and some non-fractured dolomites, a representative sample of the pore structure can be obtained in a small core plug. In this type of reservoir, it has been shown that relative permeability data obtained from laboratory experiments can be Properly applied to evaluate the flow relationships actually observed during the depletion of a reservoir. In addition, a good idea of the fraction of the reservoir which will be swept by the injected gas may be obtained from the spacing pattern used and the permeability profile of cored wells. From these data, it has been demonstrated that reliable estimates of performance and recoveries under gas injection operations can be made for reservoirs in which the Porosity is of the intergranular type. Limestone reservoirs on the other hand often present much more complex problems. One reason that the performance of limestone reservoirs is so variable is that there may be, in reality. two or more pore systems, the integrated performance of which will depend upon the physical properties of each system and their inter-relationship. For instance, in a limestone reservoir, the main pore structure, as far as void spaces for storage of hydrocarbons is concerned, may be represented by the intergranular openings of the rock. Usually in limes and dolomites, the intergranular pore openings are very small and the matrix permeabilities extremely low. Often, as in the case of Many dolomite and oolitic lime reservoirs, this is the only pore system, and performance is observed to be very similar to that of a sand stone of similar permeability. However, in the great majority of' the limestone reservoirs, there is present, in addition to the intergranular pore system, a complementary system of openings caused by fracturing or solution.',' The openings in the latter system are usually many times more permeable than the intergranular pore system due to their much larger size, and may contribute almost all the fluid carrying capacity of the formation, even though their volume may be a very small part of the total pore space in the reservoir. Obviously, the performance of such a complex under either primary or secondary control will deend upon the distribution and degree of continuity of both systems, the relative volumes and permeabilities of the

Jan 1, 1949

By R. A. Morse, W. O. Keller

The properties of limestone reservoir rocks such as the distribution and degree of continuity of the pore systems, and the relative volumes and permeabilities of the systems making up the complex cause large variations between performance of individual limestone reservoirs and their susceptibility to secondary recovery methods. The effects of these factors on the mechanism of fluid flow cannot be adequately evaluated with presently developed concepts, laboratory data, and geological information. The observance and interpretation of the performance of individual limestone reservoirs provides, at present, the most adequate approach for evaluating the integrated effects upon performance of the many. now immeasurable variations in the properties of limestone reservoirs. INTRODUCTION The development and application of techniques for increasing the efficiency of oil recovery from natural reservoirs is a problem of primary importance to the industry. During the past several years, a great deal of effort has been expended by the technical personnel of the industry toward the improvement of present known methods of oil recovery. and evaluating the factors which control the susceptibility of particular reservoir types to economic application of secondary recovery methods. Providing adequate and accurate laboratory data on the properties of the reservoir rock and its contained fluids. together with good production statistics, are available, methods have been evolved for estimating with reasonable accuracy the performance of an oil reservoir under either continued natural depletion or conditions imposed by introducing additional displacing fluid into the reservoir from an extraneous source through the injection of gas and/or water. The susceptibility of limestone reservoirs to secondary recovery by gas injection is very much dependent upon gas-oil relative permeability relationships for the reservoir in question. In a rock formation in which the porosity is of the intergranular type such as is found in sand stones and some non-fractured dolomites, a representative sample of the pore structure can be obtained in a small core plug. In this type of reservoir, it has been shown that relative permeability data obtained from laboratory experiments can be Properly applied to evaluate the flow relationships actually observed during the depletion of a reservoir. In addition, a good idea of the fraction of the reservoir which will be swept by the injected gas may be obtained from the spacing pattern used and the permeability profile of cored wells. From these data, it has been demonstrated that reliable estimates of performance and recoveries under gas injection operations can be made for reservoirs in which the Porosity is of the intergranular type. Limestone reservoirs on the other hand often present much more complex problems. One reason that the performance of limestone reservoirs is so variable is that there may be, in reality. two or more pore systems, the integrated performance of which will depend upon the physical properties of each system and their inter-relationship. For instance, in a limestone reservoir, the main pore structure, as far as void spaces for storage of hydrocarbons is concerned, may be represented by the intergranular openings of the rock. Usually in limes and dolomites, the intergranular pore openings are very small and the matrix permeabilities extremely low. Often, as in the case of Many dolomite and oolitic lime reservoirs, this is the only pore system, and performance is observed to be very similar to that of a sand stone of similar permeability. However, in the great majority of' the limestone reservoirs, there is present, in addition to the intergranular pore system, a complementary system of openings caused by fracturing or solution.',' The openings in the latter system are usually many times more permeable than the intergranular pore system due to their much larger size, and may contribute almost all the fluid carrying capacity of the formation, even though their volume may be a very small part of the total pore space in the reservoir. Obviously, the performance of such a complex under either primary or secondary control will deend upon the distribution and degree of continuity of both systems, the relative volumes and permeabilities of the

Jan 1, 1949

By Frank G. Miller

This report presents developments of fundamental equations for describing the flow and thermodynamic behavior of two-phase single-component fluids moving under steady conditions through porous media. Many of the theoretical considerations upon which these equations are premised have received little or no attention in oil-reservoir fluid-flow research. The significance of the underlying flow theory in oil-producing operations is indicated. In particular, the theoretical analysis pertains to the steady, adiabatic, macroscopically linear, two-phase flow of a single-component fluid through a horizontal column of porous medium. It is considered that the test fluid enters the upstream end of the column while entirely in the liquid state, moves downstream an appreciable distance, begins to vaporize, and then moves through the remainder of the column as a gas-liquid mixture. The problem posed is to find the total weight rate of flow and the pressure distribution along the column for a given inlet pressure and temperature, a given exit pres5ure or temperature and given characteristics of the test fluid and porous medium. In developing the theory, gas-liquid interfacial phenomena are treated. phase equilibrium is assumed and previous theoretical work of other investigators of the problem is modified. Laboratory experiments performed with specially designed apparatus. in which propane is used as the test fluid, substantiate the theory. The apparatus. materials and experimental procedure are described. Comparative experimental and theoretical results are presented and discussed. It is believed that the research findings contributed in this * paper should not only lead to a better understanding of oil-reservoir behavior, but also should be suggective in regard to future research in this field of study. INTRODUCTION In recent years much time and effort has been consumed in both theoretical and experimental studies of the static and . dvnamic behavior of oil-reservoir fluids in porous rocks. Although lack of sufficient basic oil-field data, principally concerning the properties and characteristics of reservoir rocks and fluids, largely precludes quantitative application of research results to oil-field problems, qualitative application has become common practice. In effect. oil-reservoir engineering research is serving as a firm foundation for oil-field development and production practices leading to increased economic recoveries of petroleum. This province of research. however, still poses many perplexing problems. The thermodynamic behavior of two-phase fluids moving through porous media constitutes one facet of reservoir-fluid-flow research that has not received the attention it deserves. This report embodies a theoretical discussion of this subject and a description of a series of related laboratory experiments. The significance of the problem to oil field operations is indicated but in articular the report centers around a theory and method for analyzing the steady. macroscopically linear, two-phase flow of a fluid (a single molecular species) through a horizontal column of porous medium. For simplicity in showing how the thermodynamic behavior of two-phase fluids moving through porous media affects oil-reservoir performance problems, attention is focused temporarily on a particular well producing petroleum from an idealized water-free solution-gas drive reservoir, the reservoir rock being a horizontal, thin, fairly homogeneous sandstone of large areal extent confined between two impermeable strata. The flowing hydrocarbon fluid is considered to exist entirely as a Iiquid at points in the reservoir remote from the well; however. the decline in fluid pressure in the direction of the well causes vaporization of the hydrocarbon to begin at a radial distance r from the well. Upstream from r the fluid moves entirely as a liquid and downstream from r it moves either entirely as a gas or as a gas-liquid mixture depending on the properties of the hydrocarbon and on the thermodynamic process it follows during flow. The distance r would be variable under transient flow conditions. but for purposes of analysis the flow is considered to l~e steady at the particular instant of observation during the flowing life of the well of interest. If the flow were isothermal and the hydrocarbon a pure substance, the fluid would be entirely gaseous downstream from r. Thus, this isothermal flow process for a pure substance would require that the heat of vaporization be supplied at r. over zero length of porous medium, at the precise rate necessary to maintain the constant temperature. This means that the solid matrix of the porous medium (reservoir rock) and the surroundings (impermeable strata confining the reservoir rock) would have to serve as infinite heat sources. Heat-transfer requirements would be somewhat less severe for the isothermal flow of a multicorn-ponent hydrocarbon as bubble and dew points at the same temperature correspond to different pressures. In this instance isothermal conditions would be sustained without complete vaporization of the fluid over zero length of porous medium. Nevertheless. as the flow is in the direction of decreasing

Jan 1, 1951

By Frank G. Miller

This report presents developments of fundamental equations for describing the flow and thermodynamic behavior of two-phase single-component fluids moving under steady conditions through porous media. Many of the theoretical considerations upon which these equations are premised have received little or no attention in oil-reservoir fluid-flow research. The significance of the underlying flow theory in oil-producing operations is indicated. In particular, the theoretical analysis pertains to the steady, adiabatic, macroscopically linear, two-phase flow of a single-component fluid through a horizontal column of porous medium. It is considered that the test fluid enters the upstream end of the column while entirely in the liquid state, moves downstream an appreciable distance, begins to vaporize, and then moves through the remainder of the column as a gas-liquid mixture. The problem posed is to find the total weight rate of flow and the pressure distribution along the column for a given inlet pressure and temperature, a given exit pres5ure or temperature and given characteristics of the test fluid and porous medium. In developing the theory, gas-liquid interfacial phenomena are treated. phase equilibrium is assumed and previous theoretical work of other investigators of the problem is modified. Laboratory experiments performed with specially designed apparatus. in which propane is used as the test fluid, substantiate the theory. The apparatus. materials and experimental procedure are described. Comparative experimental and theoretical results are presented and discussed. It is believed that the research findings contributed in this * paper should not only lead to a better understanding of oil-reservoir behavior, but also should be suggective in regard to future research in this field of study. INTRODUCTION In recent years much time and effort has been consumed in both theoretical and experimental studies of the static and . dvnamic behavior of oil-reservoir fluids in porous rocks. Although lack of sufficient basic oil-field data, principally concerning the properties and characteristics of reservoir rocks and fluids, largely precludes quantitative application of research results to oil-field problems, qualitative application has become common practice. In effect. oil-reservoir engineering research is serving as a firm foundation for oil-field development and production practices leading to increased economic recoveries of petroleum. This province of research. however, still poses many perplexing problems. The thermodynamic behavior of two-phase fluids moving through porous media constitutes one facet of reservoir-fluid-flow research that has not received the attention it deserves. This report embodies a theoretical discussion of this subject and a description of a series of related laboratory experiments. The significance of the problem to oil field operations is indicated but in articular the report centers around a theory and method for analyzing the steady. macroscopically linear, two-phase flow of a fluid (a single molecular species) through a horizontal column of porous medium. For simplicity in showing how the thermodynamic behavior of two-phase fluids moving through porous media affects oil-reservoir performance problems, attention is focused temporarily on a particular well producing petroleum from an idealized water-free solution-gas drive reservoir, the reservoir rock being a horizontal, thin, fairly homogeneous sandstone of large areal extent confined between two impermeable strata. The flowing hydrocarbon fluid is considered to exist entirely as a Iiquid at points in the reservoir remote from the well; however. the decline in fluid pressure in the direction of the well causes vaporization of the hydrocarbon to begin at a radial distance r from the well. Upstream from r the fluid moves entirely as a liquid and downstream from r it moves either entirely as a gas or as a gas-liquid mixture depending on the properties of the hydrocarbon and on the thermodynamic process it follows during flow. The distance r would be variable under transient flow conditions. but for purposes of analysis the flow is considered to l~e steady at the particular instant of observation during the flowing life of the well of interest. If the flow were isothermal and the hydrocarbon a pure substance, the fluid would be entirely gaseous downstream from r. Thus, this isothermal flow process for a pure substance would require that the heat of vaporization be supplied at r. over zero length of porous medium, at the precise rate necessary to maintain the constant temperature. This means that the solid matrix of the porous medium (reservoir rock) and the surroundings (impermeable strata confining the reservoir rock) would have to serve as infinite heat sources. Heat-transfer requirements would be somewhat less severe for the isothermal flow of a multicorn-ponent hydrocarbon as bubble and dew points at the same temperature correspond to different pressures. In this instance isothermal conditions would be sustained without complete vaporization of the fluid over zero length of porous medium. Nevertheless. as the flow is in the direction of decreasing

Jan 1, 1951

By C. Özgen Karacan

Degasification and methane control in underground coal mining is an important area augmenting conventional ventilation in order to prevent possible fires and explosions due to excessive methane emissions. Coal bed reservoir engineering has been applied for many years to coalbed methane production. In recent years, modeling the flow and draining methane, not only from coal bed itself but also from the overlying strata, have been recognized as important areas of expertise for supplementing underground ventilation. This paper demonstrates the applications of coal bed reservoir engineering and modeling techniques for optimizing and controlling methane emissions in longwall and continuous miner sections. Practical applications and some considerations are presented for reservoir modeling, borehole performances, and face emission predictions for various coal bed and borehole configurations. A section on field reservoir studies and well bore log analyses for methane control are also presented as a complementary way of controlling methane and optimizing ventilation.

Jan 1, 2009