Search Documents

By E. H. Lee, F. J. Fayers

Problems in reservoir analysis can usually he cupressed in terms of a system of nonlinear partial difjerential equations. A rnethod for .setting physically reasonable boundary conditions for these systerrls by ayplication of the theory of characteristics is discussed. A short description i.c. given of the technique for finding chnmcteristics of general first order and second orrler equations. Specific examples are quoted jrorn eqrrations occurring in theoretical reservoir studies. The princip1e. ossocicred wit11 the use of chnracteristic to set physically realistic boundary conditions me developed, and thr method is discussed in terms of boundary condition for the noncupillary, compressible, two-phase flow problen: jor the cupillary, ir~cowlpressible, two-phase flow problems for the tlrtergent flooding prohlenl; as for the compressible. noncapillary, two-phase flow pro11letn. INTRODUCTION In most problcms in reservoir analysis the physical situation is described by a system of nonlinear partial differential equations. The general solution to the set of equations usually comprises a wide range of functions ancl a particular solution is normally selectetl by applying suitable boundary conditions to the problem. In some cases it is clear from the physical situation what form the boundary conditions should take, but it is usually desirable to have an independent mathematical check of the suitability of the chosen conditions. From the mathematical standpoint, a system of partial differential equations with boundary conditions is said to be properly set if the combination determines a unique solution which is continuously dependent on the equation and the boundary conditions. Such a solution is sometimes termed "reasonable". In investigating the properties of the mathematical representation, it is useful to determine the structure of the system by the application of the theory of characteristics. This theory determines in the space of the independent variables, lines or surfaces which play a special role in thc system's behavior. The characteristics determine the possibility and nature of the propagation of discontinuities in the derivatives of the dependent variables. Along the characteristics certain differential relations hold which for some problems allow a so- lution to be obtained analytically. In the case of hyperbolic equations the characteristics can bc used as 3 basis for a numerical solution. If a complicated nonlinear system is to be solved, a study of behavior of the characteristics will indicate what boundary conditions should be chosen. By using the method of characteristics, several advances have rccently been made in solving reservoir engineering problems. Sheldon, Zondek and Cardwell' have discussed how the the method may be used grap-phically to solve the incompressible two-phase displaccmerit problem when the gravity term is included in the equations. Fayers and Perrine' have applied the method to study the equations representing detergent flooding. When compressibility and solution-gas effects are introduced into the two-phase flow problem, the equations and their characteristics become more complicated. It will be seen that in this problem the characteristics are useful for indicating a reasonable set of boundary conditions but cannot readily be used to compute a numerical solution. This problem has been solved on the IBM 701 using normal finite difference methods hy West, Garvin and Sheldon." The material to be prcscntcd in this report is divided into three sections. The first discusses the determination of the characteristics of first and second order partial differential equations; the second describes the selection of suitable boundary conditions by the application of the principle that the future cannot affect the past; the third discusses the use of characteristics for setting the boundary conditions of the two-phase solution-gas flow problem. A useful introduction to the method of characteristics has been given by Hildebrand.' A summary of the method will be given in this section of the paper. The First Order Equation The quasi-linear equation of first order is It is linear in the partial derivations of z. The properly linear equation has the form The characteristics of are curves along which the

By K. C. Hong, R. B. Jensen

