Search Documents

By W. B. Gogarty, H. R. Bailey

Methods are presented for determining molecular diffusion coefficients by using data from, capillary flow experiments. These methods are based on a numerical solution (presented in a previous paper) of the partial differential equation describing the combined mechanisms of flow and diffusion. Results from this numerical solution are given and compared with the approximate analytical solution of G. I. Taylor. The numerical solution is valid over a much larger time range. These methods are applied to experimental results for the fluid pairs water-potassium permanganate solution and amyl acetate-orthoxylene. Both of these fluid pairs have approximately equal densities and viscosities. Graphical and numerical techniques are presented for deter mining diffusion coefficients from the flow data. Values obtained by these techniques are compared with values obtained by other methods. INTRODUCTION The molecular diffusion coefficient is known to be a variable in determining the amount of mixing in a miscible displacement process. The effect of molecular diffusion on dispersion in longitudinal flow through porous media -has been examined by different investigators. These investigators concluded that at low velocities of flow, the amount of dispersion is approximate1y proportional to the molecular diffusion coefficient. The influence of diffusion on fingering, channeling, and overriding has been mentioned by other investigators.4 Recent studies5,6 have been made on the effects of molecular diffusion in connection with the problem of gravity segregation. Many different methods have been developed for the experimental determination of molecular diffusion coefficients? These methods differ mainly according to boundary conditions selected and analytical procedures used. Nevertheless, all of these methods have the condition in common that the bulk fluids in which diffusion is occurring are stationary with respect to each other. In connection with a series of papers on mixing in capillary flow,8-10 Taylor suggested the use of a flow method for determining molecular diffusion coefficients. Additional studies 11-13 have been conducted on miscible displacements in capillary tubes, but the data from these studies were not used for the specific purpose of determining diffusion coefficients. The flow method proposed by Taylor results in a single value of the diffusion coefficient for the fluid pair used in the displacement experiments. This single value represents the true value for the fluid pair when the diffusion coefficient is independent of concentration. If the diffusion coefficient is a function of concentration, the single value obtained by the flow method gives an average value for the coefficient of the fluid pair. These average values are based on diffusion taking place over the entire range of concentration, i.e., from 0 per cent of one fluid to 100 per cent of that same fluid. In field applications of the miscible displacement process, gradients occur over the same range of concentration as are found in the displacements in capillary tubes. Molecular diffusion coefficients obtained from the capillary flow method should, therefore, be especially relevant to field operations. This investigation was undertaken to evaluate the feasibility of obtaining molecular diffusion coefficients from capillary -flow experiments. In making this evaluation, diffusion coefficients were first determined for two systems from data obtained in capillary flow experiments. These values of the diffusion coefficient were then compared to values obtained by other methods. MIXING IN CAPILLARY FLOW—THEORETICAL The theoretical basis for determining molecular diffusion coefficients from capillary flow experiments is the partial differential equation relating the mechanisms of flow and diffusion. For the case of laminar flow and with a fluid pair having equal densities and viscosities, the equation is given by

By J. S. Levine, M. Prats

A homogeneous and uniform cylindrical reservoir containing oil and gas is fractured vertically on completion and is produced at a constant bottom-hole pressure. The fracture has an infinite flow capacity, is of limited lateral extent and is bounded above and below by the impermeable strata defining the vertical extent of the reservoir. Results show that such a fractured reservoir can be represented by a reservoir of circular symmetry having very nearly the same production history. The well radius of this circular reservoir is about 1/4 the fracture length and is essentially the same as that obtained previously for a single fluid of constant compressibility. At the same value of cumulative oil production, gas-oil ratios of fractured reservoirs producing at constant terzinal pressure are larger than those of reservoirs having no fractures. This leads to more inefficient use of the reservoir energy in fractured wells and results in lower reservoir pressures for the same cumulative oil production. The reduction in operating life due to fracturing a reservoir is not as great as that for a slightly compressible fluid. This diflerence can be accounted for by the lower reservoir pressure in the fractured reservoir and its adverse effect on the average mobility and compressibility of the oil. As anticipated, the reduction in operating life increases czs the reservoir permeability decreases. The type of results presented in this report can be used to determine the economic attractiveness of fracture treatments per se, to setect the initial spacing to be used in developing a field, and to compare the relative merits of fracturing available wells and infill drilling. INTRODUCTION The effect of vertical fractures on a reservoir producing either an incompressible or a compressible liquid has already been discussed in the 1iterature.l,2 Those results indicate that the production history of such a reservoir is essentially the same as that of a circular reservoir having an effective well radius of approximately one-fourth the fracture length. The present work reports on the effect of a vertical fracture on a reservoir producing two compressible fluids —oil and gas—by solution gas drive. Because of the empirical nature of the PVT and relative permeability data used to obtain the performance of such reservoirs, results can only be obtained numerically and with the aid of high-speed computers. Since reservoirs lose their radial symmetry when fractured vertically, pressure and saturation can no longer be given only in terms of distance from the well. Two coordinates (such as x and y) must now be used to describe the pressure and saturation within the reservoir, and, since we are dealing with compressible fluids, time is also a variable. Thus the solution of a vertically fractured reservoir requires finding two unknowns (pressure and saturation) in two space variables (say x and y) and in time (t). Since no means are readily and generally available for solving such problems at the present time, we have used the results of previous work1,2 to approximate the effect of a vertical fracture on a reservoir producing both oil and gas by depletion. The purpose of the present wmk, then, is to investigate the possibility of using available numerical techniques (limited at the moment to one space variable) to study the two-space-variable flow behavior resulting from a vertical fracture. Results obtained in the course of this investigation are also reported and discussed. Input and output data of the numerical methods used are given in practical units: BOPD, feet, psi, cp, and md. Results are discussed fist in terms of specific reservoir and crude properties and geometries. Later, dimensionless parameters are introduced in order to extend results to different values of some of the reservoir and fracture properties. IDEALIZATION AND DESCRIPTION OF THE FRACTURED SYSTEM It is assumed that a horizontal oil-producing layer of constant thickness and of uniform porosity and permeability is bounded above and below by impermeable strata. The reservoir has an impermeable, circular, cylindrical outer boundary of radius r,. The fracture system is represented by a single, plane, vertical fracture of limited radial extent, bounded by the impermeable matrix above and below' the producing layer (reservoir). It is assumed that there is no pressure drop in the fracture due to fluid flow. 1 indicates the general three-dimensional geometry of the fractured reservoir. Gravity effects and the effects of differential depletion resulting from variations in hydrostatic head (pressure) will be neglected. Thus, the flow behavior in the fractured reservoir is described by the

By D. R. Wieland, H. T. Kennedy

