Search Documents

By S. S. Marsden, S. A. Khan

Externally generated foam was injected continuously into short porous media. Both flow rate and pressure drop were measured. Liquid saturation was determined by electrical conductivity. Foam yuality I, expressed as the ratio of gas volume to total volume, was varied from 0.70 to 0.96. As measured with a modified fann VG Meter, apparent viscosity of this foam µa decreases with increasing shear rate but usually falls within the range of 50 to 500 cp. At a given shear rate, µa increases almost linearly with I-. When measured with a Bendix Ultraviscoson, kinematic pa is independent of V but absolute pa increases with T from about 3 to 8 cp. The effective permeability -apparent viscosity ratio ke/µa decreases almost linearly with V for porous media of high permeability, but the rate of decrease becomes less for tighter ones. The relative permeability-apparent viscosity ratio kr/µa us V data does not fall on a single line. The kr/µa ratio increases with liquid saturation in the porous medium and with surfactant concentration. Estimates of pa for foam in porous media vary from 30 to 100 cp. INTRODUCTION Although research on the development of a foam-drive, oil recovery process has been going on for almost a decade, most of the significant publications have appeared within the last several years. This illustrates well the rate at which interest in this process is accelerating. Bond and Holbrook 1 were the first to describe the use of foam to improve oil recovery in their patent of 1958. They proposed that an aqueous foaming agent slug be injected into the formation and that this be followed by gas to produce a foam in situ. Fried 2 studied the injection of foam into porous media which had already been subjected to conventional gas or water drives and found that gas could be used to drive a foam bank which would, in turn, displace additional oil in the form of an oil bank. He attributed the increased oil recovery to the high effective viscosity of foam flowing in porous media. His microscopic observations showed the importance of foam generation and regeneration within the porous medium. By injecting both air and aqueous surfactant solution, Bernard3 generated foams within the porous medium in which oil displacement was being studied. In a separate empirical test, he also measured the dynamic foaming characteristics of the same surfactants in water and/or oil. With some exceptions and for the seyen surfactants studied, there seems to be a qualitative relationship between the efficiency of liquid rlisplacrment and the dynamic foaming testtsed. This relationship was not consistent enough to eliminate the necessity of actual foam flood tests in porous media for surfactant selection. In a study basic to gas storage in aquifers, Bennett4 described the displacement of brine by foam in consolidated porous media. Among other things, he stated that the ability of a surfactant solution to foam is more important than the stability of its foam. The presence of a foam bank between the displacing air and the displaced brine improved both breakthrough and ultimate recovery. In a continuation of this work Kolb5 attributed the great reduction in surfactant solution production rate as displacement by air progressed to a decrease in relative permeability to gas. These several effects reported by both Bennett and Kolb can all be attributed to the high apparent viscosity of foam which was obviously flowing in the porous media. Bernard and Holm6 observed that gas permeability in the presence of foam was much less than in the absence of foam when both were measured at the same gas saturations. This reduction in permeability was proportionately greater for more permeable sands and sandstones than for the less permeable. Viscosities of the surfactant solution and air were used for their calculations rather than apparent foam viscosities. The foam stability increased with decreasing permeability, which was probably due to the natural regenerative properties of porous media, and the permeability reduction increased

Jan 1, 1967

By D. L. Archer, W. W. Owens