The problem of determining the optimum set of steam volumes and cycle lengths for a single well undergoing multicycle steam stimulation in order to maximize the cumulative discounted net income has been formulated mathematically and programmed for a digital computer. The mathematical fonnulation of the problem and the method for its solution are discussed in this paper. The oil production performance during each stimulation cycle was simulated by either a constant percentage or a harmonic decline. Using simple analytical expressions for production performance, the cumulative discounted net income and cumulative time of operation were related to the pertinent process and cost parameters and two principal process control variables (cycle length and steam injection volume). The discrete maximum principles was used to transform the equations for cumulative time of operation and cumulative discounted net income into a set of simultaneous equations. The simultaneous equations then were solved by trial and error on a digital computer to determine the set of cycle lengths and stem injection volumes that gives maximum cumulative discounted net income over the project life. INTRODUCTION Steam stimulation is a process for improving the oil recovery rate from wells producing high-viscosity crudes. The process is applied on an individual well basis and is executed in a series of cycles, each consisting of three phases: steam injection, soaking (steam condensation), and production. Significant increases in production rate following the stimulation operation result from heating the reservoir around the wellbore. As heat is removed with produced fluids and dissipated into nonproductive formations, the production rate declines, usually to near the prestimulation value. Typical production responses are given in case histories reported by Owens and Suter,l and are depicted in Fig. 1. Duration of the production phase is equal to the time for the oil production rate to decline to some specified value and is called the cycle length. (Cycle length is defined as the producing portion of the cycle and does not include the downtime required for steam injection and soaking operations.) Termination of the production phase coincides with the start of steaming for the next cycle. The process is continued, cycle by cycle, until it becomes unprofitable. Models 2-4 have been developed to simulate behavior of a single well during one cycle of steam stimulation and have been used to investigate the effects of system and operating conditions on the production responses. However, these investigations have not shown how the models can be used to determine a most profitable set of operating conditions for a multicycle stimulation project. Perhaps these models are too complicated mathematically to be adapted readily to an optimization study. This paper presents a method for optimizing a multicycle steam stimulation project based on simple production performance models. First, the performance of a steam-stimulated well during each cycle was simulated by either a constant percentage or a harmonic decline model. Using these simple analytical expressions for production performance, the cumulative discounted net income (before tax) and cumulative time of operation for each cycle were related to the pertinent process and cost parameters and two process control variables—the cycle length and the steam injection volume. Finally, the profit function to be maximized was formed. The analogy between the above problem and problems of optimizing multistage decision processes that have been reported extensively in the chemical engineering literatures5-10 in recent years has led this investigation to the use of the discrete maximum principle.5-7 This principle was used to transform the equations for cumulative discounted net income and cumulative time of operation into a set of simultaneous equations characterizing the optimum operating conditions in terms of the two control variables. These equations

Jan 1, 1970

By J. C. Martin

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

Jan 1, 1969

By H. J. de Haan, L. Schenk

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

Jan 1, 1970

By 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 B. H. Caudle, B. M. Hickman, I. H. Silberberg

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

Jan 1, 1969

By R. W. Parsons

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

Jan 1, 1967

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

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

Jan 1, 1970

By 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 G. Thodos, R. J. Smialek