The object of the work reported here was to determine the content of elemental sulfur in gaseous methane, carbon dioxide, hydrogen sulfide, and in mixtures of these gases, at pressures and temperatwes encountered in natural gas reservoirs. Sulfur content at equilibrium is reported for pure methane, carbon dioxide, hydrogen sulfide, and on three binary mixtures of each of the three pairs of gases at pressures of 1,000, 2,000, 3,000, 4,000, 5,000 and 6,000 psia and at temperatures of 150°, 200° and 250°F. In addition, the sulfur content of three ternary mixtures at the same temperatures and pressures are reported. The results indicate that the sulfur content is higher in the gases at higher temperatures and pressures. The content is highest in hydrogen sulfide, intermediate in carbon dioxide and lowest in methane. INTRODUCTION At ordinary pressures and temperatures, the concentration of a nonvolatile material, such as sulfur in a gas at equilibrium, is a function of the vapor pressure of the material and is independent of the nature of the gas. As the pressure and temperature increase, the gas assumes some of the properties of liquids, including the power to dissolve other liquids and solids, to an extent dependent on the nature of both the gas and the material dissolved. At equilibrium, the content of sulfur in the several gases considered here may thus be considered as solubilities which are fixed for a given composition of gas, temperature and pressure. The study of the solubility of elemental sulfur in gases is of interest because sulfur is sometimes present in reservoirs producing natural gas and must be present in the vapor phase. Upon reduction of pressure and temperature, the sulfur precipitates from solution in the reservoir and in the tubing and fittings. Due to the fact that the greatest pressure drop in the reservoir is around the wellbore, the volume of sulfur precipitated will be greatest in this locality and can cause a substantial reduction in the permeability of the formation in this area. As the natural gas flows up the tubing string, the pressure and temperature are further decreased, causing further sulfur precipitation. If the volume of free sulfur is large and remedial measures are not taken, complete plugging of the tubing can occur. Knowledge of the content of sulfur as a function of composition, temperature and pressure will enable the operator to determine what changes in pressure and/or temperature are necessary to keep the sulfur in solution. If the sulfur solubility is large and can be controlled, the operator may desire to produce the well at a high wellhead pressure and temperature and reclaim the sulfur at the surface. The high solubility of sulfur in hydrogen sulfide (5.6 per cent by weight at 5,000 psi and 200°F) suggests the feasibility of recovering sulfur by injecting this gas into sulfur-bearing reservoirs and expanding the produced gas to precipitate out the sulfur. The present project was designed to measure the solubility of sulfur in carbon dioxide, methane, hydrogen sulfide and mixtures of these three gases at various temperatures and pressures. These gases were the major constituents of the gas well where plugging by sulfur was encountered. A study of the published literature shows no data on the solubility of sulfur in any gas, although Hannay and Hogarth,' in 1880 reported that it was soluble in carbon disdlide above the critical temperature of carbon disulfide. EQUIPMENT AND PROCEDURE Fig. 1 shows the layout of equipment employed to bring gases to equilibrium with sulfur and to isolate the quantity of sulfur contained in a known volume of gas. Pure gases were measured in the charging bomb at known temperatures and pressures and displaced into the reservoir bomb to make up the desired mixtures. When equilibrium was achieved

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. W. Tracy, R. D. Carter

Numerical calculations were made to determine the behavior of reservoirs with high-pressure drawdown and wide well spacing where the initial productivity is low and the wells are completed by hydraulic fracturing. The two-phase flow equations were solved for the flow into a single well. This well was assumed to be producing from a reservoir with hydraulically created horizontal fractures (four different systems with fractures were studied). For comparison purposes, additional two-phase flow calculations were made assuming a reservoir with uniform rock properties. The two-phase flow results were also compared with the conventional calculation methods, which do not include the effect of saturation gradients resulting from a simultaneous flow of oil and gas which are normal to this type reservoir. It was found that the conventional methods predict (1) a high and too optimistic value of ultimate recovery, (2) a high producing rate and a high reservoir pressure at a given oil recovery and (3) a low trend of gas-oil ratio with oil recovery. Included in the two-phase flow calculations were provisions to control the oil production rate by an allowable rate and, also, by a gas-oil ratio penalty rule. For the systems with hydraulic fractures, the producing rate was controlled by the gas-oil ratio penalty rule for most of the life. This is in contrast to the system with uniform rock properties which went "on decline" almost immediately. An unexpected characteristic of the systems which included fractures was the early rise in producing gas-oil ratio from 730 cu ft/bbl to approximately 1,200 cu ft/bbl, followed by a "leveling off" before the normally expected gas-oil ratio rise began. Additional features which are a result of hydraulic fracturing are (I) greater ultimate recovery, (2) higher average producing rates and (3) a lower average reservoir pressure at a given oil recovery. INTRODUCTION Some oil fields discovered during the past few years are producing from certain volumetric ally controlled reservoirs (often referred to as solution or internal gas- drive reservoirs) which are characterized by high-pressure drawdown at the wells. Since the available pressure drawdown at a well is limited by the static reservoir pressure and the producing rate is controlled by the available drawdown, wells completed in this type of reservoir usually produce at a rate less than the allowable from the time of completion. Because of this, this type of reservoir is referred to as a low productivity reservoir. Economic considerations require the use of wide well spacing and well stimulation by hydraulic fracturing to make commercial wells in this type of reservoir. Performance predictions for volumetrically controlled reservoirs have been made using a combination of two standard equations. 1. The "Schilthuis" or "Muskat" type material balance equation is used to relate the average reservoir pressure and the cumulative oil recovery. 2. The results from the material balance equation and the productivity factor as described by Pirson' are used to relate the cumulative recovery with producing rate and time. The material balance equation assumed uniform pressure and liquid saturation conditions throughout a reservoir. The steady-state radial flow formula allows for a pressure gradient toward a well but assumes uniform liquid saturation. These calculation methods are adequate for application to reservoirs wherein the drawdown at the well to realize satisfactory producing rates is small compared to the total pressure. In low productivity, volumetrically controlled reservoirs, the pressure drawdown at the well is large cornpared to the total pressure. Although a precise number cannot be given for the magnitude of a large pressure drawdown, values in excess of 1,000 psi would definitely be included. For practical considerations, this usually occurs when the formation flow capacity is less than about 100 md-ft. However, this limit of formation flow capacity will vary with the well producing rate. The low pressure in the neighborhood of the well which results from a high drawdown causes evolution of large volumes of gas. This causes the gas saturation to be higher near the well than at a greater distance—-hence, a non-uniform gas saturation. Also, the relationship between the relative permeability to oil (K/K) and gas saturation is nonlinear but decreases approximately in an exponential way with increases in gas saturation. Because of this, the following chain reaction is established.

By F. Selig, A. S. Odeh