Conventional waterflooding often is uneconomic in highly fractured reservoirs because of the gross bypassing of the reservoir oil by injected water. Imbibition and pressure pulse flooding have been used in at least one fractured reservoir in an attempt to achieve better oil recovery performance. This paper presents the results of laboratory flow tests conducted on large cores to evaluate the possible applicability of these methods (particularly pressure pulse flooding) to different types of reservoir systems. Test data were obtained on both water-wet and oil-wet systems, and on systems having two widely different levels of compressibility and flow capacity. Results indicate that pressure pulse recoveries from fractured reservoirs will likely not exceed 5 to 10 per cent pore space with maximum response achieved during the first pulse cycle. Improved recovery by this method is possible from both oil-wet and water-wet reservoirs. Comparable saturation distributions during imbibition and pressure pulse production suggest that an initial pressure pulse cycle to speed production response would not interfere with subsequent imbibition flooding in water-wet reservoirs. INTRODUCTION There is a steady increase in the number of fluid injection projects being initiated each year to improve oil and gas recovery over that obtained by primary production mechanisms. Experience gained in different geographical areas and with different recovery methods is teaching that reservoir anatomy is one of the more important factors controlling the success or failure of such projects. Fractured formations appear to be the rule rather than the exception, especially in carbonate reservoirs. Conventional methods of fluid injection, whether it is waterflooding, gas injection or miscible flooding, have limited applicability to highly fractured reservoirs because of the severe bypassing of reservoir fluids. The result is early breakthrough of the injected fluid and rapidly increasing ratios of injected to in-place fluid with an undesirable effect on return on investment. Thus, innovations to conventional methods must be developed if recovery from fractured reservoirs is to be optimized. One new method proposed for application to highly fractured reservoirs is pressure pulse waterflooding. Although pressure pulsing with gas has been suggested' and rested in at least one reservoir; the Sohio Oil Co. is generally credited with the development and testing of pressure pulsing in conjunction with waterflooding. Publications" of Sohio's waterflood operations in the Spra-berry Trend indicate that the idea for pressure pulse flooding evolved from a critical analysis of the disappointing early performance of their flood initiated in a portion of the trend in April, 1961. This flood was planned to take advtantage of imbibiltion and high pressure gradients across the reservoir matrix blocks (unfractured Mocks of the reservoir rock which are visualized to be generally surrounded 'by the fracture system) evaluated in pilot tests by the Atlantic Richfield Co.5 and Humble Oil & Refining Co.' However, the early recovery and pressure performance of this flood indicalted that high pressure gradients induced between the fracture system and the centers of the matrix blocks were forcing water into the periphery of the blocks and temporarily interfering with the counter-current flow of oil due to imbibition. Fill-up in the blocks was occurring with re-solution of the free gas phase as pressures climbed above the original reservoir bubble point. It appeared that cessation of water injection to permit the capillary forces to become dominant and that expansion of the rock and its contained fluids during pressure reduction might aid in expulsion of oil from the rock matrix into the fractures. Subsequent production performance. aflter cessation of injection, proved this hypotheslis correct. Thus, pressure pulse waterflooding was born. Pressure pulse flooding appears to have several advan'tages over imbibition type flooding in highly fractured reservoirs. First, all wells may be used for water injection which should hasten fill-up and achieve a more rapid increase in reservoir pressure. Second, because of increased injection pressures and rates, flooding gradients will force water into the reservoir matrix more rapidly than would be achieved by imbibition alone. Third, during the pressure depletion cycle of the flood, the compressibility of the system, which has been increased by resolution of the free gas phase, provides energy for the displacement of oil at a rate greater than would correspond to countercurrent flow during imbibition. Fourth, during depletion all wells may be put on production, thus contributing to higher oil withdrawal rates on a reservoir-wide basis. One possible disadvantage of pressure pulse flooding as compared to imbibition floodling, however, is that the outer periphery of each matrix block is flooded to a residual or near residual oil saturation during the Water lnjection stage. This zone of reduced permeability to oil m'ay offer greater restriction to the production of

Jan 1, 1967

By B. A. Eaton, R. H. Jacoby

This paper presents the depletion-performance correlations developed using data from PVT studies and well-completion tests of 27 rich gas-condensate and volatile-oil reservoir fluids. The PVT behavior for each fluid war determined experimentally in the laboratory, and the depletion performance of each fluid was then calculated. Well-completion test data consisting of reservoir pressure and temperature, initial total gas-oil ratio, and stock-tank oil gravity were used with the depletion-performance data to develop general performance correlations. Three correlations are presented. The first equation shows the empirical relation of gas in place per barrel of hydrocarbon pore space with reservoir pressure and temperature, and initial total gas-oil ratio. The second equation is the correlation of stock-tank oil in place per barrel of hydrocarbon pore space. The third equation relates stock-tank oil production (to 500-psia depletion pressure) with pressure, temperature and oil gravity. All three correlations are also presented in convenient graphical form. Gas production to 500-psia depletion pressure is approximately 92.6 per cent of the gas in place for each fluid, so no correlation was developed for gas production. INTRODUCTION The rate of discovery of shallow, low-pressure and low-temperature oil reservoirs in the United States has been declining in recent years. For this reason, the search for oil has gone to greater depths where high pressures and temperatures are encountered. A high-pressure and high-temperature reservoir is the environment in which many rich gas condensates and highly volatile oils are found. Because of this trend, the performance prediction of such reservoirs has become an important task of the reservoir engineer. The need and usefulness of a performance correlation based on field data that are available from the usual well-completion tests have previously been established. These earlier attempts have shown that the phase behavior of rich gas-condensate systems and volatile-oil systems can be correlated with fluid properties, pressure and temperature. However, very little information on gas-condensate depletion performance5 was available at the time of these earlier attempts. The availability of more recent data has made possible the development of an improved correlation. The data, as listed in Tables I and 2, are for 25 gas-condensate and two volatile-oil systems that were studied. The problem was to correlate the data for in-place fluids and the production data with the initial total gas-oil ratio, initial reservoir pressure and temperature, and stock-tank oil gravity. These correlations were obtained by multiple regression analysis, whereby the data were fitted to an equation by the method of least squared error. BASIS OF THE CORRELATIONS 1. The data used in the development of these correlations were obtained from depletion-performance calculations6 and well-test or production measurements. 2. All of the independent variables used in the correlations are available from the usual well-test completion data. 3. The reservoir volume considered is 1 bbl of hydrocarbon pore space (HCPS). 4. The gas-oil ratio used is the total standard cubic feet of gas (at 14.7 psia, 60F) coming off the entire separator train per barrel of stock-tank oil. This minimizes the effect of separator conditions which might vary from system to system. 5. The correlated pressures are the actual reservoir pressures which might or might not be the saturation pressure of the reservoir fluid.

Jan 1, 1966

By H. O. Jahns

This paper describes the application of regression analysis for obtaining a two-dimensional areal description of heterogeneous reservoirs from short-term pressure-time data such as that obtained in interference tests. The method replaces the time-consuming trial-and-error procedure commonly used to match field data on an electric analyzer or digital computer with a systematic search which is programmed for a computer. The computer program adjusts the properties of a reservoir model automatically until a least-squares fit is obtained between observed and calculated pressure data. The reservoir is simulated by a single-phase, compressible, two-dimensional model. It is divided into a number of homogeneous blocks whose transmissibility (kh/F) and storage (Fch) values are varied to obtain the least-squares fit. The reliability of these values is determined from their standard deviations and correlation coefficients. Although the method is rigorously applicable to single-phase flow only, multiphase flow can be handled provided saturation changes are small during the test. Possibly the method can also be used to obtain a reservoir description from pressure-production history, but this application is outside the scope of this work. The paper includes, in addition to a description of the numerical procedure, a discussion of some of the problems associated with the method. Rules are given to help in selecting the number of homogeneous blocks and deciding upon their arrangement. The uniqueness of a reservoir description is considered. Finally, the use of the method is illustrated by the interpretation of field data from two interference tests. INTRODUCTION Pressure data from short-term transient tests, such as single-well and interference tests, are widely used to obtain reservoir properties. These tests are usually analyzed by assuming a simple reservoir model; very often, a homogeneous one is used. As a result, analysis of the transient data from each well frequently gives different values for reservoir properties. The problem then arises to combine all these differing results into a more detailed picture of the reservoir. One technique is to simulate the reservoir with a digital computer or with an electrical analyzer and to adjust the reservoir parameters by trial and error until the simulated pressure data are in reasonable agreement with the observed pressure data for all wells. Although this method has been used for both transient tests and pressure-history data, it is time-consuming and subjective. A second technique uses regression analysis to replace the trial-and-error procedure with a systematic search that can be programmed for a digital computer. Use of regression analysis in reservoir description was proposed recently by Jacquard and Jain.1 They divided the reservoir into a number of homogeneous blocks whose properties are varied until a least-squares fit is obtained between observed and calculated pressures. However, they did not consider their technique to be operational, mainly because of ".. .(I) the lack of experience in using the method .. . notably for the improvement of convergence; and (2) limitations imposed by the insufficiency of available computers". 1 While the analysis presented in this paper applies the same general principle used by Jacquard and Jain, the specific method is significantly different. Some differences are (1) the regression problem is solved in a different way which requires less computer time in most cases; (2) a stepwise solution, in which the detail in the reservoir description is increased from step to step, is used to improve convergence; and (3) the reliability of the estimated reservoir .properties.. as measured by their standard deviation and correlation coefficient, is estimated'

Jan 1, 1967

By G. W. Thomas

A mathematical tnodel of forward combustion in an oil reservoir is treated in this paper. The model describes a radial system having a vertical section of essentially infinite thickness, all of which is permeable to gas flow. Combustion, however, is presumed initiated over a limited thickness of the total vertical section. In the interval supporting cotnbustion, the mechanisms of radial conduction, convection and heat generation are taken into account. Above and below the burning interval, heat transport in the radial direction is by cottduction and convection. Vertical heat losses from the ignited interval are accounted for by conduction alone. A general solution is presented for the temperature distribution caused by radial movement of the combustion front. The results show that no feedback of heat occurs into the ignited interval when convection and conduction are acting in the bounding media. Peak temperatures are also 5 to 10 per cent higher than in the case where heat transport in the bounding media is by conduction alone. We arbitrarily define vertical coverage to be that fraction of the total ignited interval which is at 600F above atnbient, or greater, at any given time. The radial distance at which the vertical coverage becomes zero is the propagation range of the combustion front. It was found that an increase in vertical coverage results when the oxygen concentration, fuel concentration or gas-injection rate is increased. Moreover, the combustion front can be propagated 10 to 15 per cent further than in the case where only conduction is acting above and below the ignited interval. INTRODUCTION In the theoretical treatment of forward combustion in a radial system, one of the problems encountered is the determination of the transient temperature distributions caused by an expanding cylindrical heat source. Bailey and Larkin' and Ramey' simultaneously presented analytical solutions to the problem assuming heat transport by conduction alone. In a subsequent publication, Bailey and Larkin3 included the effects of both conduction and convection while treating linear and radial models. In this latter work, however, vertical heat losses were largely neglected. Selig and Couch' dealt with a radial model in which both conduction and convection were acting. Only a limiting case involving vertical heat losses was considered, however. Namely, temperatures on the boundary of the bed of interest were set equal to zero. Solutions thus obtained were representative of a system having a maximum vertical heat flux. Chu5 recently treated a more general case in which a permeable bed was considered bounded by impermeable media. Conduction and convection took place within the bed, and only conduction outside of the bed. The effects of vertical heat losses were included in his study. Solutions were obtained by numerical techniques. This paper is an extension of the theoretical work of other authors pertaining to forward combustion in a radial system. In particular, a mathematical model of the process is treated in which heat generation occurs over a small vertical interval of a larger permeable section. In the interval supporting heat generation, and above and below this interval, the mechanisms of radial conduction and convection are also presumed acting. Heat losses from the ignited interval are accounted for by vertical conduction. An analytical solution for the temperature distribution caused by radial movement of the burning front is presented. The effects of certain process variables are indicated and comparisons with Chu's results are made. THEORY To render the mechanism of forward combustion tractable to mathematical treatment, we idealize the problem to the extent of assuming continuous reservoir media possessing homogeneous and isotropic properties. The following additional assumptions are implicit in this analysis. 1. The thermal parameters, i.e., heat capacities, thermal conductivities and thermal diffusivities are invariant with temperature and pressure. Moreover, the bounding media possess the same thermal properties as the bed of interest. 2. The temperatures of the porous media and its contained fluids at any point and at any time are equal. 3. The reaction rate between the oxidant gas and the fuel is infinite. This assumption implies that the incoming oxygen concentration instantaneously goes to zero within an infinitesimal distance, i.e., the width of the combustion zone is negligible. 4. The rate of gas injection is constant and corresponds to the average rate throughout the lifetime of the project. 5. The fuel concentration is constant throughout the volume of rock swept out by the burning zone. 6. There is complete burnoff of fuel. This assumption demands that the rate of propagation of the burning front equals the rate of fuel burnoff. In a radial system, with a

By R. L. Perrine

This paper presents results obtained by consider ation of unstable miscible displacement within a porous medium as a quasi-turbulent displacement process, that is, a process characterized by fluctuations around a simpler, mean flow behavior. The instability usually results from the miscible displacing fluid having a viscosity lower than that of the displaced fluid. Such instability is evidenced in suitable experiments by the presence of viscous fingering. A unified theoretical relationship applicable to both stable and unstable displacement is derived. Parameters in the theory are obtained by empirical means. Although the theory describes an idealized conceptual model, agreement with experiment is quite reasonable up to moderately adverse viscosity ratios, such as 17:1. INTRODUCTION Many experimental studies of miscible displacement in porous media have been conducted with prediction of reservoir displacement efficiency as the ultimate objective. Most such studies have utilized lower displacing than displaced fluid viscosity, scaled to potential reservoir fluid pairs. Theoretical approaches have been largely limited to unit viscosity ratios, however, in spite of the necessity for an understanding of the mechanism of the unstable, adverse viscosity ratio displacement process. An obvious reason is the difficulty of describing in mathematical form the viscous fingering characteristic of these conditions. Observation of experiments conducted with dyed fluids in transparent systems suggests an analogy between unstable miscible displacement and turbulent flow in a pipe. In both cases there are fluctuations around a simpler, mean flow behavior. An important difference is that flow behavior of interest in the porous medium is entirely transient, contained within a transition zone between displacing and displaced fluids. Transient behavior complicates description in that coefficients in the equations become variables rather than constants. In the study reported here, the analogy with turbulent flow has been used in treating unstable miscible displacement as a quasi-turbulent displacement process. The purpose has been to derive, even if restricted to an idealized conceptual model, a unified theoretical relationship applicable to both stable and unstable displacement. A relationship meeting these specifications up to moderately adverse viscosity ratios, such as 17:1, is presented. One fluctuation parameter in the theory and dispersion coefficients are obtained by empirical means. The idealized theory is compared with experimental results. UNSTABLE MISCIBLE DISPLACEMENT AS QUASI-TURBULENT DARCY FLOW The miscible displacement of one fluid by another within a porous medium is considered in this study. Flow conditions are such that Darcy's law is applicable. It is further assumed that, by the stability criterion proposed by perrine,1 initial flow conditions lie well within the regime of instability. Thus substantial viscous fingering is assured. We wish to show how this flow regime can be represented as quasi-turbulent. That is, the Reynolds number for the established flow conditions is below that for inertial or turbulent flow within a porous medium, and lies within the regime for which Darcy's law is valid. Yet flow behavior can be described as the combination of some relatively simple average result, and characteristic fluctuations that are superimposed on the simpler behavior. Stated another way, flow behavior includes the movement of layers of fluids with different velocities past or over one another. Such descriptions are characteristic of turbulent pipe flow. The fluctuations in the present case are a direct consequence of the local viscous fingering which accompanies the unstable displacement process. Should displacement become stable, fluctuations would die out. It is of particular importance to note that fluctuations such as these can interact in a way that contributes to material transport by the basic flow. The basic transport equations required for engineering calculations must be modified to reflect this fact.

By W. B. Braden

This paper presents a suitable method for predicting gas-free oil viscosities at temperatures up to 500F knowing only the API gravity of the oil at 60F and the viscosity of the oil measured at any relatively low temperature. The API pravity and the one viscosity value are used as parameters to determine the slope of a straight line on the ASTM slanaord viscosity-temperature chart. Then, knowing the slope of the line and one point on the line, the vrscosities at higher temperatures can be determined. The line slope correlations were developed at I00 and 210F since viscosity data are frequently measured at these temperatures. A procedure is given for predicting line slopes from measurements at other tetnperatures. A nomogram is furnished for solving the relationship. The correlation has been evaluated at temperatures up to 5OOF for oils varyzng in gravity from 10 to 33 " API. The correiution is applicable only to Newtonian fluids. Comparison at 500F of true viscosities and those predicted from values at 100F shows an average deviation of 3.0 per cent (maximum deviation of 6.0 per cent). Predictions from the values at 21 0F for the same oils how an average deviation of 1.5 per cent (maximum deviation of 3.4 per cent). INTRODUCTION Correlations have been developed by Beal' and by Chew and Connally' for predicting viscosities of gas-saturated oils at reservoir conditions. Each of these correlations requires a knowledge of the solution gas-oil ratio and the viscosity of the gas-free oil at the reservoir temperature. For temperatures below 350F, measurements of the gas-free oil viscosities can be made easily using commercially available equipment. In thermal recovery processes, however, reservoir temperatures well in excess of 350F are encountered. Viscosity measurements at such conditions are more difficult and time consuming and require modification of existing equipment or the construction of new equipment. Measurements are further complicated by the difficulty of handling highly viscous oils associated with thermal recovery processes. Therefore, it is desirable to have a correlation which allows accurate prediction of viscosities at high temperatures. A commonly used technique for predicting viscosities at high temperatures is to measure the viscosities at two lower temperatures, plot the values on ASTM standard viscosity-temperature charts and extrapolate to the temperatures desired. If either of the values is slightly in error, the extrapolated value can be significantly in error. To justify an extrapolation, three points are actually necessary. This procedure can consume much time, particularly with heavy oils. Considering the cost of viscosity measurements, it would be desirable to eliminate the need for direct measurements by having correlations which would allow viscosity predictions from other physical or chemical properties. Beal1 investigated the possibility of correlating viscosity with oil gravity at temperatures from 100 to 220F. While showing that a general relationship exists, he also found significant deviations. It is possible that correlations will be developed based on oil composition as more information becomes available. While not eliminating the need for viscosity rneasurements, the method presented herein requires that only one viscosity measurement be made. The API gravity must also be known. The theory is based on the fact that the viscosity of paraffins (high gravity) changes less with temperature than does the viscosity of naph-thenes or aromatics (low gravity). The gravity. therefore, is used as a parameter to determine the slope of a straight line on the ASTM standard viscosity-temperature charts. The correlation is applicable only to Newtonian oils, and deviations due to thermal decomposition and nonhomo-geneity cannot be predicted. Oils containing additives have not been evaluated. PROCEDURE Fifteen oils were used in developing the correlation; eight were crudes and seven were processed oils. Oil gravities varied from 9.9" API (naphthene base) to 32.7' API (paraffin base). The temperature range studied was 81 to 516F. Each oil used had a minimum of three viscosity measurements and each plotted essentially as a straight line on the ASTM charts. In all, 91 viscosity measurements were used in the correlation. Saybolt, rolling ball and capillary tube viscometers were used for the measurements. Viscosity data for Samples 1, 2, 4, 7, 10, 11 and 14 were obtained in Texaco, Inc. laboratories. The data for Samples 3, 5, 6, 8, 9, 12 and 15 were from Fortsch and Wilson,3 and data for Sample 13 were from Dean and Lane.' All data points used in the correlation are plotted in Fig. 1. It is seen that some of the viscosity values deviated slightly from the straight-line plots at the higher temperatures. Properties of the oils after exposure to the

Jan 1, 1967

By L. T. Stanley

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

By W. E. Showalter

This paper discusses some of the results of combstion-drive tests which were made in a test cell using a sand bed 10 in. in diameter x 10-ft long. The test method is illustrated and described. The relationship between the API gravity of the in situ oil and the amount of air required for combustion drive is discussed in detail. Other things constant, the air requirement for combustion drive increases as the APl gravity of the in situ oil decreases. If the test results apply to actual reservoirs, the lowest-priced oils may cost the most to recover by this method. Information is shown which indicates that the effect of pressure on the amount of hydrocarbon burned is not large. A method of predicting air requirements from the API gravity of the in situ oil is presented. INTRODUCTION Combustion drive is the term used to identify the process of interstitial or in situ burning as an oil recovery method. Part of the in situ oil is burned to generate the energy needed to produce the remainder of the oil. Combustion drive as an oil recovery mechanism remains an economic uncertainty in spite of all the work that has been done by the industry in both laboratory and field?-6 This paper will show some of the results of tests which were made in a test cell for the purpose of studying the nature of the combustion-drive process. It will present data which indicate that the API gravity of the in situ oil is a significant indicator of the amount of air required to drive a burning front through oil sand. Air requirement varies inversely as the API gravity of the in situ oil. EXPERIMENTAL The tests were performed in a cell which utilized a cylindrical sand section 10 in. in diameter x 10-ft long. The thin-walled metal pipe which held the sand was wound with twenty external electrical resistance heaters which, by means of an automatic controller, maintained adjacent sections of the wall of the pipe at temperatures equal to the temperatures of the contained sand. Each heater covered 6 linear in. of the pipe. By this means lateral heat loss from the sand section was minimized, thereby causing the sand section to simulate more closely a horizontal increment of a combustion-drive reservoir. Fig. 1 shows a schematic diagram of the test assembly. Thermocouples to measure the temperature in the sand were located every 6 in. along the length of the sand section. The pipe containing the sand was enclosed in a cell designed for an operating pressure of 500 psig. The inlet air pressure was controlled at the inlet, and the gas flow rate was controlled and measured at the outlet of the cell. The oil sand used for the tests was prepared by mixing first water and then oil with the noncon-solidated sand using a closed mixer similar to a cement mixer. Table 1 shows a screen analysis of the sand. Ninety percent of the sand was 100 mesh or finer. This sand was a mixture of 80 per cent No. 120 Nevada White Sand and 20 per cent Tennessee Hi-Fusion Moulding Sand No. 3. The Nevada sand was a clean silica sand. The Tennessee sand con-

By K. H. Coats, M. H. Terhune

Analysis and example applications have been performed to compare the accuracy and computing speed of alternating-direction explicit and implicit procedures (ADEP and ADIP) in numerical solution of reservoir fluid flow problems. ADIP yields significantly greater accuracy and requires about 60 per cent more computing time than ADEP, not 300 or 500 per cent more as reported elsewhere. l,2 INTRODUCTION Several recent papers 1-3 discuss an alternating-direction explicit difference approximation (ADEP) to the diffusion equation. Example applications of ADEP and ADIP4 were reported to support conclusions that ADEP is comparable in accuracy to ADIP and requires one-fifth to one-third the computing time of ADIP. Applications of ADEP in calculation of two-phase flow in reservoirs was also proposed. 3 This study was performed to compare further the relative merits of ADEP and ADIP in simulation of two-dimensional flow of one and two fluid phases in reservoirs. Since two-phase flow equations are often essentially elliptic rather than parabolic, the efficiency of ADEP in solving the elliptic equation was also examined.

Jan 1, 1967

By D. G. Russell

A simple method has been developed with which flowing bottom-hole pressure data from two-rate flow tests in oil or gas wells can be analyzed to estimate the formation permeability, skin factor and average reservoir pressure. The required pressure data are obtained by observation of the transient bottom-hole pressure behavior after the stabilized producing rate of the well is changed to another, higher or lower, rate. The new method yields the same information as a conventional pressure buildup analysis, but eliminates the need for closing in the well. The analysis of a two-rate flow test is of the same degree of difficulty and requires about the same engineering time for application as a conventional pressure buildup analysrs. The extended closed-in periods experienced with conventional buildups because of long, low-rate afterproduc-tion periods are eliminated by fIow tests. Other anomalous pressure buildup effects, such as "humping" due to weli-bore phase segregation, can be successfully eliminated with the new method. Generally, flow tests of about 24 hours' duration run with conventional pressure measurement equipment are sufficient for interpretation purposes. Field testing of the two-rate flow test method has established that it is a reliable and economical method which can be used in many instances to complement or even replace conventional pressure buildup methods. INTRODUCTION The principal method for estimating formation permeability and well damage, or skin factor, in a producing oil or gas well is the analysis of shut-in bottom-hole pressure buildup data.' This familiar method has been used quite successfully by reservoir engineers for many years. It is based on the solution of the radial flow equation for constant rate conditions, and requires that the well be closed in for a sufficient period of time to obtain a clearly defined linear portion on the plot of observed t + ?t bottomhole pressure vs log t + ?/?t(where At is shut-in A? time, and t is producing time to the instant of shut-in). From the slope of the plot and other normally obtainable data, the permeability, skin factor, and reservoir pressure at infinite shut-in time (if the reservoir were infinite) can be estimated. Over the years several drawbacks have become apparent in the use of conventional shut-in pressure buildups for determining permeability and skin factor. The conventional pressure buildup interpretation theory assumes that a well is closed in at the sand face and that no production into the well occurs after shut-in. In practice, of course, the well is closed in at the surface, and inflow into the well continues until the well fills sufficiently to transmit the effect of closing in to the formation. This adjustment period is commonly referred to as the "after-production" portion of the pressure buildup. In the tight reservoirs, long, low-rate after-production periods frequently occur, and the well must be shut in for several days or, in some instances, even weeks to obtain an interpretable buildup curve.' Obviously, such long shut-in times can cause loss in current income, both from reduced oil production and from the fact that personnel and pressure measurement equipment are occupied with a single well for too long a time. In other cases, even long shut-in periods do not seem to be of much aid in obtaining an interpretable buildup." If there is considerable phase redistribution (liquid fallout or bubble rise) after a well is shut in, then curves with no interpretable portion are often obtained during the buildup. In addition to instances in which wellbore effects cause trouble, there are also cases in which the major objection to use of the closed-in pressure buildup is simply the fact that the well must be shut in. When there is no proration and when the well has limited producing capacity. closing in the well means loss of income. From the foregoing discussion it is apparent that it ib desirable to have an alternative method of obtaining the same information as that derived from a conventional buildup without the need of closing in the well. One possible solution which has been offered for this problem is the use of a bottom-hole shut-in tool' which isolates the major portion of the flow string from the formation face during the buildup. In this paper an alternative method which is frequently successful in avoiding wellbore effects and which does not require the use of special equipment is presented. A new, simple method has been developed with which the flowing bottom-hole pressure data from flow tests in oil or gas wells can be used to estimate permeability, skin factor and the average reservoir pressure. The required pressure data are obtained by observation of the transient bottom-hole pressure behavior after the stabilized production rate of the well is changed to another, higher or lower rate. The need for closing in the well is eliminated, and pressure measurement periods of only 18 to 24 hours arc usually sufficient, even in tight reservoirs. Thus, the new

By J. K. Elliott, P. H. Kelly

Many engineering functions such as surface metering work and laboratory compressibility check points involve the use of liquid densities of light hydrocarbon mixtures at various pressures and temperatures. However, at the present time, no simple reliable method exists for determining density variation, particularly if the composition of the liquid is unknown. Consequently, a study was undertaken to develop and present a simple and accurate method of predicting density variation of a light hydrocarbon liquid with pressure and temperature, knowing only the density of the liquid at some condition. The experimental liquid compressibility data from API Project 37 by Sage and Lacey' have been considered to be accurate within 0.5 per cent and cover a wide range of pressure (14.7 to 10,000 psia), temperature (100" to 400°F) and molecular weight (up to 150). From these data, a set of liquid density curves, which relate density to pressure, temperature and molecular weight, was developed. These curves make it possible to predict density variation with pressure and temperature. Compared to extensive laboratory compressibility data on a complex, light hydrocarbon liquid, the use of the liquid density curves resulted in an average error of less than 0.5 per cent. Based on the results of this analysis, it is concluded that the set of liquid density curves developed from the data of Sage and Lacey provides an accurate and simple method for predicting the density variation of light hydrocarbon liquids when the density at some condition is known. These curves should be very helpful in many engineering calculations, particularly in the surface metering of light hydrocarbon liquids. INTRODUCTION Many situations arise in field and engineering laboratory work, such as reservoir engineering studies, check of experimentally determined laboratory data and orifice flow-meter formulas, where liquid density factors at various pressure-temperature conditions are required. Also, the need for accurate light hydrocarbon liquid information has become more important with the advent of miscible-type displacements for secondary recovery purposes in oilfield operations. Several reliable methods are available1 - "or determining the density of liquid hydrocarbons if the composition of the liquid is known. However, there is a definite lack of methods for accurately determining the variation of density when the composition of the liquid is unknown. The purpose of this study is to review the various methods for determining hydrocarbon liquid densities and to develop a simple and reliable method of determining variation in density of light hydrocarbon liquids with pressure and temperature when the compositio~n of the liquid is unknown. METHODS FOR DETERMINING DENSITY OF LIQUIDS OF KNOWN COMPOSITION Sage, Lacey and Hicks' have proposed a method to predict the density of light liquid hydrocarbons by using partial molal volumes. Data are available on experimentally developed partial liquid volumes of hydrocarbons over a rather limited range of temperature, pressure and composition. The partial mold volume method has proved satisfactory for determining the density of some hydrocarbon liquids when the composition is known. Within the range covered in the Sage, Lacey and Hicks1 data, the results agree within about 3 per cent of the experimental values. Hanson mentions the limitation of this method to a composition range of approximately 10 per cent by weight of methane, which will not allow this correction to cover most low molecular weight-light hydrocarbon liquids. Standing and Katz2 studied data on light hydrocarbon-liquid systems containing methane and ethane at high temperature and pressure and have presented a method for determining liquid densities, assuming additive volumes for all components less volatile than ethane and using apparent densities for methane and ethane. The compressibility and thermal-expansion curves used by Standing are based on assumptions that compressibility of a hydrocarbon liquid at temperatures below 300°F is a function of the liquid density at 60°F and that thermal expansion of the liquid is affected little by pressure. The information required to use this technique with an example problem is furnished by Standing.' Hanson eports an average error of - 0.5 per cent using the method of apparent densities in calculating liquid densities of several volatile hydrocarbon mixtures. However, as implied, the apparent density method is not applicable for liquid density calculations when the composition of the liquid is unknown. Watson- as presented a method

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

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

By M. R. Tek, R. L. Nielsen

The scaling laws as formulated by Rapoport relate dynamically similar flow systems in porous media each involving two immiscible, incompressible fluids. A two-dimensional numerical technique for solving the differential equations describing systems of this type has been employed to assess the practical value of the scaling laws in light of the virtually unscalable nature of relative permeability and capillary pressure curves and boundary conditions. Two hypothetical systems — a gas reservoir subject to water drive and the laboratory scaled model of that reservoir — were investigated with emphasis placed on water coning near a production well. Comparison of the computed behavior of these par ticular systems shows that water coning in the reservoir would be more severe than one would expect from an experimental study of a laboratory model scaled within practical limits to the reservoir system. This paper also presents modifications of the scaling laws which are available for systems that can be described adequately in two-dimensional Cartesian coordinates. INTRODUCTION Present day digital computing equipment and methods of numerical analysis allow realistic and quantitative studies to be carried out for many two-phase flow systems in porous media. Before these tools became available the anticipated behavior of systems of this type could be inferred only from analytical solutions of simplified mathematical models or from experimental studies performed on laboratory models. To reproduce the behavior of a reservoir system on the laboratory scale, certain relationships must be satisfied between physical and geometric properties of the reservoir and laboratory systems. Where the reservoir fluids may be considered as two immiscible and incompressible phases, the necessary relationships have been formulated by Rapoportl and others.2-5 Rapoport's scaling laws follow from inspectional analysis of the differential equation describing phase saturation distribution in such systems. It will be recalled that these scaling laws presuppose three conditions: (1) the relative permeability curves must be identical for the model and prototype; (2) the capillary pressure curve (function of phase saturation) for the model must be linearly related to that of the prototype; and (3) boundary conditions imposed on the model must duplicate those existing at the boundaries of the prototype. These three requirements seldom if ever can be satisfied in scaling an actual reservoir to the laboratory system because: (1) The laboratory medium normally will be unconsolidated (glass beads or sand) while the reservoir usually is consolidated. Relative permeability and capillary pressure curves are usually quite different for consolidated and unconsolidated porous media. (2) The reservoir usually will be surrounded by a large aquifer which could be simulated in the laboratory only to a limited extent. (3) Wells present in the reservoir would scale to microscopic dimensions in the laboratory if geometric similarity is to be maintained. In view of these considerations, rigorous scaling of even a totally defined reservoir probably would never be possible. The purpose of this paper is to assess the practical value of the scaling laws in the light of the unscalable variables. This has been done by carrying out numerical solutions in two dimensions to the differential equations describing the flow of two immiscible, incompressible fluids in porous media for a field scale reservoir and a laboratory model of that reservoir. While both the reservoir and the laboratory model were purely fictional, each has been made as realistic and representative as possible. The field problem selected as the basis for the investigation was an inhomogeneous, layered gas reservoir initially at capillary gravitational equilibrium and subsequently produced in the presence of water drive. The laboratory model of this reservoir was designed to utilize oil and water in a glass bead pack. The numerical treatment employed was similar to that of Douglas, Peaceman and Rachford6 and it included both capillary and gravitational forces as

By D. G. Russell

Two techniques have been developed with which the applicability of pressure build-up analyses can be extended to include pressure data which previously have been considered virtually unusable. One of the interpretation methods makes possible the analysis of pressure build-up performance during the wellbore fill-up or after production period which occurs soon after a well is closed in. The other technique is an extension of a method for analyzing pressure build-up performance during the late-time portion of the pressure build-up which occurs after boundary eflects first begin to alter the shape of a conventional pressure build-up curve. With both of these methods it is possible to obtain estitnates of the kh product, the skin factor and the reservoir pressure. In addition, with the late-time analysis technique it is possible to obtain an estimate of the contributory drainage volume of the well being tested. This means that in some cases a check on reservoir limit test and (or) material-balance calculations can now be obtained from pressure build-ups. Both methods are slightly more time-consuming than conventional pressure build-up analysis methods because trial-and-error plots of pressure data must be made. The late-time method for analysis of pressure build-ups is in principle applicable to the late-time portion of a two-rate flow test or a pressure drawdown test. The interpretation formulas and procedures for these types of tests are also outlined. In these cases, as with pressure build-ups, it is significant that an estimate of the contributory pore volume is also obtained. On the basis of limited experience with the new techniques, it appears that satisfactory estimates of the kh product, skin factor, reservoir pressure and, for late-time analysis, contributory drainage volume can be obtained. INTRODUCTION The analysis of bottom-hole pressure build-up behavior in closed-in wells has been a a subject of interest in petroleum engineering circles for many years. In fact, few other subjects have received as much attention as pressure buildup analysis methods have. The cause for this interest is essentially twofold in nature. First, the pressure behavior of a well can normally be measured with a reasonably high degree of accuracy so that good data for analysis can be obtained. Secondly, over a fairly wide range of operating conditions, valuable information as to the quality of the reservoir rock and completion efficiency of the well can be obtained at a nominal cost. In recent years, numerous papers have been prepared on the effects of various operating conditions and reservoir heterogeneities on pressure buildup behavior. Very little work has been done, however, on extension of pressure build-up analysis methods to those pressure data which are not amenable to analysis by the present methods. The theory upon which the analysis of shut-in bottom-hole pressure build-up data is based is derived from the solution of the radial flow equation for a slightly compressible fluid for constant-rate conditions. It requires that the well be closed in for a sufficient period of time to obtain a clearly defined linear portion on the plot of observed bottom-hole pressure vs log (t + ?t)/ ?t (where At is shut-in time, and t is producing time to the instant of shut-in). From the slope of the plot and other normally obtainable data, the formation permeability, the well damage or skin factor, and the reservoir pressure at infinite shut-in time (if the reservoir were infinite) can be estimated. The successful application of this procedure depends on being able to recognize the straight-line section on the basic pressure build-up plot. The presently used pressure build-up interpretation theory also assumes that a well is closed in at the sand face and that no production into the well occurs after shut-in. In practice, of course, the well is closed in at the surface, and inflow into the well continues until the well fills sufficiently to transmit the effect of closing-in to the formation. This adjustment period is commonly referred to as the "afterproduction" or "fill-up" portion of the pressure build-up. During the period that the well fill-up effect is most pronounced, the basic pressure build-up plot is nonlinear. At later shut-in times after the effects of a drainage boundary have been felt at the well, deviation from the straight-line behavior of the pressure build-up plot also results. In many cases either of these effects or a combination of both can make the straight-line portion on the pressure build-up plot difficult to recognize. Obviously, an extension of pressure build-up analysis methods to include the afterproduction period and the period in which boundary effects are being felt would be desirable and might render valuable pressure data which for years have been considered virtually unusable. The principal reference of note concerning pressure buildup analysis during the afterproduction period is a paper by Gladfelter, Tracy and Wilsey.1 In the approach of these authors it is necessary to measure the rate of influx into the well during the afterproduction period. This is done through sonic measurements or through measurement of tubing-head and casing-head pressures simultaneously with

Jan 1, 1967

By C. R. Knopp

Gas-oil relative permeability ratio is an important relationship in oil reservoir predictive calculations. A correlation has been developed from 107 gas-flood k/k tests on Venezuelan core samples. The correlating parameter is based on restored-state water saturation tests and' is applicable to both consolidated and poorly consolidated sandstone reservoirs. Data of the correlation show that there are no distinguishah1e differences between the mass-data groupings for the two c1assifications A procedure is recommended for running .sufficient relative. permeability analyses to compute a geometric mean of the sample group. The geometric mean is more representative of the total core, and probably the entire reservoir. For example, while only one in four of the k,,/k,., test curves agreed closely with the resultant correlation of this report, the geometric mean curves of the 16 suites (three samples or more). showed good agreetment ill three cases out of four. INTRODUCTION The gas-oil relative permeability ratio is an important, fundamental relationship in most oil reservoir predictive calculations. Predictive calculations are made to estimate future reservoir production characteristics and ultimate oil recovery. The k1,/k2, relationship is specifically needed to relate the surface gas-oil ratio to the reservoir oil and gas saturation, and to calculate the relative movement of these phases within the reservoir whenever some of the more complex driving mechanisms are present. Laboratory k1/k2, tests are not generally run as a routine analysis. Consequently, k1/k2 data often are not available when needed because the cost of laboratory work could not be justified or the need for such data had not been properly anticipated. When laboratory k1/k2, data are available, they are often very difficult to interpret. For example, wide divergence is sometimes shown in a family of k1,/k1, tests representative of the producing horizon in a single well. With these considerations in mind, a study was made to determine if a relationship might exist between the k1,/k2, curve and some other simple laboratory test criteria. The most probable k1/k2, curve correlation for Venezuela described in this paper is the result of the investigation. The presented correlation defines the most probable gas-flood k,,/k,, curve through the medium of air-water capillary displacement and centrifuge water saturation tests. The laboratory procedures of these tests are. relatively simple, and inexpensive; test data should be. widely available- from routine analysis. DATA AVAILABLE, LABORATORY METHODS The report correlation utilized 107 gas-Hood k1/k2, tests run on sandstone cores of Venezuelan reservoirs. Table 1 is a general tabulation of data pertinent to the tests, while Table 2 summarizes the data. Thetests include 96 from Western Venezuela and 11 from Eastern Venezuela. Eighty-two- of the 107 test samples were sandstones that varied from poorly consolidated to-unconsolidated; 25 were consolidated. The average sample porosity was 26.7 per cent and the average permeability was 1,121 md; these values typify the better sandstone reservoirs of' Venezuela. The Welge gas-flood technique,' based on fundamental Buckley-Leverett frontal displacement theory, was introduced in about 1952 and is widely accepted in the industry. The laboratory procedure is relatively simple, rapid, and can be performed on small core samples. While there have been some minor variations in sample preparation and laboratory procedure in the tests used for the correlation, these tests can be generally summarized as follows. The core sample was first sol vent-extracted and dried. Connate-water saturation was restored by the oil-flushing or evaporation-blow down methods. At the beginning of gas flood the hydrocarbon pore volume was completeiy saturated with the test oil phase. Unsteady-state gas-oil displacement then began with the injection of nitrogen or helium. while the displaced oil and gas phases were incrementally metered at the out-flow face. From the test data, the k,,/k,, curve was calculated by the Welge method.' The individual oil and gas relative permeabilities were also calculated." CORRELATING PROCEDURES In attempting to establish a basis of correlation, we found that broad mid-range sections of 105 of the 107 k,,/k,, test curves could be closely duplicated by a straight line. Only two curves did not show a degree of linearity in this region. Correlation-curve definition parameters were subsequently developed from this observation of consistent mid-range linearity. Possible correlating variables were limited to the physical properties measured on core samples that (1) were widely available as common test data and (2) could be easily and cheaply obtained through future laboratory work. The more obvious possibilities were porosity, permeability and

Jan 1, 1966

By J. M. Smith, G. P. Willhite, J. S. Dranoff

Heat transfer rates were measured in sandstones with flow of gases perpendicular to the direction of energy transfer. Effective thermal conductivities ker ranged from 0.7 to 1.7 Btu/(br)(ft)(°F). The contribution of the solid phase appeared to be the most important in these consolidated materials, although the thermal conductivity of the gas had some effect. The velocity of the gas through the pores of the sandstones had no influence upon ker up to values of 168 lb/(br) (sq ft) in agreement with data obtained for unconsolidated beds of glass beads. The present results indicated that gas mixing, and hence heat transfer by convection in the pores, is less for perpendicular transfer of energy than when fluid flow and energy transfer are in the same direction. HEAT TRANSFER PERPENDICULAR TO FLUID FLOW IN POROUS ROCKS Heat transfer in porous media with pore sizes in the micron range depends upon the fluid in the pores and the geometry of the solid phase. The best characterized system is a bed of solid spherical particles. Heat transfer in this system has been studied extensive1y (Ref. 8 summarizes the literature up to 1959) when the pores contain stagnant fluid. When the fluid is in motion the directions of flow and energy transfer have an effect on the heat transfer rates, as demonstrated by comparing the work of Willhite, et al l3 for perpendicular flow and that of Kunii and Smith9 for parallel flow of energy and fluid. For perpendicular flow, no increase in effective thermal conductivity ke was noted up to mass velocities G of 77 lb/(hr) (sq ft). In contrast, for parallel flow ke, increased with G in the same range of flow rates. These higher values of ke in the direction of flow also have been observed3,4,10 in beds of larger particles, 0.1- to 0.5-in. diameter. For beds of consolidated materials, such as porous rocks, data are not available for these compar- isons, although Adivarahan1 reported results for the parallel case. Hence the primary objective of this work was to measure ke values for perpendicular flow of fluid and energy in porous rocks. Of interest also was the variation in effective conductivity with fluid velocity. APPARATUS AND PROCEDURE In the experimental method, a constant heat flux was applied to the inner surface of an annular section of the porous rock. By cooling the outer wall, a temperature gradient through the annular sample was established and measured with thermocouples placed within the sample at various radial positions, and at three elevations (A,B.C). The location of the 2-in. O.D., 3.75-in. long sample in the apparatus is shown in Fig. 1. Fluid entered the bottom (1) of the 3-in. I.D. (approximate) steel shell, flowed upwards through the sample and out at the top (3). pressure taps (2,4) were used to check the permeability of the sample. The energy flowed radially from the centrally-located electric heater through the sample and was absorbed in the water-cooled jacket. The samples studied were naturally occurring sandstones from different locations with the properties given in Table 1. These materials are identical with those used by Adivarahan 1 for the parallel flow of energy and fluid. Prior to use they were refluxed with toluene to remove hydrocarbons and leached with distilled water to remove soluble salts. Each sample was visually examined and discarded if large nonhomogeneities, such as cracks or stone particles, were noted. Eight copper-constantan thermocouples were inserted in holes drilled radially into the sample with

By W. Hurst

What is entailed here is the extension of the sinzplified material balance formulas to encompass interference between oil fields. As previously reported, the ex-plicitness as so revealed for the cunzulative pressure drop as a function of all factors contributing to its change in the material balance equation, is now transcribed through the inter-conzmunicating aquifer to effect an increased pressure drop on an adjoining field by interference. Such is performed by mathematical analyses and the application of the Laplace transformations. What is accomplished is that the reiteration problems previously associated with interference studies are nullified, since volumetric changes for the fluids in situ are automatically adjusted by the explicitness so expressed, and what pertains to the superposition principle applies only to the impeded water drive upon a su.bject field, which is likewise incorporated in the over-all pressure drop that results. The mathematics treats with the rigorous solution of the problem, as well as with methods easily amenable to numerical interpretation by the practicing engineer. In all cases, however, what is deduced for areal extents, reveal interference between wells when time becomes large that further substantiate the analyses. Probably of equal significance, all variables and for/nation characteristics are accounted for. This applies to the differences in PVT analyses that occur from field to field and the physical parameters associated with the lithology of the reservoirs. The latter is deduced from a rigorous interpretation of the unsteady-state flow problem for sands of different permeabilities in series. Thus, what purports to be a trial-and-error calculation to include variations in sand conditions within the intercommunicating aquifer to define water drive has little importance compared to the designation of such paranzeters at the reservoir to constitute the essential criteria. INTRODUCTION In a recent publication, the reader has been introduced to the simplified material balance formulas,' It has been shown that the complexities formerly associated for identifying and determining reservoir pressure in the material balance relationship can be resolved by treating with the Laplace transformations, the inversion of which reveals pressure as an explicit function of all factors contributing to its change. For an under saturated oil reservoir, this constitutes an integrated effect from the inception of production; and for a saturated oil reservoir, such represents a survey traverse. However, what is most important, for the first time we are able to express the pressure change in a reservoir directly to the oil, gas and water produced, and to the physical parameters of the formation, in the same manner that reservoir pressure has served in the past to define well interference — which opens up this avenue of endeavor for investigating interference between oil fields. The occurrence of interference is by virtue of the juxtaposition of two or more oil reservoirs producing in the same aquifer. Thus, in any one field, the extended production so incurred creates a pressure lowering in the field, as well as in the aquifer, which can now be transcribed by these simplified material balance formulas to the adjoining reservoir to observe the increased pressure change so induced. These effects are retroactive from field to field, to reveal a pressure drop far in excess of normal depletion to evidence the interference pattern that can result. Such is the purpose for the present undertaking — to develop the mathematical physics treating with this phenomenon of interference between oil fields, employing the simplified material balance formulas. Since Ref. 1 often will be referred to in the text, it is identified as "The Simplification", in the same connotation that another paper has been referred to as the "Laplace Trans-formation. However, the amount of material to be covered in these passages is substantial; therefore, the total emphasis will be placed upon developing these fundamental ideas treating with interference rather than attempting to offer work curves, which always can be prepared. For the practical application of these simplified material balance formulas that form the framework of the present development, the reader will find a comprehensive treatment in the Appendix of Ref. 1, illustrating a factual field example. In essence, what is proposed by this presentation is to make available to the practicing engineer a simplified means to treat with interference between oil fields, independent of the complex machine computations that have attended such undertakings in the past. As the problem unfolds, a common pattern will be discernible, whereby the many reiteration processes that previously required computing devices are eliminated.

By G. A. Zeito

Detailed visua1 examination of outcrops was used to ob-tain data on the lateral extent of shale breaks. Thirty vertical exposures belonging to maritie, deltaic and channel depositiorral environrrrents were exatmind, surveyed and photographed. The dimensions of the outcrops ranged from 356- to 8,240-ft long and 25- to 265-ft thick. Shale breaks were found to extend laterally for significant distances. and in some sands terminates by joining other break v much more frequently than by disappearance. Consequently with regard to flaw, a gross sand consisted of both continuous and discontinuous subunits. The degree of continuity of shale breaks as well as the occurrence and spatial distribution of discontinuities were different for the three depositional environments. Statistical eva1uations were performed to determine the confidence level with which estimates derived from outcrops can be applied to reservoir sands. Results of these evaluations revealed that: (I) the lateral continuity of shale breaks in marine. sands is si~nificatit, and the estimates of lateral extent can he applied to reservoir sands with a high degree of confidence (80 to 99 per cent of the shale breaks continued more than 500 ft, with a confidence of 86 per cent); and (2) the tendency for adjacent shale breaks to converge upon each other over small distances in deltaic and channel sands is highly significant (62 to 70 per cent of the shale breaks converged in less than 250 ft, with a confidence of 50 per cent), hut the probable magnitude of the resulting sand discontinuities cannot yet he predicted with adequate confidence. INTRODUCTION Almost all of the efforts devoted to characterization of the variable nature of reservoir sands have been focussed on permeability variations. Among the widely used concepts that have emerged from these efforts are those of stratified permeabilities, random permeabilities, and communicating and noncommunicating layers of different permeabilities. This study is concerned with the presence of interbedded shales and silt laminations. These features are impermeable or only slightly permeable to flow. Therefore, knowledge of the extent to which they continue laterally and the manner in which they terminate within the bodies of gross sands is important for proper description of reservoir flow. Initial field observations made on outcrops revealed that shale breaks and the relatively thinner silt laminae have impressive lateral continuity. They appeared to divide sand sections into separate individual sand layers. Although most of the layers were continuous across the total lengths of the outcrops, some were discontinuous because the- bounding shale breaks converged. Furthermore, the discontinuous layers appeared more prevalent in channel and deltaic sands than in marine sands. Based on these initial findings, a detailed investigation was carried out to determine, quantitatively: (1) the degree of continuity of shale breaks in marine. deltaic and channel sands; and (2) the frequency and spatial distribution of discontinuities in the three environments. PROCEDURE The procedure used to obtain field data from outcrops included visual examination, surveying and photographing each outcrop. The photographs were examined carefully and important outcrop features were traced, measured and recorded. The selection of outcrops for this study was made on the basis that each outcrop should be exposed clearly to permit detailed visual examination of vertical lithology. and it should also be sufficiently long (over 200 ft) to provide useful data on the lateral continuity of lithology. Identification of the depositional environment for each outcrop was made on the basis of bedding characteristics, vertical sequence of lithology and the presence of indicative sedimentary features. Whenever possible, hand specimens of associated shales were collected to determine depositional origin. Almost one-half of the outcrops used in this study required environmental identification; the remainder had already been identified by previous investigators. Several photographs of each outcrop were usually required to cover the entire length of the outcrop. These photographs were taken from one station or several, depending on the terrain, size of the outcrop and distance to the outcrop. A Hasselblad camera, with a standard 80-mm lens and a 250-mm telephoto lens, was used. The telephoto lens permitted photographing outcrops as far as two miles away. Slow-speed films were used. either Panatomic-X or Plus-X. The final operation conducted in the field was that of surveying the outcrops. The distance of an outcrop from a point of observation was determined by a triangulation method using the plane table. The measured distance was then combined with the angle of view of the camera lens to establish a scale to be used on the photographs. Films were processed using standard processing techniques and 4.5X enlargements made. The enlargements of each outcrop were butted together to form a single panorama. Slides were also prepared on several outcrops; these were used whenever greater magnification (wall projection) was required to bring out maximum lithologic detail. The shale breaks and bedding planes in each outcrop were traced on transparent acetate film superimposed on

Jan 1, 1966

By E. L. Dougherty

A phenomenological theory for a one-dimensional unstable miscible displacement similar in type to the Buckley-Leverett model but including the effects of mixing is proposed An equation giving the fractional flow of pure solvent through an oil phase containing dissolved solvent is derived including the effects of gravity and nonbomogeneities. Numer ical integration of the resulting pair of nonlinear hyperbolic partial differential equations gives volumes of dissolved and undissolved solvent as functions of space and time. Parameters in the model which characterize mixing due to dispersion were determined from experimental data; the parmeters correlated with viscosity ratio. The results indicate that the rate of dispersive mixing is proportional to volumetric flow rate Incorporated into our equations were Koval's heterogeneity factor H. which characterizes mixing due to channeling. This allowed satisfactory predictions of displacement behavior observed in horizontal floods in nonhomo-geneous cores; but in most cases the values of H which we used were considerably larger than those used by Koval. INTRODUCTION It has been shownl,2,3 that for favorable viscosity ratios, the diffusion equation with convection satisfactorily describes the behavior of a miscible displacement in a porous media. Numerous attempts to develop a satisfactory mathematical description for the case of unfavorable viscosity ratio have been reported, but they have met with only partial success. These attempts, which are reviewed by oval! have taken two approaches: (1) development of a flow model akin to the Buckley-Leverett approach for immiscible displacement neglecting the effects of mixing5,6 and (2) simultaneous application of Darcy's Law for flow and the diffusion equation with a convection term for mass transfer.7,8,9 We set out to construct a mathematical model of Type 1 including, though, the effects of mixing. Based upon a set of hypotheses, we derived a pair of nonlinear hyperbolic partial differential equations which describe in a one-dimensional system the combined effects of flow and mixing. To work with these equations it was necessary to develop a fractional flow formula. The combined system was integrated numerically using the method of characteristics. In our equations the mixing process is characterized by four parameters. Three of these, labeled ß, P1 and P2, account for dispersive type mixing. The fourth is the heterogeneity factor H, proposed by Koval to account for channeling due to nonhomogeneities in the porous media. Values of the parameters which would cause agreement between theory and experiment were determined by trial-and-error for miscible floods conducted in horizontal cores of both homogeneous and heterogeneous materials. Calculations were also performed for vertical floods in homogeneous cores. The purpose of this paper is (1) to present the details of the mathematical analysis, (2) to present the results of the calculations, (3) to consider what light the results shed on the mixing process in an unstable miscible displacement, and (4) to provide a firmer foundation for the correlation technique developed by Koval for predicting the behavior of unstable miscible floods. STATEMENT OF PROBLEM The problem is to construct a mathematical model which describes in one dimension the observed behavior in an unstable miscible displacement. The approach is phenomenological in that the equations are based on assumptions which violate certain physical precepts known to apply to the displacement phenomenon. However, the assumptions do allow us to account mathematically for the more essential phenomena. The work is of value if we can quantitatively predict experimental results which heretofore could not be predicted, even though the synthetic nature of the model belies complete explanation of the observations. We assume that the system is comprised of two contiguous flowing phases. One phase, which we call free solvent, has the properties of pure solvent. We call the fraction of the pore volume occupied by this phase the solvent saturation, designated s. The oil phase consists of oil containing dissolved