A method has been developed that permits the correlations of Mehra and Thodos8 and Dastur and thodos2 for binary systems to be used for the prediction of the vapor-liquid equilibrium constants of the light and heavy constituents, Kland Kh, of a ternary hydrocarbon mixture. Furthermore, the unique relationship observed in this study, i.e., a plot of log Ke/Ki vs log Kp/K h is linear and passes through the origin, permits the calculation of the K-value for the intermediate component. This approach has been used to establish K-values for the constituents of the three ternary hydrocarbon systems used to develop this method and has also been applied to two additional ternary hydrocarbon systems to yield an average deviation of 7.7 per cent resulting from a total of 294 points for the five ternary systems considered. INTRODUCTION The rigorous calculation of vapor-liquid equilibrium constants for systems containing two or more components is a difficult problem and becomes extremely complex at elevated pressures and for conditions approaching the critical point. Although thermodynamic relationships can be developed to account for the vapor-liquid equilibrium behavior of simple systems, their use is limited because the necessary PVT data for these systems are not generally available. To establish a method for calculating K constants directly, without any recourse to the PVT behavior of a system, Mehra and Thodos8 applied the principle of corresponding states to the vapor-liquid equilibrium constants of binary hydrocarbon systems. They developed a series of charts which made possible the prediction of K-values for all binary hydrocarbon systems not containing methane. This work has been extended by Dastur and Thodos2 to include binary hydrocarbon systems containing methane. The charts developed in these two studies cover the high-pressure region and are applicable to conditions in the vicinity of the critical point. The ability of this approach to predict vapor-liquid equilibrium constants of binary hydrocarbon systems suggests that this method may be extended to the prediction of the K-values of the components of ternary hydrocarbon systems. In this investigation, the charts developed for both the methane and non - methane binary systems,2,8 together with mixture rules developed in this study, have been used to predict the K-values of the light and heavy components of a ternary hydrocarbon system. These calculated values are then used to establish the K-value of the intermediate component through the use of a novel relationship developed in the present study. DEVELOPMENT OF METHOD FOR PREDICTING K-VALUES OF TERNARY SYSTEMS The application of the corresponding state charts developed for binary hydrocarbon systems is restricted to those systems containing two components. The extension of these correlations to three-component systems requires that the ternary system be reduced to an equivalent binary, which can be treated in a manner consistent with the corresponding state development. To establish such an approach, the vapor-liquid equilibrium constant of the light component K is first calculated by assuming that the ternary mixture consists of the light component and a pseudo-heavy component of the intermediate and heavy constituents of the original ternary. To establish the vapor-liquid equilibrium constant of the heavy component Kh, the light and intermediate constituents of the ternary are combined to form a pseudo-light component of another equivalent binary system. The calculated values Kl and Kb are then utilized to calculate the vapor-liquid equilibrium constant of the intermediate component EQUILIBRIUM CONSTANT FOR LIGHT COMPONENT K To calculate the K value of the light (most

Jan 1, 1965

By J. J. Cosgrove, J. E. Warren

A general model which approximates the effect of cross- flow has been developed to give a practical method for predicting the waterflood behavior of a stratified reservoir. The model is based on a modification of Dietz's theory and allows for variations in both the permeability and hydrocarbon pore volume. In the particular cases considered, it is assumed that the permeability can be characterized by a log-normal distribution and the hydrocarbon pore volume by a normal distribution. A simple graphical method which enables the practicing engineer to predict the behavior of a stratified system is presented The results obtained by the proposed method are compared with those obtained by the Dykstra-Parsons method. As a result of this study the following conclusions have been drawn: (1) The effect of cross-flow in a stratified system can be appreciable, particularly at very favorable or very unfavorable mobility ratios; (2) Under normal conditions, the effect of variations in the hydrocarbon pore volume can be neglected; and (3) The failure to use all of the available permeability data can lead to large errors in the prediction of the behavior of a stratified reservoir. INTRODUCTION Various methodshavebeenproposedto characterize and to predict the waterflood behavior of a stratified reservoir. Most of the methods assume that the reservoir is composed of discrete homogeneous continuous layers. With this model, the degree of stratification can be measured by several parameters based on core-analysis data. Among these are the Lorenz coefficient1 and the variation of a log-normal permeability distribution? The behavior of stratified systems is usually predicted by the Stiles, Dykstra-Parsons method or some modification of these. In both of these methods the reservoir is divided into discrete homogeneous layers with no cross-flow between layers. In the Stiles method the mobility ratio is assumed to be equal to unity while in the Dykstra-Parsons method it is allowed to take on any value. Other predictive methods have been proposed; method allows for cross-flow between the beds and a method by Schmalz and Rahme6 correlate.s the recovery directly with the Lorenz coefficient. In this new approach, which is essentially a continuous analog of Hiatt's method, the effects of both mobility ratio and cross-flow between the beds have been included. It is assumed that the permeability can be represented by a log-normal distribution and the hydrocarbon pore volume or porosity by a normal distribution. These types of distributions have been observed by several authors, 1,2,4 and most field data seem to confirm their observations; e.g., if these distributions are assumed and there is a 1:l correspondence between porosity and permeability samples, the commonly encountered linear relationship between porosity and the log of permeability is obtained. Since, in many field cases the permeability and poro sity data are truncated or cut-off at predetermined upper and/or lower values, the effect of discarding part of the data was investigated with the proposed method. A better technique for truncating the data is suggested. THEORY In the derivation of the method to be described in this paper the following assumptions are made: 1. Capillary forces are negligible per se; their effects are only manifest in the relative permeability curves. 2. The fluids are immiscible, incompressible and homogeneous. 3. The reservoir is horizontal and uniformly thick; it is initially liquid saturated i.e., connate water and oil. 4. The reservoir is composed of discrete layers, each having its own permeability, porosity and connate water saturation; each layer is homogeneous,

Jan 1, 1965

By W. C. Hauber

Techniques previously published to predict the production performance of a water flood when the mobility ratio is not unity have been primarily restricted to a five-spot well pattern. When other types of well patterns are considered, the prediction of performance has required either a model study or a complex numerical solution of a multiple-fluid system with a high-speed digital computer. Simple methods are developed in this report for predicting the areal sweep efficiency and injectivity history of a water flood for arbitrary well patterns and mobility ratios. The effect of an initial gas saturation can be incorporated into this method of prediction. With the layer superposition concepts published by Prats, et al.,1 these methods can be used to predict the production performance of a water flood in a reservoir having a wide range in permeability. A reasonably good agreement is shown by comparisons of the results calculated by this method with those obtained either by model studies or by complex numerical solution of multiple-fluid systems. Therefore, it has been concluded that a reasonably accurate method has been developed for predicting the production performance of a water flood employing an arbitrary type of well pattern for any given mobility ratio. A medium-speed digital computer can be used in applying this method. INTRODUCTION Production performance of multi-fluid five-spot water floods can be predicted by the techniques published by Prats, et al.1 This technique takes into account the effects of initial gas saturation, mobility ratio and vertical variations in horizontal permeability. A basic assumption in applying this technique is that unit displacement efficiency does not vary with time after the passage of the flood front. Prats, et al. demonstrated that this assumption is valid for normal waterflood applications with small mobility ratios. To extend Prats' technique to other types of well patterns has required either an experimental model study or a numerical solution of a multiple-fluid system to describe the performance of an individual homogeneous layer. Sheldon and Dougherty2 presented an exact numerical solution employing any arbitrary well pattern. This technique is complex to apply, requiring (1) numerical solution of the single-fluid system, (2) conformal mapping of the streamlines and iso-potential lines to form a new grid network, and (3) a numerical solution of the multi-fluid system on the new grid network. A method was developed in this paper to approximate areal sweep efficiency and injectivity as functions of time for an individual homogeneous layer for arbitrary well patterns and mobility ratios. Analytic expressions were derived to directly apply this method to five-spot and direct line drive well patterns. For other types of well patterns, this method is similar to the approach of Higgids and Leighton,6 in that it utilizes the stream tubes and iso-potential lines obtained by the solution of differential equations when the mobility ratio is unity. However, this method is easier to apply than the Higgins-Leighton method since it does not require knowledge of the variation between oil and water relative permeabilities and saturation. Using the relationship between areal sweep efficiency, injectivity and the volume of water injected into an individual homogeneous layer as obtained by the method presented in this paper, and Prats' technique for the superposition of individual layers, one can predict the production performance of any water flood. ASSUMPTIONS GENERAL Basic assumptions made in this report are that fluids are incompressible and the layers are horizontal, homogeneous and uniformly thick. It is assumed that as water is injected into the layer, three distinct regions are formed —a water bank, an oil bank, and an unflooded region. Each region is assumed to displace other regions in a pis-tonlike manner, with only one fluid flowing in each region, so that there is no change in unit displacement efficiency with time after the flood front has passed in the reservoir. Therefore, the residual oil saturation in the water bank and the oil saturations in the unflooded region are constant throughout the layer. The ratio (F) of the bulk volume of a layer occupied by the oil and water banks at any time before oil breakthrough, to the corresponding bulk volume of a layer occupied by the water bank, is given' by the equation BREAKTHROUGH AREAL SWEEP EFFICIENCY To predict the behavior of any type of waterflood pattern, one must know the breakthrough areal sweep efficiency Eb(Mc,o, F) for the pattern for any water:oil mo-

Jan 1, 1965

By T. D. Mueller, F. C. Miller, Jr. Earlougher R. C.

In analyzing pressure buildup tests for field wells producing both oil and gas, the common practice is to use a modification of single-phase flow theory. Validity of such an approximation has been demonstrated for single-well solution-gas-drive systems. This paper indicates that such approximations are also valid for two-well solution-gas'-drive systems. From this it is inferred that the technique can be used for multiple-well systems. A computer was used to simulate the behavior of a two-well solution-gas-drive reservoir to test the validity of the above type of analysis. Simulation results indicate that pressure buildup tests in such a system can be analyzed within engineering accuracy for formation permeability and pressure. A rule of thumb is given for estimating the length of time a well must be shut in for the pressure in a nearby producing well to increase significantly. To observe such an increase, the shut-in well would have to be left shut in much longer than normal. INTRODUCTION An important technique for obtaining data concerning a producing petroleum reservoir is the pressure buildup test. Such a test, when properly conducted and analyzed, provides information on the average reservoir pressure and the permeability in the major drainage area of a well. The theory upon which the analysis of a pressure buildup test is based assumes that the behavior of the fluid in the reservoir is adequately described by the diffusivity equation.1 Use of the diffusivity equation implies that a single fluid of small and constant compressibility is flowing. To analyze a pressure buildup test in situations which involve more than one fluid phase, the single-fluid analysis is extended. This paper verifies the extension of pressure buildup analyses to two-phase, two-well systems. Miller, Dyes and Hutchinson1 have presented a technique for analyzing pressure buildup tests in circular bounded reservoirs. Their method requires that the reservoir be producing at pseudo-steady state (constant pressure gradients) prior to shut-in. An alternate reservoir model, presented by Horner,2 assumes that the reservoir is infinite. The Homer model does not assume a Pseudo-steady state prior to shut-in, but the assumption that the reservoir is infinite implies that the average pressure can be determined only if the total production prior to shut-in is small compared to the total fluid originally present. Limitations imposed by the pseudo-steady state and infinite reservoir assumptions are avoided by the technique proposed by Matthews, Brons and Hazebroek.3 Their method can be used to calculate both the Permeability and the average pressure in bounded single-well systems which are not necessarily at steady state. They also suggested determining the average pressure of a multi-well reservoir producing from pseudo-steady state by volumetrically averaging the static pressures of the individual wells. Matthews and Lefkovits4 used numerical simulation techniques to verify that this method does produce adequate results for single-phase, multiple-well systems. Perrine5 proposed modifications to the single-phase theory so that pressure buildup tests in multiple-phase systems could be analyzed. He suggested replacing the single-phase mobility (k/µ) by the sum of the mobilities of the individual phases, and replacing the fluid compressibility by an average compressibility weighted by the saturations of the separate phases. By using numerical techniques to solve the two-phase flow problem for a single-well system, he concluded that this approximation gives results within engineering accuracy. More recently, Weller6 has shown that, for a gas-oil system with a single well at the center of a circular reservoir, pressure buildup tests can be analyzed adequately by using the modification proposed by Perrine. Feller's results also indicate that, as the gas saturation increases, this analysis becomes less accurate. To date, no evidence is available to show that the Perrine approximations can be combined with the Matthews-Brons-Hazebroek technique to calculate the average pressure of a multi-well,

By R. C. Earlougher, F. G. Miller, T. D. Mueller, H. J. Ramey

There are many studies of flow in radial systems that can be used to interpret unsteady rerervoir flow problems. Although solutions for systems of infinite extent can be used to generate solutions fur finite ow systems by superposition, application is tedious. In this paper a step is made toward simplifying calcu1ations of such solutions for finite flow systems. Superposition is used to produce a tabulu-lion of the dimension1e.s.s pressure drop function at several locations within a bounded square that has a well at its center. The square system provides a useful building block that may be used to generate flow behavior for any rectangular shape whose sides are in integral ratios. Values of the tabulated dimensionless pressure drop function are simply added to obtain the dimensionless pressure drop function for the desired rectangulm system. The rectangular system may contain any number of wells producing at any rates. Furthermore, the outer boundaries of the rectangular system may be closed (no-flow) or they may he at constant pressure. Mixed conditions also may be conyidered. Tables of the dimensionless pressure drop function for the square system are presented and various applications of the technique are illustrated. Introduction In 1945, van Everdingen and Hurst' published solutions for the problem of water influx into a cylindrical reservoir. Since this problem is mathematically identical with the depletion of a cylindrical reservoir with a well at the origin, the van Everdingen-Hurst solution may be used to study the depletion problem. In their analysis, they assumed that the fluid had a small, constant compressibility such that flow was governed by the diffusivity equation For a constant production rate q starting at time zero, van Everdingen and Hurst showed that the unsteady pressurc distribution for both finite and infinite systems could be expressed in terms of a dimensionless pressure 'Sabulations of the dimensionless pressure drop for a unit value of r, were provided by van Everdingen and Hurst,' and later by Chatas.' Others also presented values in graphical or tabular form." If the radius of the well becomes vanishingly small, r,+ O, the line source solution may be used for Eq. 2 when infinite systems are considered. Eqs. 5 and 6 are excellent approximations for Eq. 2 under certain conditions: In 1954, Matthews, Brons and HazebroekV emonstrated that solutions such as Eq. 5 can be superposed to generate the behavior of bounded geometric shapes; i.e., the behavior of a bounded single-well system can be calculated by adding together the pressure disturbances caused by the appropriate array of an infinite number of wells producing from an infinite system. These wells are referred to as image wells. Matthews, Brons and Hazebroek considered systems containing a single well producing at a constant rate. This superposition can be represented analytically as

Jan 1, 1969

By H. K. van Poollen, H. C. Bixel

A treatment is given of the transient pressure behavior of a well located at the center of a circular region surrounded by a radial discontinuity. On either side of the discontinuity, the values of permeability, viscosity, compressibility and porosity are uniform but may be different from those on the other side of the discontinuity. The results are obtained by solving a pair of finite difference equations. The numerical solution to these equations is obtained using a digital computer. The results are compared with previously published, approximated analytical solutions to the same problem. Solutions, presented graphically, show pressure decline for constant rate fluid production. The range of variables studied include dimensionless time from '0.001 to 100 and storage capacity ratio from 0.001 to 1,000. The well behaves as if it were in an infinite reservoir for dimensionless times less than 0.25. Reservoir properties near the well can be estimated in the usual manner. An overlay technique is used to match an experimental curve with one of the theoretical curves. It is possible to estimate the distance to a discontinuity by substituting the actual time t and the corresponding dimensionless time tD at which a match occurs into the equation a = (0.000264 tk1/tD , µ1 cI)I/2 where kI/ I µI cI is the diffusivity near the well, and may be estimated from data taken at early time. Several buildup curves are computed. These curves show that for early shut-in times, correct values for transmissibility are obtained from conventional analysis. However, erroneous values of static reservoir Pressure are obtained unless data at large shut-in times are taken. INTRODUCTION A mathematical treatment of the transient pressure response of a well located at the center of a region bounded by a circular discontinuity is given. Within a region (Fig. 20) the properties of both the rock and the fluid are considered to be constant on either side of the discontinuity. While these properties are considered constant, they may be different on opposite sides of the discontinuity. A discontinuity of this type could be a fluid-fluid contact or a sudden change in rock characteristics such as thickness, porosity or permeability. Analytical solutions to this problem are available. 1-3 However, they are so involved that they are of little practical use. An approximate solutionl is available but the range of times over which it is valid has not been specified. A numerical solution to a pair of finite difference equations is used to obtain the solution given in this paper. PREVIOUS WORK One of the first solutions to the problem was published by William Hurst.1 Hurst considered unsteady-state flow of fluids through two sands in series with different mobilities in each sand. He used Laplace transforms to obtain a solution for a single well located at the center of the circle enclosing the first of the two sands. In this case, the solution for the pressure change in Region I (Fig. 20) is2

By D. M. Beeson, G. D. Ortloff

Water-propelled banks of carbon dioxide recovered both high- and low-viscosity crude oil substantially in excess of that recovered by water flood in linear flow model experiments. The increase in oil recovery is explained by oil swelling and viscosity reduction. Additional oil was recovered by solution gas drive at the end of the carbon dioxide floods when the model systems were depres-sured. This study was an investigation of the use of a water-driven carbon dioxide bank to recover crude oil. Techniques involving the injection of carbon dioxide have aroused con-siderable interest in recent years among methods considered for increasing oil recoveries. A number of investigators have studied the process in the laboratory1-' and in the field.7,5 Laboratory studies referenced indicate that carbon dioxide in conjunction with water flooding recovers more oil than does water flooding alone. These re- ports diverge widely, however, in stating magnitude of improvement. We believed that a water-driven carbon dioxide bank would increase recovery of crude oils by two mechanisms: viscosity reduction and oil swelling. Due to reduction in the oil viscosity caused by solution of carbon dioxide, displacement efficiency should be improved. In addition, the residual oil should be highly swollen by the carbon dioxide held in solution. Thus, the final residual oil saturation, in terms of stock-tank oil, would be lower than in the absence of carbon dioxide. To test this hypothesis, experi-

By A. B. Dyes

In the operation of a water-driven reservoir, a free gas saturation can he established by maintaining production rates fast enough to cause the reservoir pressure to decline below the bubble point. The benefit of such a procedure on the displacement efficiency of the oil by water is illustrated for two types of sandstone samples from one reservoir. Those rock samples showing the poorest recovery to water drive in the absence of a free gas saturation give the most improvement in the presence of a free gas saturation. Employing one type of reservoir fluid (low .shrinkage) the benefit of evolved gas on oil displacement is calculated at several reservoir pressures below the bubble point. This demonstrates the presence of air optimum pressure for water displace,rnent which is .several hundred pounds per .square inch below the bubble point. Displacement of oil at the optimum pressures results in the removal of 7 to I? per cent more oil than would be obtained by displacement at the bubble point pressure. In some instance. this magnitude of increase may be worthy of consideration in establishing the maximum efficient rate of production from many water-driven fields. INTRODUCTION For many years the belief that the use of gas could be combined with that of water to result in a better recovery operation than the use of either of these phases alone has been prevalent in the oil industry. This is evidenced by the several field trials which have employed the injection of gas and water in combination.1.2.v In most cases the small amount of any additional recovery, in combination with the limited knowledge available on the older reservoirs, has left uncertain the influence of gas on recovery in these earlier attempts to combine the use of gas with water. Research studies in recent years have demonstrated through laboratory. floods the magnitude of the benefit of the presence of gas on the displacement of oil from sandstones by water.5,6,7,8,9,10,11 In addition, studies of the distributions of the phases which are present during the displacement process have developed concepts concerning the mechanism by which this benefit is obtained."." It now remains for the engineer to incorporate these findings into the general technology of reservoir operations. It is the purpose of this paper to illustrate that the production of a water-driven reservoir at pressures below the bubble point is one such operation whereby improved oil displacement efficiency can be obtained as a result of gas evolution. The establishment of a gas saturation by evolution of dissolved gas was one of the techniques used by several investigators in the laboratory studies. The possibility of obtaining the benefit of gas evolution in water-driven reservoirs by producing at pressures less than the bubble point has 3150 been pointed out previously." The desirability of this operation has undoubtedly

Jan 1, 1955

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 L. H. Reiss, J. Groult, L. Montadert

Many fields, where the reservoir is composed of sandy layers, show great complexity because of the lack of continuity which results from a particular type of seditnentation. This complexity may be a factor of prime importance in studying reservoir drainage, especially in the displacement of oil by water or gas. This paper describes work done to study this heterogeneity and its influence in the exploitation of the fields of the lower Devonian sandy shale in the Fort Polignac Basin in Algeria. Two lines of research were followed. 1. On the border of the hasin, a detailed sedimentological study established rules for the distribution of the sandy zones, which furnished hypotheses necessary for devising a geologic model of the reservoir. A certain number of sandy zones were recognized and characterized by their sedimentation. 2. In the interior of the basin, production logs deterrnined the productivity profile. Study of these point data led to a description of the distribution of oil and gas phases in the reservoir so that the continuity of different zones between the wells could be deduced. These techniques were applied to the study of the lower Devonian sandy shale reservoir of Zarzaitine. The sedimentological study guiding the correlation work, with the aid of the conventional static data (logs and cores), led to the establishment of a useful geologic model of the reservoir. Production logs, thanks to the contribution of dynamic elements, permitted corrobora-tion of the model. In particular, the logs confirmed separation, as proposed by the sedimenfological study, of two comparable sandy groups by revealing the continuity of a thin impermeable layer, invisible on conventional logs. INTRODUCTION The heterogeneity of sandstone reservoirs considerably influences the distribution of fluids, as well as their displacement in the reservoir. Rational development of a field and maximum recovery of hydrocarbons are closely tied to an exact and precise knowledge of the intern architecture. We will try to show how it is possible I refine the knowledge of the reservoir by seeking qui different and varied sources of information. Productic geology, the group of techniques which leads to idei concerning reservoir architecture and its geometric descris tion, furnishes the framework of this undertaking. Of these techniques. those used so far are qui varied. Descriptive as well as quantitative, the technique partake as much of classical geology as of reservo science. They include (1) geologic studies of the enti basin, especially a stratigraphic and sedimentologic study of the reservoir member at its outcrop, detail1 observations of cores taken from wells and, finall oriented cores; (2) interpretation of geophysical dat (3) interpretation of electric logs, nuclear and sonic lo and dipmeters: (4) petrophysical study of the reserve on samples taken from wells; (5) reservoir dynamics (stu of the production history and, in particular, interpretati of production logs); and (6) analysis of the fluids a thermodynamic studies. This incomplete list could also include (1) study the migration of the fluids, (2) geochemistry (chemi characterization of the oils), (3) regional hydrodynamii and (4) injection and detection of tracers. If most of these techniques are classic, attention i nevertheless be called to applying two of the techniql in conjunction: the study of the reservoir member its outcrop and the study of dynamic distribution fluid phases in the reservoir through use of product logs, used in the study of the sandy shale reservoirs the lower Devonian of the Fort Polignac Basin, Alge in particular, their application to the study of the 2 zaitine field where they have led to a new and satisfact model of the reservoir will be described. THE ZARZAITlNE FIELD Situated in the interior of the Fort Polignac 8 about 150 kilometers (km) north of the lower Devol outcrop, the Zarzaitine Devonian oil field was discovc On Ian. 28, 19'' (Figs. 1 and l8). It is in the fern a monocline truncated by two perpendicular fault: great throw, with a downthrown section to the north\: Erosion at the top of the structure has caused one of it to disappear. It is a reservoir with a gas cap an m coefficient on the order of 1:5, open to the aql On the northeast, with the oil zone covering an Of 106 sq km and with initial oil reserves of 222 mi

Jan 1, 1967