A second-order approximation to the exact solution of the diffusivity equation corresponding to the pressure build-up of a well producing at a variable rate is derived. This approximation is applicable when the well's shut-in time is larger than the total time elapsed since the well was first produced. The resulting equations are compact in form and easy to use. Thus, the need for Horner's' theoretically precise but rather laborious solution to the above problem is eliminated. In addition, these equations apply where the use of Horner's widely known approximate method is questionable. From a practical point of view, the reported method is best suited for analysis of drill-stem tests and short production tests conducted on new wells. INTRODUCTION The utility of drill-stem and short production tests in reservoir studies has long been recognized by the reservoir engineer. If interpreted correctly they could lead to a wealth of information upon which may depend the success or failure of reservoirs' analyses. Initial reservoir pressure and the average flow capacity are two quantities that are normally sought from a drill-stem and/or a short production test analysis. Pressures are the most valuable and useful data in reservoir engineering. Directly or indirectly, they enter into all phases of reservoir engineering calculations. Therefore, their accurate determination is of utmost importance. The flow capacity kh of the reservoir is indicative of its commercial capability. In addition, it can indicate the presence of a darnaged zone around the wellbore and, thus, the necessity for remedial measures. Of the several methods used to analyze drill-stem and short production tests, Horner's' method is by and large the most common. It applies to an infinite reservoir and or a limited reservoir where the effect of production has not been felt by the boundary. Horner's method makes use of the so-called "point-source" solution of the diffusivity equation. The point-source solution is approximated by a logarithmic function and the superposition theorem is utilized to give the familiar pressure build-up equation where is the shut-in time, q is in reservoir barrels per day and the rest of the symbols conform with AIME nomenclature. Eq. 1 was derived for a well which produced at a constant rate q from time zero to time t and was then shut in. In actuality, such a constant rate of production does not normally obtain. Therefore, a correction must be applied to Eq. 1 to account for the varying rates of production. Horner suggested two methods. The first, which results in a theoretically accurate solution, is rather lengthy and laborious and, thus, it is not suited for routine analysis. The second which has been termed a "good working approximation" is the one used by the majority of the reservoir engineers. In the second method, Eq. 1 is modified by simply introducing a corrected time t, and writing where q is the last established production rate prior to shut-in, and t, is obtained by dividing the total cumulative production by the last established rate. Horner's original paper does not give any indication that this method of correction is based on any theoretical justification. In addition, there is a question as to what constitutes the last established rate. In case of a drill-stem test some engineers use the average rate obtained by dividing the total fluid produced by the total flow time, while others calculate the average rate by dividing the total fluid produced by the last flow-period time. Obviously, different results obtain for the different flow rates used. Because of this, a simple method to the varying-rate case was developed which is theoretically sound and which defines clearly the flow rate and its associated time to be used in the calculations. The final equation arrived at is where q* and t* are a modified rate and time, respectively, and can be easily calculated. In addition, it is shown theoretically that Horner's approximate method, if used for a variable-rate case, gives the correct pressure but would not be expected to give the correct flow capacity. MATHEMATICAL ANALYSIS The general equation governing the flow of slightly compressible fluid in porous media may be written as The elementary solution to Eq. 4, representing an instantaneous withdrawal of Q units volume of fluid at the origin at t = 0, is known as the instantaneous sink

By L. G. Jones, A. S. Odeh

A theoretical developtment is presented which provides a straightforward method of handling the drawdown analysis for both oil and gas wells flowing at variable rates. In the past our inability to analyze variable-rate drawdowns has been a major obstacle in using them for formation evaluation. This method should permit wider use of these tests in the future. INTRODUCTION Pressure build-up analyses are widely used to obtain reservoir data such as the effective flow capacity of the formation. On the other hand. although pressure drawdowns have been run for a long time, they have not been used extensively to evaluate reservoirs because they are difficult to analyze. One of the difficulties in analyzing pressure drawdowns has been that the theory developed until now has required the test to be run at a constant rate. For both oil and gas wells this is a difficult requirement to satisfy. especially during the early period of a well's life when drawdown tests are usually run. In spite of this limitation, drawdown tests are often run with build-ups as a check. In this paper a theory, which handles the drawdown analysis for variable rates, is presented for both oil and gas wells. Thus, one of the biggest deterrents to running drawdowns is removed, and they should become a much more useful engineering tool in the future. THEORETICAL CONSIDERATIONS 011, RESERVOIRS The general equation used for describing the unsteady-state radial flow of slightly compressible fluids in homogeneous porous media can be written in cylindrical coordinates as where r is the radial distance in cm, t is the time in seconds, p is the pressure at r and t in atmospheres, + is the fractional porosity, p is the viscosity in cp, c is the compressibility in vol/vol/atm, and k is the permeability in darcies. There are a number of solutions of Eq. 1 for various boundary conditions. One of these is the so-called "point-source" solution, which approximates the case of an infinite oil reservoir with a well located at r = 0 and produced at a constant rate. The point-source solution of Eq. 1 is

Jan 1, 1966

By T. F. Mueller, P. A. Witherspoon

For the case of an infinite radial system operating at constant terminal rate, the reservoir engineer often uses the "point source" solution of the diffiusivity equation to study pressure interference effects. At early times and at short distances from the inner boundary these solutions are invalid. The amount of error is often not precisely defined because of mathematical difficulties. Work presented here shows that, if dimensionless time is defined appropriately, previous solutions of the pressure equation can be displayed .as a family of curves on one chart. These solutions include the point source solution (referred to in the field of hydrology as the Theis solution) and other solutions obtained with digital computer methods. With these curves, an exact evaluation of the pressure drop within a reservoir or an aquifer can be made by the engineer. Examples of field problem solutions are presented. In most reservoirs the error involved when the Theis solution is employed is often negligible; whereas, in the calculation of interference effects in an aquifer, a substantial error can occur through such an approach. INTRODUCTION Flow equations are used in petroleum engineering to study the behavior of individual wells and reservoirs. In the case of wells, the pressure response at the wellbore face is the major point of interest; whereas, in the case of reservoirs, the pressure response at the interface of the aquifer boundary is sought. To aid in such studies, the flow equations have been solved in terms of the behavior at these two inner boundaries. Only limited work has been published in regard to the pressure conditions away from these points, i.e., within the reservoir or aquifer. Theis and Mortada are among the few who have reported on this problem. The Theis approach employs the exponential integral and is valid for pressure conditions that occur some distance away from the flow disturbance. It is derived from the concept of a point source, as opposed to a flow across some finite area. The Mortada results, on the other hand, are valid at all points within the reservoir or aquifer. They are presented in terms of dimension less ratios of the radius where the pressure is desired to the radius where the flow rate is measured. Their main use, in the past, has been in aquifer studies. The published results are presented in the form of graphs that are limited to a maximum radius ratio of 64. These graphical results are cumbersome to interpolate at non-integral radius ratios, so that one may be forced to utilize a rather involved analytical expression presented by Mortada. BASIC EQUATIONS The solutions of Mortada and Theis are both based on the diffusivity equation as applied to the case of an infinite radial system subject to a constant terminal rate. The equation is obtained by combining the material balance equation with Darcy's flow equation. The assumptions implicit in the use of this equation are as follows: (1) a single fluid is present that occupies the entire pore volume; (2) the reservoir is horizontal, homogeneous, uniform in thickness, and of infinite radial extent; (3) compressibility and viscosity of the fluid remain constant at all pressures; and (4) fluid density obeys the equation Using the diffusivity equation in situations where the above conditions do not hold will result in errors. These errors (not discussed here) but only the errors which arise in the solution of the equation itself. The diffusivity equation for the homogeneous reservoir conditions cited above can be written in cylindrical coordinates, as To obtain a dimensionless equation, so a single solution can be used for applications of different porosity, permeability and fluid properties, the following transformations are usually made: After these transformations are made, Eq. 3 .can be written in dimensionless form, as MORTADA SOLUTION One solution of Eq. 6 has been given by Mortada.

Jan 1, 1966

By A. N. Crownover, H. T. Kennedy, E. P. Miesch, C. H. Bowman

The paper presents correlations of (l) molar volume of gaseous hydrocarbon mixtures with pressure, temperature, composition and properties of the C,-plus fraction; (2) shrinkage of oils during flash and differential liberation of gas, including the calculation of formation volume factor under various conditions; and (3) bubble-point pressure with temperature, composition and characteristics of C, plus. The data on which the correlations are based comprise 1,615 measurements on 900 hydrocarbon systems, including numerous systems containing nitrogen, hydrogen sul-fide and carbon dioxide. In each correlation, the number of data points covered and the accuracy is substantially greater than in previously available work. Thus, the equation yielding molar volumes of gases has an average deviation of 2.04 per cent, applied to mixtures having temperatures up to 313F and pressures up to 9,800 psis, compared to 2.37 per cent for the Benedict-Webb-Rubin equation applied to the same data, and 4.53 per cent for the method based on the law of corresponding states. he equations resented are all explicit in the dependent variable, and require no iteration on the digital cow puter. INTRODUCTION The ease and accuracy of determining the composition of hydrocarbon mixtures, compared to the difficulty of measuring their properties under reservoir conditions, makes it desirable to utilize composition as the key to physical behavior to the greatest possible extent. As a result, there are available correlations between composition, or easily measured characteristics dependent on composition, and practically every important engineering property of reservoir fluids. The task confronting us is one of finding more exact relationships between important variables rather than extending correlations to new properties. This paper describes new correlations of molar volumes of gases, formation volume factors, and bubble-point pressures with composition, temperature and where possible, pressure. Each correlation is obtained by employing a sufficiently large amount of data so the calculated properties are probably as least as accurate as the measurement on which they are based. MOLAR VOLUME OF GASES Although many equations of state have been proposed for pure gases, only a few methods are applicable to hydrocarbons at conditions comparable to those in petroleum reservoirs. Still fewer arc useful in describing the behavior of mixtures, with which the petroleum engineer is largely concerned. The correlation presented here involves procedures similar to those of Alani and Kennedy The van der Waals equation is modified to make a and b functions of temperature instead of being constants for each material.This is a cubic equation, which may have either one or three real roots. The lowest root corresponds to liquid volume, while the highest applies to gas. When the procedure was applied to 703 pressure-temperature points of 164 gases of known composition and volume, the average deviation was 12.08 per cent, and the standard deviation 8.15. The above calculations were made using the constants derived for liquids by Alani and Kennedy and the mixture equation developed by them. A closer approximation to measured molar volumes is obtained by employing different sets of constants for different areas on the pressure-temperature chart, and by changing the relationship between a,,, and b,,, for mixtures and the a, and b, for individual hydrocarbons. The various sets are shown in Table 1, and the areas for which they are recommended are plotted in Figs. 1 through 7. In these figures, the area designated by zero is in the critical region for pure materials, and the values obtained for them may be unreliable. For normal mixtures, however, the area may be added to any of the bounding areas without substantial error. The new relationships employed are

Jan 1, 1966

By L. A. Wilson, N. H. Harrison, R. R. Clay, R. L. Reed. D. W. Reed

Experiments have been performed in a linear near-adiabatic system for the purpose of extending data on reverse and forward combustion from atmospheric pressure to 1,000 psig. Results obtained from reverse combustion appear to conform qualitatively with the existing description of the process when reaction kinetics are suitably modified to account for the increase in pressure, i.e., increasing tbe pressure decreases the peak temperature and increases the combustion zone velocity. Forward combustion appears to be a fuel dominated process wherein peak temperature and combustion zone velocity are not very sensitive to changes in pressure. The moderate effects of pressure that do exist at low flux virtually disappear at high flux providing all oxygen is consumed. With this provision. increasing the pressure decreases the frontal velocity and increases the peak temperature. Results me shown graphicaIIy which demonstrate the effects of presswe on peak temperature, rate of advance, oil recovery, air/oil ratio, carbon oxides produced and temperature distributions. INTRODUCTION Numerous field tests of oil recovery using the technique of underground combustion are now in progress. Operating pressures used are always substantially in excess of atmospheric pressure. Nevertheless, the literature contains no laboratory data pertaining to the effects of pressure except for those of Martin, Alexander and Dew.1 Unfortunately, these data appear to reflect heat losses which may have obscured some of the effects of increased pressure. The underground combustion processes are exceedingly complex, and general concepts relating to them have necessarily been described in relatively simple terms. Particularly with regard to forward combustion, the interaction between multiphase fluid flow, heat transport and various rate mechanisms is so involved that to date a rigorous treatment is not available. In this paper a linear near-adiabatic system is described which was found capable of operating at pressures up to 1,000 psig and temperatures of at least 1,100 F. Using this equipment, experiments on both forward and reverse combustion were accomplished at elevated pressures and the data presented as a function of air flux or temperature. It will be assumed that the reader is familiar with the papers of Wilson, Wygal, Reed, Gergins and Henderson on forward combustion and that of Reed, Reed and Tracht3 on reverse combustion, so that it is not necessary to repeat the details of the two processes. EQUIPMENT AND PROCEDURE It was decided to study the effects of pressure over the interval 0-1,000 psig using a near-adiabatic system. Since satisfactory operation of such a system requires very thin tubing walls, it proved necessary to insert the entire combustion tube and heater assembly into a high pressure jacket in such a way that the internal tube pressure could be counterbalanced by pressurizing the annulus. This meant that the numerous heater and thermocouple wires had to pass through gas-tight seals that required considerable effort to maintain effective. An unexpected difficulty arose in connection with the nature of the annulus pressuring gas. It turned out that the composition of this gas critically affected the functioning of the adiabatic control system. Operation of the equipment at elevated pressures was not particularly difficult. The problems that did arise were a direct consequence of the large numbers of connections and seals involved. EQUIPMENT The equipment consisted of the combustion tube which fitted within the pressure jacket. This jacket was in turn mounted on a wheeled carriage which facilitated movement and provided for rotation of the

By F. E. Crane, H. A. Kendall, G. H. F. Gardner

Two effects of a difference in density may be observed in miscible displacement. As the difference is increased from zero there is first a change in mode from multiple fingers to a single finger. A further increase results in a more rapid growth of the single finger. Miscible displacement with multiple fingers has been discussed by Blackwell, et al, whereas miscible displacement with a single finger has been discussed by Craig, et al. In the present paper experiments are described which bridge the gap between these extreme cases. INTRODUCTION Craig, et al,1 investigated the effect of gravity on frontal displacements for circumstances in which the gravity forces were significant. Blackwell, et a1,2 dealt with the other extreme in which gravity forces were negligible. The considerable difference between their results is illustrated in Fig. 1 where the breakthrough recoveries are plotted against the mobility ratio. The experiments described in the present paper were conducted to examine how the transition from the one extreme case to the other occurs when the gravity forces are gradually reduced. The transition is well illustrated by the production history after breakthrough. Fig. 2 shows how the number of pore volumes which must be injected at a constant rate to displace 98 per cent of the original fluid depends on the gravity force. For each mobility ratio there is a transition range, and a very significant change in the number of pore volumes occurs when this range is traversed. This range of the gravity force is physically linked with the formation of viscous fingers. A single finger is predominant for larger gravity forces; a multiplicity of fingers is formed for smaller gravity forces. Improved recovery and the occurrence of multiple fingers is concomitant with a large increase in the amount by which the fluids are mixed within the porous medium. The transition range depends on all the dimension-less groups which characterize the problem; but this dependence has not been exhaustively studied as yet and only fragmentary results are discussed in this paper.

By D. H. Gray, G. Bergamini, I. Fatt

A method developed by soil scientists for measuring permeability anisotropy of soils has been adapted for use in studies of sandstone permeability aniso-tmpy. Studies on Berea and Bandera sandstone show that permeability of these materials decreases upon application of simulated overburden pressure. Permeability reduction of Berea and Grubb sandstone cores subjected to simulated overburden pressure is shown to be a function of the ratio of radial to axial stress. INTRODUCTION The effects of mechanical stress on reservoir rock parameters are now generally accepted as being important in reservoir engineering calculations. In several oil fields the unique behavior of the reservoir has been attributed to the influence of overburden pressure or internal fluid pressure on rock permeability or fluid content. l Although the gross effects of overburden pressure on permeability, porosity, waterflood efficiency, electrical resistance and sonic velocity have now been established,3-9 there are some very interesting, and perhaps important, details of rock behavior under stress that have not heretofore been studied. In this paper we will show that permeability anistropy of several sandstones is a function of overburden pressure. We will also show that permeability reduction of cylindrical samples of two sandstones, when subjected to mechanical stress, is a function of the ratio of radial to axial stress. Most studies of overburden pressure effects made in the past have attempted to duplicate reservoir conditions by applying hydrostatic external pressure on all surfaces of the sample. Assuming perfect elasticity, that is, no viscoelastic flow during geological time, reservoir rock in situ cannot be under hydrostatic stress as is apparent from the value of 0.25 for Poisson's ratio for most rock. In the absence of tectonic forces the horizontal stress component If there are directional permeability effects, as has been reported from field and laboratory studies, l0,11 then one would expect hydrostatic stress as used in the laboratory to give a different permeability reduction from that occurring under reservoir stress conditions. The influence of stress on directional permeability was measured by a modification of the method first described by Maasland and Kirkham for use on soil samples.12 Their method permits calculation of the permeabilities along each of the three mutually perpendicular axes of a cubic sample. The modification reported in this paper extends the method to observing these permeabilities when the sample is under hydrostatic stress. The Hassler sleeve type of core holder is satisfactory for studying overburden pressure effects on permeability if a known force is applied to the end plugs. Wyble6 used a sleeve core holder but did not load the end plugs. In this paper we combine sleeve pressure and axial force to give permeability changes which we believe are close to those occurring in reservoirs when reservoir pressure changes. The combination of sleeve pressure and axial force assumes, however, that the sleeve pressure simulates both vertical stress and horizontal stress in one direction. METHOD PERMEABILITY ANISOTROPY The Maasland and Kirkhaml2 method of measuring permeability anisotropy uses a cube of porous material which is Open On three faces and has discs with a Center hole for an air inlet cemented to the other three mutually perpendicular faces (Fig. 1). The rate at which air pressure falls in a bottle of known volume, and pressure slightly above atmospheric, when the bottle discharges through each disc, is

By F. W. Jessen, N. Lamont

The effect of an existing fracture or joint plane, which may exist in a rock, on the extension of a hydraulically induced fracture through the rock has been investigated in the laboratory. By use of a series of models made from various outcrop rocks, the width and orientation of existing fractures do not alter the extension or direction of the hydraulic fracture. The results, which conform quite well with the Grifith theory of failure, are illustrated by a series of photographs of the various models. INTRODUCTION The process of hydraulic fracturing has been widely used in the oil industry since its introduction in 1948 and has made possible the production of many reservoirs which would have been uneconomical prior to this process. Numerous studies of the mechanics of the process2-h nd of the effects of oriented fractures on recovery'.' have appeared in the literature. The increased recoveries led to attempts to predict the orientation of hydraulic fractures at the wellbore and to the development of methods to control this orientation. Field tests have been devised"." which indicate the validity of theoretical predictions of fracture orientation at the borehole. The direction of extension of hydraulic fractures from the borehole has not received much attention since most studies have predicted, at least implicitly, that a fracture will continue in the same plane unless a change of state of stress in the rock occurs. Since such predictions are based on the assumption of rock homogeneity, the effects of rock heterogeneities are left unknown. Most sedimentary rocks are composed of layers which reflect the changing depositional conditions of geological time. In addition, the more competent rocks frequently are fractured and jointed as the result of structural deformations or tectonic movements. Evidence of such joint systems exists in many surface outcrops, and it can be assumed that similar systems occur in many subsurface rocks, although the individual joints may be rather tightly closed due to overburden forces. This paper investigates the effect that an existing fracture in a rock would have on the direction or orientation of an extending or advancing hydraulic fracture when it intersects the existing fracture. The study was conducted on small rock models under triaxial stress conditions. The results indicate that an existing fracture will have little effect on the hydraulic fracture. THE MODEL STUDY The rock model represents an elemental portion of the earth containing an existing fracture which is located at some distance from a borehole from which hydraulic fracture is being extended. A generalized view of a model is presented in Fig. 1. The models were constructed from cement blocks and natural rocks. The various materials tested and some of the physical properties are listed in Table 1. The rocks were cut into rectangular blocks with dimensions of 11/2 X3Y2 X4 to 8 in. An initial slot was cut into the end face of each block along the longitudinal axis normal to the top face. The slot was filled with plastic aluminum along its outer edge to a depth of 1/16 in. A hole of 3/8-in. diameter was drilled into the sealed slot at the center of the end face to serve as an entry port for the fracturing fluid. Existing fractures of various types were simulated in the models. A "hairline" fracture with essentially zero width was created by cutting the model into two parts with a diamond saw and replacing the two parts together along the cut. A finite-width fracture was created by placing a layer of sand grains between the two faces of the cut just described (see Fig. 12). A large open fracture (similar, perhaps, to a vug) was created by cutting the model in two parts and removing a portion of the downstream block to a depth of % in. See Fig. 13. For a fourth

By C. Chu

An investigation has been made of the radial heat-wave process using a mathematical model in two-dimensional cylindrical coordinates. This model considers combustion, convection and conduction inside the reservoir, but only conduction in the bounding formations. From a study of the general features of the process, an important phenomenon has been revealed, namely, the feedback of heat into the reservoir on the trailing edge of the heat wave. The effects of various process variables on the performance characteristics of the process have also been investigated. It was found that up to the time when the combustion front reaches a given point, the per cent heat loss, provided it is not higher than 40 per cent, is approximately directly proportional to the square root of thermal conductivity arid fuel content, but inversely proportional to the square root of gas-injection rate and oxygen concentration. The effecr of reservoir thickness is more pronounced, since halving the thickness doubles the per cent heat loss. The most decisive factor in determining the center-plane peak temperature is the fuel content of the reservoir. Within the temperature range investigated, doubling the fuel content doubles the peak temperature in the early stage, but the rate of decline of the peak temperature is high. Reservoir thickness is also a very influential factor. The peak temperature is lowered when the thickness is reduced; however, the effect of thickness becomes less pronounced when the thickness is high. Reduction of oxygen concentration increases the peak temperature in the early stage but lowers it afterwards because of the higher rate of decline of the peak temperature. Increase in gas injection rate or decredse in thermal conductivity geives a higher peak temperature which stays high for a longer period. The propagation range of the heat wave is chiefly governed by the fuel content of the reservoir. An increase of 0.2 1b/cu ft in the fuel content increases the propagation range by 100 per cent. The propagation range is more than doubled by doubling the gas injection rate, or reservoir thickness, or by reducing the thermal conductivity by 50 per cent. Comparatively, oxygen concentration has less effect on the propagation range. INTRODUCTION Several investigators have conducted theoretical studies of a radial heat wave. Vogel and Krueger1 studied a system with a moving cylindrical heat source of constant temper- ature, considering conduction in the radial direction only. Ramey2 included conduction in the vertical direction in his studies. Bailey and Larkin2 attacked a more general problem where initial well heating, vertical heat losses and arbitrary frontal velocities were included. In all these studies, however, conduction was considered to be the only means of heat transfer. Bailey and Larkin in a later paper included the effects of convection in a study involving both linear and radial geometries. Vertical heat losses were neglected in the radial case. Katz5 studied a similar problem in a one-dimensional radial model, using a heat-loss coefficient to account for vertical heat losses. Selig and Couch6 mployed a cylindrical model and investigated two limiting cases. In one case they considered no heat loss from the reservoir whereas in the other they assumed a constant temperature at the interface between the reservoir and its bounding formations. Thomas' studied a more general case but assumed a permeable bounding formation so that the convection effect is not confined to the reservoir. In the present work a more realistic and more generalized model is used. It involves a two-dimensional cylindrical system with combustion, convection and conduction inside the reservoir, but only conduction in the bounding formations. The purpose is to establish the temperature distribution both inside and outside the reservoir, to study the general features of the radial heat wave process, and to investigate the effects of various process variables on the performance characteristics of the process. THEORY We fist consider a circular porous reservoir of thickness H extending vertically from y = — H/2 to y = + H/2. The reservoir extends from a well bore radius r, to an external radius re. A stream of oxygen-containing gas is introduced into the reservoir through the wellbore. The oxygen-containing gas reacts with the fuel contained in the reservoir and forms a combustion front wherever the prevailing temperature can support the combustion. It is here assumed that this combustion front constitutes a cylindrical surface source of heat having an infinitesimal thickness in the radial direction and extending vertically throughout the whole thickness H. This is designated as Region I. We next consider a Region II corresponding to the upper and lower formations bounding the reservoir, extending from y = — - to y = — H/2 and from y = + H/2 to y = + m. Since r, is very small, we may assume that the two bounding formations have the same dimensions, symmetric with respect to the center plane of the reservoir. In this way, we may take the upper half of the system alone into consideration. In contrast with

By R. G. Hawthorne

This report is concerned with fluid displacement in porous media, in those cases where viscous and gravitational forces control the displacement. Such a system would usually be found in a sand body of large physical dimensions such as an oil reservoir, although it is possible to create such a system in the laboratory. It is shown that the position of the fluid interface can be predicted by numerical calculations using a basic idea presented by Dietz. Fluid flow is considered in a vertical plane in a homogeneous, porous medium of sufficient thickness that the capillary transition zone is small in comparison with the total reservoir. A theory developed by Dietz' is used to make numerical calculations of the position of the fluid interface. The results for several conditions are compared with scaled model experiments. The results show that, for gas drive in a reservoir of steep dip, a relatively low flow rate can displace large volumes of oil before gas breakthrough. On the other hand, water injection at favorable mobility ratio and low dip may show best performance at high rates. Water tends to underride the oil and, given sufficient rime, will break through without much oil displacement. For certain conditions, which include relatively low flow rate, the interface is a straight line and its behavior is simple to calculate. At higher flow rates, the interface is unstable, and a numerical solution was programed for an automatic computer. In general, good agreement is shown between the fluid model and the computed results so long as gravitational forces have control. For a water drive at very unfavorable mobility ratio, many small water fingers appear. These viscous fingers are not controlled by the relatively small gravitational forces. When viscous fingering becomes the controlling factor, the mathematical model is oversimplified, and results do not check the fluid flow model. INTRODUCTION Present methods of reservoir analysis depend upon certain simplifying assumptions to obtain mathematical descriptions of practical use. Material-balance methods (Muskat2 or Tarner2) assume uniform fluid saturations in the entire reservoir, or in a few subdivisions of the reservoir. An unsteady-state flow calculation by West, er at considered pressure and saturation changes in flow to a well during solution gas drive, and neglected gravity effects. Results showed only a 4 per cent difference for ultimate oil recovery by the Muskat method, even though the case chosen for study was one in which unsteady-state effects should be high. The Buckley-Leverett5 method commonly assumes a one-dimensional flow system. It is applicable at high flow rates where viscous forces predominate over gravity forces. Simultaneous, parallel flow of the two fluids is assumed, and the concept of a fluid interface is not introduced. Permeabilities to each fluid for a given saturation must be known. The method is not applicable for a two-dimensional system where cross flow becomes possible. Less well known is the displacement equation derived by Dietz. This method is designed for two-dimensional flow systems and assumes a definable fluid interface within the porous medium. Dietz showed that, for a range of low flow rates, the interface would be stable. straight and at an angle of inclination which could be simply calculated. At a certain critical flow rate, the calculated interface tilt would equal the formation dip. For higher flow rates, a finger of displacing fluid would invade the displaced fluid. Dietz indicated that his method applied only to macroscopic reservoir behavior, while the Buckley-Leverett method applied to the small transition zone at the fluid interface. The examples worked out in this report are based on the fluid-displacement theory of Dietz. It is shown that the Dietz theory may be used to derive equations analogous to the Buckley-Leverett equations. In contrast to the Buckley-Leverett method, flow is considered in a plane rather than being limited to a line. Rather than a frontal advance, the movement of a fluid interface is followed. For flow rates substantially exceeding the critical rate and for high viscosity ratio, many fingers of invading fluid occur-—rather than the single finger assumed by Dietz. On the other hand, so long as some gravitational influence remains, the flow is not entirely parallel to the bedding planes as assumed by Buckley and Leverett; therefore, both methods fail to give an adequate descrip-

By M. A. Jones

Performance of a pilot flood in eastern Kansas indicates improved recovery and accelerated production resulting from mobility ratio control obtained by adding a high moleculur weight polymer to injected water. Ultimate pilot flood recovery of the 75-cp reservoir oil is anticipated to he 350 bbllacre-ft during the six-year flood life. Cumulative recovery to date totals 242 bbl/acre-ft after 30 months of injection. Reservoir parameters and flood performance are disc~rssed in detail to provide a case history relating to the feasibility of polymer flooding. INTRODUCTION The Vernon Polymer flood of the Brazos Oil & Gas Co., located in Sec. 2, 24 S., 16 E., Woodson County, Kans., was begun in late Oct., 1963. Injection of a 500 ppm solution of Pusher* polymer was initiated in a 15-acre pilot to determine feasibility of improving oil recovery by mobility ratio control. Performance to date and the anticipated recovery of 115,000 bbl from the 329 acre-ft pilot indicate efficient displacement of the viscous reservoir oil. Dependence of flood recovery on the displacing fluid-displaced oil mobility ratio has long been recognized. Displacement efficiency of the reservoir oil decreases with increasing mobility ratio. Therefore, economic reduction of this parameter is an attractive approach to improving oil recovery. particularly from viscous oil reservoirs where mobility ratios are often quite high. Waterflood mobility ratios can be controlled by dilute solutions of certain polymers'; which provide a displacing fluid with considerably reduced mobility. Reservoir parameters and flood performance are emphasized in this paper to provide data relating to polymer flood applicability. THEORY The polymer used in the Vernon pilot exhibits the property of unusually high resistance to flow through porous media.'.' A small concentration in injected water effects a mobility reduction considerably beyond that predicted from increases in bulk viscosity as measured in a capillary or rotating cup viscosimeter. Water mobility in the subject flood is reduced by a factor of 7.55 with the addition of 500 ppm Pusher chemical to injected water while the bulk viscosity is increased from 1.03 to only 1.40 cp at 75F. The exact flow mechanism causing resistance to flow has not been established; however, it is observed only in tortuous flow systems and appears to be influenced by brine properties, rock properties and oil saturation. Consequently, mobility reduction is estimated from flow studies for a given sample of reservoir rock and is expressed in terms of a "resistance factor".' This parameter is defined as the ratio of brine mobility to polymer solution mobility measured at residual oil saturation: RF= ?m/?p= (km/µm)/(kp/µp)M.......(1) where RF = resistance factor (dimensionless) k*. and A, = mobility of water and polymer solution, respectively k,,, and k,, = water and polymer permeability, respectively (darcies) µie and µp = viscosity of water and polymer, respectively (poise). Eq. 1 is similar to the expression commonly used to relate the relative mobility of the water bank at residual oil saturation to oil-bank mobility at irreducible water saturation: where M =mobility ratio of water to oil (dimension-less) A, = mobility of oil k, = permeability to oil (darcies) )i* = oil viscosity (poise). From Eqs. 1 and 2 it is apparent that the polymer-solution-to-oil-mobility ratio can be approximated by dividing the brine-to-oil ratio by the resistance factor RF: Mpn=M10/RF...........(3) where M,-, = mobility ratio of polymer solution to oil (dimensionless). This simplified presentation assumes no dilution of the polymer solution by connate water. Also, the presence of a connate water bank separating the polymer solution from the displaced oil is not considered. Displacement of oil by waterflooding is increasingly dominated by viscous forces as the viscosity of the oil displaced increases. When the oil mobility is less than that of the injected water, the oil-water interface is not a piston-like front.3,' Instead, instabilities occur at the flood front in the form of viscous fingers which tend to penetrate the less mobile oil bank. The magnitude of these protrusions, or the degree of channeling, increases

Jan 1, 1967

By A. J. Conelius, D. P. Sobocinski

This paper presents a correlation for predicting the behavior of a water cone as it builds from the static water-oil contact to breakthrough conditions. The correlation is partly empirical and involves dimensionless groups of reservoir and fluid properties and of production and well characteristics. The groups were deduced from the scaling criteria for the immiscible displacement of oil by water. The correlation is based on a limited amount of experimenfal data from a laboratory sand-packed model and on results from a computer program for a two-dimensional, incompressible system. Because the correlating groups are dimensionless, they can be used to estimate the performance of water coning cases not specifically considered in the correlation. However, despite its dimensionless nature, the correlation is not completely general and will not provide meaningful estimates of cone behavior in many situations. INTRODUCTION AND BACKGROUND The production of water from oil wells is a common occurrence which increases the cost of producing operations and may reduce the efficiency of the depletion mechanism and the recovery of reserves. We will deal with one cause of this water production, namely, coning. The coning of water into a producing well is caused by pressure gradients established around the wellbore by the production of fluids from the well. These pressure gradients can raise the water-oil contact near the well where gradients are most severe. Gravity forces that arise from fluid density differences counterbalance the flowing pressure gradients and tend to keep the water out of the oil zone. Therefore, at any given time, there is a balance between gravitational and viscous forces at points on and away from the completion interval. When the dynamic forces at the wellbore exceed gravitational forces, a cone of water will ultimately break into the well to produce water along with the oil. We can expand on this basic visualization of coning by introducing the concepts of stable cone, unstable cone and critical production rate. For instance, if a well is produced at a constant rate and the pressure gradients in the drainage system have become constant, a steady-state condition is reached. If, at this condition, the dynamic forces at the well are less than the gravity forces, then the water or gas cone that has formed will not extend to the well. Moreover, the cone will neither advance nor recede, thus establishing what is known as a stable cone. Conversely, if the pressure in the system is in an unsteady-state condition, then an unstable cone will continue to advance until steady-state conditions prevail. If the flowing pressure drop at the well is suffcient to overcome the gravity forces, the unstable cone will grow and ultimately break into the well. It is important to note that in a realistic sense, stable cones may only be "pseudostable" because the drainage system and pressure distribution generally change. For example, with reservoir depletion, the water-oil contact may advance toward the completion interval, thereby increasing chances for coning. As another example, reduced productivity due to well damage requires a corresponding increase in the flowing pressure drop LO maintain a given production rate. This increase in pressure drop may force an otherwise stable cone into a well. The critical production rate, well known in the literature, is the rate above which the flowing pressure gradient at the well causes water (or gas) to cone into the well. It is, therefore, the maximum rate of oil production without concurrent production of the displacing phase by coning. At the critical rate, a built-up cone is stable but is at a position of incipient breakthrough. Numerous papers have been published on critical rates. Some of the better known of these include the work of (1) Muskat and Wyckoff,' who first dealt with the coning problem; (2) Chaney, et al,' who developed expressions similar to those of Muskat but who presented results in a conven ient-to-use graphical form (the "Sun" method); and (3) Meyer and Garder,b hose analysis is based on radial-flow formulas. One assumption in critical production rate analyses is that the cone has built-up to just before its breakthrough into the well. But, these analyses reveal nothing directly about the time it takes for the cone to build up to this incipient breakthrough position. Thus, water-free oil can be produced from a well for prolonged periods at rates above the critical rate before the well reaches the condition to which the critical rate applies. The published literature contains little on the rate of growth of a cone. Experimentally, Meyer and Searcy studied the rate of rise of a cone in a Hele-Shaw model.' Additional related work on water breakthrough and produced water-oil ratios in water driven reservoirs was reported by Muskat," Hutchinson and Kemp," Henley, et al,' and Stevens, et a1.B Theoretically, the basic coning equations for a water-oil system can be developed by applying the conservation of mass to each of the phases, relating flow velocities with pressure by Darcy's law, and relating pressures across water-oil interfaces by capillary pressure. With the usual boundaries at the well and reservoir limits, the solution of the resulting equations for the time behavior of a water-oil interface constitutes a free-surface, boundary-value problem. The shape of the boundary depends on the internal potential distribution, while this

Jan 1, 1966

By G. L. Chierici

Starting from Muskat's theory of water and gas coning, maximum permissible oil production rates without water and/or free-gas production have been determined, in a broad range of reservoir and well parameters, using the potentiometric model technique. The main assumptions made are as follows: (1) the reservoir rock is homogeneeous (either isotropic or anisotropic); (2) the volume of the aquifer underlying the oil zone is very small, so that it does not contribute to reservoir energy; and (3) the gas cap expands at a very low rate, so that it can be assumed to be in quasi-static conditions. The results obtained are presented in the form of diagrams which can be used for solving two types of problems: (1) given the reservoir and fluids characteristics, as well as the position and length of the perforated interval, determine the maximum oil production rate without water and/or free-gas production; and (2) given the reservoir and fluids characteristics only, determine the position and length of the perforated interval which optimize the maximum permissible oil production rate, without water and/or free-gas production. INTRODUCTION In oil reservoirs where the oil-bearing formation is underlain by an aquifer which does not participate in the production mechanism, water-coning is a limiting factor to the flow rates of producing wells. Production rates are usually kept to a value that will prevent the water from entering the wells. The entry of water into a well lowers its productivity by increasing the weigbt of the fluid column; moreover, the separation of water from the effluent, at the surface, may constitute a very difficult problem in cases of heavy viscous oils. A similar situation is encountered in oil reservoirs with a gas cap overlying the oil-saturated zone; here a downward gas cone is induced by the flow of oil towards the producing wells. Production rates must be low enough to prevent the gas from being produced; producing gas from the gas cap would be a waste of energy. Of course, water-coning and gas-coning phenomena can occur at the same time in the same reservoir if the oil-producing formation is both overlain by a gas zone and underlain by a water zone. Due to its relevant practical importance, the mechanism of coning was studied by many people.2,3,5-8 Defining the conditions for getting the maximum water-free and/or gas-free oil production rate is a difficult problem, often encountered under one of the following aspects: 1. Predict the maximum flow rate that can be assigned to a completed well without the simultaneous production of water and/or free-gas. 2. Define the optimum length and position of the interval to be perforated in a well, in order to obtain the maximum water and gas-free production rate. A systematic study of these problems was made by means of the electrical analog technique. The results of this study are presented here, under the form of a set of curves providing solutions for the above stated problems. These curves are valid only for homogeneous forrnations, either isotropic or anisotropic. Should the formation be non-homogeneous (by horizontal or vertical variation of permeability, shale diaphragms, fractures, etc.), a specific potentiometric study would be required for each specific case. Especially when shale diaphragms of some radial extension are present, the critical rates observed are much larger than would be expected from the diagrams. STATEMENT OF THE PROBLEM In the present study the aquifer is supposed to be of such limited volume that it does not contribute to the energy of the reservoir. Moreover, the gas cap is supposed to expand at such a low rate that the potential gradient in the gas cap is negligible. Under static conditions water-oil and gas-oil interfaces (T1 and T2) are both horizontal. When the reservoir production starts, below each well these interfaces take a cone-like shape (Fig. 1) having as an axis the axis of the well. This shape results from the equilibrium between potential gradients in the oil zone and gravitational forces due to density differences between oil and water and between oil and gas. Assuming the oil-bearing formation to be homogeneous and the oil to be incompressible, the analysis of the problem (see Appendix) shows that the oil-water and gas-oil interfaces are stable only if the oil production rate of the well is not higher than the following values.

Jan 1, 1965