Search Documents

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 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 D. A. Flanagan, P. B. Crawford

A study has been made of the feasibility of storing liquid meihane at low pressures in undergrohd caverns. Methane liquefies at — 258°F at atmospheric pressure. It is shown that the methane evaporation rates will rapidly decrease and cool the surrounding rock so that at the end of one month they would be between 50 Mcf/hr for a cavern of 25-ft radius and 700 Mcf/hr for a cavern of 100-ft radius. At the end of 10 years, the evaporation rates would be 18 and 100 Mcf/hr, respectively, for caverns of the same radii. The evaporation rates may be reduced by a factor of 2 to 10 by the application of insulation. The cost for mining the caverns is estimated to be $.75 to $1.25/Mcf of storage. This is substantially less than surface storage; it is believed to be safer and to result in lower maintenance, savings in space and savings in strategic materials. INTRODUCTION During the past few years, there has been an increasing interest in the economic feasibility of liquefying methane. Methane liquefaction is being considered for tanker transport; in addition, liquefaction is being reconsidered for shaving peak gas demands. Several articles have described natural gas liquefaction and the progress of the tanker in making trial runs from the United States to Great Britain to determine the feasibility of tanker transport. At the present time, pipelines are not designed to supply peak gas loads during extremely cold periods such as are oiten encountered in the North and Northeast. Gas is being stored in underground reservoirs en-route to its destination, but in many instances satisfactory storage in porous reservoirs has not been practical, especially along the Eastern seaboard where few petroleum reservoirs have been found. In England and other foreign countries, it is unlikely that satisfactory porous structures could be found, and it may be desirable to mine or excavate the rock to obtain storage. By storing the methane near the consumption point, product availability can be increased during periods of need. Methane liquefies at - 258OF at atmospheric pressure, and 1 cu ft of liquid methane will make about 600 cu ft of gaseous methane. Methane liquefaction for shaving peak gas demands was conducted in Cleveland, Ohio in the early 1940's. The Cleveland plant was designed to liquefy 4 MMcf/D of natural gas.' Development of several low-temperature processes has resulted in some improvements in the design of low-temperature storage vessels. Liquefied-methane storage vessels can be designed for various evaporation rates. The desired rate will depend on the reason for liquefaction. For shaving peak gas demands, the desired evaporation rates would be low, possibly between 0.2 and 1.5 per cent/day. For methane tanker transport and supply, the evaporation rates must be high and must equal the total daily rate of supply. It is apparent that the cost of storage facilities will depend on the objective and will vary appreciably. Surface storage of methane may cost between $3,000 and $20,000/MMscf, depending on the size and design. The failure of the surface storage vessels at Cleveland, with the resulting injury and loss of many lives, has caused considerable emphasis to be placed upon developing new, safer and improved methods of storing liquid methane. Storage of liquefied propane and butane in underground salt domes has been very satisfactory. The per-barrel cost of underground storage for any sizeable capacity is very small compared to the cost of surface storage for LPG mixtures. The larger the capacity of the underground storage, the less is the per-barrel cost. Numerous advantages appear to exist if liquefied methane could be satisfactorily stored underground, the principal advantages being increased safety, lower initial cost, lower maintenance cost, savings in space and savings in strategic materials. The purpose of this paper is to determine the feasibility of storing liquid methane underground in mined caverns. If a large spherical cavern is dug underground and is filled instantaneously with liquid methane, the surface of the sphere may approach the liquid methane temperature almost instantaneously. The temperature distribution for a sphere of radius r. in an infinite medium, initially at one temperature with spherical surface kept at T, from time t = 0, is given by Eq. 1.'

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 A. Madrazo

The Standing-Katz method for predicting liquid densities of reservoir fluids has been tested using experimental data of 154 bottom-hole or recombined reservoir fluid samples. New pressure- and temperature-correction curves for the Standing-Katz type correlation are presented which improve the accuracy of prediction. INTRODUCTION The present study was undertaken to investigate the known methods of determining liquid densities and to select one of them which might yield a better correlation if additional data were employed. An appropriate correlation would incorporate the following characteristics: (1) ease of handling in the computations, (2) results within engineering accuracy and (3) data easily accessible. After careful investigation of the different methods, the Standing-Katz correlation1 appeared to have these characteristics. The Standing-Katz correlation for liquid densities of hydrocarbon systems is based upon limited data. The original work, which employs apparent densities for methane and ethane, was based on data of 15 saturated crude oils in equilibrium with natural gas. Therefore, a test of the correlation using additional experimental data would serve either to validate the Standing-Katz correlation or to form a basis for corrections if such were necessary. PROCEDURE AND GRAPHICAL SOLUTION The data of 154 bottom-hole or recombined reservoir fluid samples were employed to calculate the densities at 14.7 psia and 60°F by the method described in Table 1. Volume contributions for methane and ethane were assigned from the pseudo liquid density plot (Fig. I), but volume contributions for N2, CO2 and HzS were not considered. The densities computed at reference conditions of 14.7 psia and 60°F were then elevated to their respective temperatures and pressures by the use of the correction curves proposed by Standing.' The percentage differences between the experimental values and the calculated values were determined; these samples were grouped into temperature, pressure and density ranges to determine a possible trend. It appeared that, at temperatures above 160°F, the calculated density values were consistently larger than the experimental density values; however, no obvious trend was found for the density and pressure ranges. There is a relationship between the pressure-correction and temperature-correction curves. This relationship, embodied in the temperature-correction curve, can be expressed as follows. The experimental densities of 125 samples above their bubble points were known. From these isothermal data, it was possible to determine an experimental-density difference between two pressures; and (from Fig. 2*), a calculated density difference between the same two pressures was determined. The relative effect of temperature must be considered in the calculated-density values. The effect of temperature can be determined by entering Fig. 3* at the two densities in consideration and obtaining the change in density at constant temperature. This change, or delta density, must be added

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

By G. R. Countryman, G. H. Thomas, I. Fatt

Miscible displacement in both thewetting and non-wetting phase has been studied in two-phase systems. Experimental data show that dispersion is a function of saturation and cannot be predicted from dispersion theory developed for single-phase miscible displacement. Some of the effect of saturation on dispersion is a result of the development of dead-end pores or a dendritic pore structure in multiphase systems. The presence of dead-end pores or dendritic structure is confirmed by pressure transient studies. INTRODUCTION Displacement by flooding with a miscible liquid is a possible means for recovering the estimated two-thirds of the oil that remains behind after primary production. The immense economic importance of a process that can recover such large quantities of oil has led to extensive laboratory studies of miscible displacement. As usual in production research, most laboratory studies of miscible displacement have not attemped to reproduce all of the conditions existing in a petroleum reservoir. In the early stages of investigation the need has been for information on the broad, basic principles of the phenomena involved in miscible displacement. The effects of overburden pressure, reservoir temperature, and wettability have been considered of secondary importance. However, one property of a petroleum reservoir which is expected to be of major importance and yet has been omitted from many laboratory studies is the presence of interstitial water. Two possible effects of interstitial water on the displacement mechanism in the hydrocarbon phase immediately come to mind. First, from the general1y accepted theory that capillarity governs the distribution of oil and water in porous rock one would expect that for water-wet rock the water will be in the small pores and oil in the large pores. A miscible displacement of oil carried out in the presence of water is operating in a pore size distribution shown in Fig. lb, whereas if the same test had been performed with only one phase present the pore size distribution is as shown in Fig. la. Although there is not yet a theory of miscible displacement which explains in detail the effect of pore size distribution, one would expect the differences between Figs. la and lb to influence the displacement efficiency. A second factor which may make a multiphase system different from a single-phase system is the presence in the multiphase system of dead-end pores or dendritic structure. Experiments of various kinds on reservoir rock have led to the belief that all pores in the network structure of a porous rock take part in conducting fluid during single phase fluid flow.1,2 There are then no dead-end pores and no fingers or dendritic structures containing stagnant fluid. In a multiphase system, however, the second phase may

By L. L. Handy

Solvent floods using slugs of solvent have been found to show continuity in behavior from the vapor pressure of the solvent to the critical pressure for the two-component driving gas-solvent system. In the pressure region between the solvent vapor pressure and the critical pressure for the gas-solvent system, the gar and solvent are only partially miscible. Although complete miscibility cannot be obtained at these pressures, complete oil recovery is possible in principle. In two-phase solvent floods the solvent is propagated tllrough the reservoir, primarily, in the vapor phase. The carrier gas requirements constitute a significant factor in the economics of the process. A qualitative theory is proposed for estimating the amount of dry gas required to move the solvent through the reservoir. The theory shows that for two-phase solvent floods the total gas needed is a minimum at the vapor pressure of the solvent and at the critical pressure for the gas-solvent system, and is a maximum at some intermediate pressure. The predictions of the theory are supported by experimental studies using methane, butane and decane or methane, propane and decane in a natural sandstone core. INTRODUCTION Previously, solvent slug processes have been found effective for oil recovery in two pressure ranges. First, conventional miscible displacements are possible at pressures greater than the critical pressure for the gas-solvent system. Second, Jenks, et al,' have shown that, at pressures slightly in excess of the vapor pressure of the solvent, a solvent slug can be propagated through a reservoir by a gas essentially insoluble in the liquid solvent. The solvent bank displaces the oil ahead of it. Both of these processes, at least ideally, are capable of recovering all of the oil in the swept regions. Slug processes for which the gas and solvent are partially miscible have not been considered; that is, those systems for which the solvent and driving gas form two equilibrium phases in which the vapor phase contains a significant amount of solvent and the liquid phase an appreciable amount of the driving gas. Welge and Johnson' have shown that the gas needed to movc a solvent slug through the reservoir increases with increasing pressure above the vapor pressure of the solvent. It will be shown that solvent slug processes can, theoretically, recover all of the oil at any pressure greater than the vapor pressure of the solvent. But the amount of gas required to move the solvent through the reservoir depends very much on the pressure and temperature. In the present study a maximum in the gas requirements was both predicted theoretically and observed experimentally. This result has not been reported previously, and would not have been predicted from the Welge and Johnson model. The gas requirements are a minimum at the pressures corresponding to the vapor pressure of the solvent and again at the critical pressure for the gas-solvent system, and are a maximum at some intermediate pressure. AN APPROXIMATE THEORY FOR TWO-PHASE SOLVENT FLOODING The differences and similarities between conventional solvent floods and two-phase solvent floods are best understood by referring to concepts developed for miscible displacement in which miscibility is generated in the reservoir. In Fig. 1(A), a ternary diagram is shown for a hypothetical gas-solvent-oil system. To be rigorous the three components should each consist of a single molecular species. The pressure for Fig. 1(A) is greater than the critical pressure for the binary gas-solvent system at the specified temperature. Diagrams of this type are the ones most frequently referred to in discussions of enriched-gas drive and miscible displacement. A limiting tie line is shown tangent to the two-phase envelope and intersecting the gas-solvent line at Point A. To obtain generated miscibility with this type system, others have shown that, for an oil of Composition D, a mixture of gas and solvent must be injected which is richer in solvent than that composition indicated by A. An oil repeatedly contacted with a gas phase richer than A changes toward a composition which would be at equilibrium with the injected mixture, that is, a composition lying on a tie line which passes through the injected-gas composition. Since no such tie line exists, the oil is enriched to the point at which it becomes directly miscible with the injected mixture. At pressures lower than the critical pressure for the gas-solvent system, other types of phase diagrams are observed. The ones of interest in this paper are for pressures greater than the vapor pressure of the solvent, but less than the critical pressure of the gas-solvent system. Such a ternary diagram is shown on Fig. 1(B). In this case, two-phase behavior is observed not only for gas-oil mixtures, but also for certain compositions in the gas-solvent system. If a gas of Composition A (a dew-point vapor) is injected, once again the original oil is enriched by successive

By A. K. Csazar, L. W. Holm

This paper describes our investigation of a post-water-flood, oil recovery process which consists of injecting a slug of propane followed by water. Also described are the results obtained by applying a modification of the process in which gas was injected ahead of the water. Under the conditions of the latter experiments, misci-bility was not achieved between the propane and gas. Preliminary experiments or) uniform, watered-out sandstone cores showed that an oil bank could be formed and produced by applying this recovery process. However, since reservoirs are not uniform in structure, the process also was applied to porous media containing irregular porosity and to stratified sand systems. As a supplenzerrt to the experinlental work, a mathernatical procedure was developed for calculating the performance of the recovery process in a bounded, layered, porous system with crossflow between layers. As a specific example, the method was applied to predict the perforrnance of the recovery process in a 6-ft long, two-layer, stratified, unconsolidated sand model for comparison with experinlental data. The calculations were programed for the ZBM 704 computer. The equations and calcula-tional procedure presented can be extended to systems containing any number of randomly distributed permeability variations or any number of parallel layers. INTRODUCTION The problem of recovering the oil that remains in a reservoir which has been waterflooded is receiving considerable attention now as an increasing number of water floods reach an economic limit. A large number of the waterflood projects are in shallow reservoirs which are at pressures below 1,000 psi. It has been demonstrated in the laboratory that post-waterflood oil can be recover-ered by miscible displacement, but the LPG-gas, miscible flood and the enriched gas drive cannot be applied effectively at pressures below 1,000 psi. Only a few reports have appeared in the literature2-4 on low pressure, partially miscible recovery methods. However, it is possible to use LPG in a partially miscible displacement process in a reservoir where pressures of 200 to 1,000 psi can be achieved. Under these Pressures and at normal reservoir temperatures, propane is miscible with the oil; but, of course, gas or water used to drive the propane slug would not be miscible with the propane. Because of the lack of complete miscibility, it has generally been concluded that excessive amounts of propane would be required to recover oil and that such a recovery method would not be economical; however, we have found that under conditions present in certain reservoirs, an imrniscible recovery process can be applied effectively. The oil saturation in reservoirs at the economic limit of waterflood projects is usually in the range of 20 to 35 per cent of the pore space." A certain portion of this oil is left trapped by water in various size pores of the rock, but a good part of this so-called "residual" oil can be present in the less permeable lenses or layers of the reservoir rock which were by-passed to some degree by the water. The oil in these permeability traps can be produced only if favorable pressure gradients are formed in the reservoirs between adjacent zones of high and low permeabilities. A low viscosity liquid, miscible with the oil in place, which is driven by water through a stratified or heterogeneous porous system can aid in the development of these favorable pressure gradients. The oil that is released thereby from the permeability traps can be recovered by the subsequent water flood. Studies were made to determine how much oil could be recovered from homogeneous and stratified cores and models, which had been water flooded, by injecting a slug of propane and driving it with water. The effect of injecting a slug of gas ahead of the water was also determined. Most of the work described herein was done with the propane-water combination; unless otherwise specified, no gas was injected. The principal objectives of the investigation were to determine (1) if an oil bank could be formed and (2) what ratio of oil recovered to propane injected would be obtained. A further objective was to develop a method for calculating fluid-flow performance in stratified systems which would account for fluid transfer between zones in hydrodynamic communication but of different permeabilities. THEORETICAL ANALYSIS In a theoretical study of the recovery process, analytical expressions were derived to calculate the pressure distribution, the fluid flux in longitudinal (parallel to layers) and transversal (across the layers) directions, and the fluid distribution at any point in the system. The equations were developed for a two-layer porous system in which it was assumed that the fluids in the system were incompressible and that capillary and gravity effects were

By M. Mortada

Nonsteady-state flow of slightly compressible liquids in porous media of non-uniform properties has been the subject of a number of recent studies. Most of these studies considered one-dimensional flow in systems with cylindrical geometry and with the properties of the porous medium varying with radial distance only. A. Houpeurt' studied the problem in its general form by examining the solutions of the differential equation, In Eq. 1, the permeability k and the porosity + are arbitrary functions of radial distance. Houpeurt concluded that, except for special forms of the functions k and the problem is not well suited for mathematical treatment, but lends itself to high-speed digital computation. Carslaw and Jaeger' discussed the solution of Eq. 1 for the case in which the permeability varies as a power of the radial distance. P. Albert3 studied analytically the pseudo steady-state solutions in systems of finite radial extent and non-uniform properties. He also obtained some numerical solutions for nonsteady-state flow using a high-speed computer. William Hurst4 examined the interference pressure drop due to a point sink located in an infinite system consisting of two permeability regions in series. In this case, the permeability changed at a radial distance r from the point sink. INTERFERENCE BETWEEN OIL FIELDS IN A NON-UNIFORM AQUIFER The interference pressure drop for oil fields located in an extensive aquifer of uniform properties is discussed in a number of places in the literature. This work considers the interference pressure drop for oil fields located in a non-uniform extensive aquifer comprising two regions of different properties as shown in Fig. 1. Region I of the aquifer extends from which is the radius of Oil Field A, to ar,,, where a>1. It has permeability k,, porosity +, and effective compressibility c,. Oil Field B is located in Region II which extends from ar, to infinity and has different properties from those of Region I, but the same formation thickness h. Region II has permeability k2, porosity and effective compressibility c2. The choice of this type of aquifer should not imply that it is of frequent occurrence. It is chosen primarily to illustrate the effects of non-uniform aquifer properties on oilfield interference. It is of interest to point out that a non-uniform radial aquifer in which both the permeability and the porosity are inversely proportional to the radial distance can be treated in a manner similar to a linear aquifer of uniform properties. The situation where only the permeability varies as the inverse of the radial distance is also amenable to analytical treatment. Both aquifers behave quite differently from a radial aquifer with uniform properties so far as the water influx and the interference pressure drop are concerned. METHOD OF SOLUTION The expression for the interference pressure drop in Oil Field B due to a constant rate of water influx in Oil Field A is developed in the Appendix. This is accomplished by the simultaneous solution of two partial differential equations describing the pressure in Regions I and 11. The two differential equations are solved subject to the following conditions: (1) at r,, the pressure gradient is constant; (2) at ar the pressure and

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 D. D. Dunlop, J. R. Welker

The growing interest in the use of CO, in crude oil recovery increases the need for data on the effect of CO, on hydrocarbon physical properties. Data are presented on the solubility of CO, in various dead oils, the swelling changes in CO2-oil solutions and the effect of CO, on dead oil viscosity. This latter property shows the most pronounced effect, with viscosity reductions up to 98 per cent of the uncarbonated viscosities. An empirical method of estimating the viscosity of carbonated oils is presented. The apparatus and procedures used are described in sufficient detail to allow others to make similar studies. INTRODUCTION The effect of dissolved carbon dioxide on the swelling and viscosity reduction of specific hydrocarbon oils has been observed and recorded by a number of investigators.'- me object of this paper is to offer a means of predicting these effects for crude oils free from natural gas, using the dead state viscosity and gravity of the crude oils. The CO, solubility and swelling of numerous crude oils were determined in a visual cell at various pressure levels. The viscosity of the oils carbonated to various pressure levels was then determined by measuring the pressure drop across a capillary tube. From these data, the physical properties were correlated empirically. The resulting correlations allow the prediction of CO, solubility, swelling and viscosity reduction if the dead state gravity and viscosity of the oils are known. SOLUBILITY AND SWELLING MEASUREMENT EQUIPMENT AND PROCEDURE A high pressure visual cell was installed in a constant temperature cabinet. A test gauge was attached at the top of the cell for pressure measurement, and a line was run through the cabinet wall to a wet test meter which was used for volumetric measurement of the gas. The first step in making a test run was to put the oil in the cell up to a level about half to two-thirds of the total volume. This required about 50 to 65 ml of oil. carbon dioxide was then bubbled up through the oil for a time during which the pressure of CO2 in the cell was kept above 800 psia. Saturation of the oil with CO2 at this pressure and ambient temperature was confirmed by slowly bleeding CO2 through a valve to the atmosphere. If the oil was completely saturated with CO2, bubbles of gas would form in the oil at the first small decrease in pressure. If the oil was under-saturated, no bubbles formed until the pressure was decreased to the saturation pressure existing in the oil. If this saturation pressure was lower than that desired, more CO2 was bubbled through the oil until the desired level was reached. After saturation at ambient temperature was completed, the cabinet temperature was adjusted to the desired level and the cell was allowed to reach temperature equilibrium. After temperature equilibrium was reached, the pressure was again decreased slightly, and the oil again checked for full CO2 saturation at the cell pressure. The pressure now had changed because of the difference in solubility of the CO, in the oil at higher temperatures and the expansion of CO2 as the temperature increased. The outlet tube from the cell was then connected to the wet test meter and the CO2 was allowed to flow slowly out of the cell and through the wet test meter at ambient temperature and pressure. The water in the wet test meter had previously been saturated with CO2 at ambient temperature and pressure by allowing CO2 to flow continuously through it lor a period of several hours. The gas flow was stopped at several pressures during the run and the cell was allowed to come to equilibrium; this made possible the measurement of solubility and swelling data at the intermediate pressures. The volume of the oil in the cell was recorded at each of the equilibrium pressures in order to obtain swelling data. DATA AND RESULTS The solubility of CO2 in the oil was calculated by the relationship V — V, where R. = solubility of CO, in crude oil, cubic feet of CO, measured at 60F and 1.0 atm/ bbl of dead state oil at the temperature under which solubility was measured, V, = volume of gas released from the cell between the saturation pressure and zero pressure, corrected to 60F and 1.0 atm, cu ft, V, = volume of CO: contained in the gas space above the oil, corrected to 60F and 1.0 atm, cu ft, and V, = volume of the dead oil in the cell in bbl at the temperature of the run. The volumetric data of Sage and Lacey' were used to calculate V., from the volume of CO2 at high pressures. The swelling factor was calculated as where V, is the volume of the C0,-

By K. E. Brown, G. H. Fancher

An 8,000-ft experimental field well was utilized to conduct flowing pressure gradient tests under conditions of continuous, multiphase flow through 2 3/8-in. OD tubing. The well was equipped with 10 gas-lift valves and 10 Maihak electronic pressure recorders, as well as instruments to accurately measure the surface pressure, temperature, volume of injected gas and fluid production. These tests were conducted for flow rates ranging from 75 to 936 B/D at various gas-liquid ratios from 105 to 9,433 scf/bbl. An expanding-orifice gas-lift valve allowed each flow rate to be produced with a range of controlled gas-liquid ratios. From these data an accurate pressure traverse has been constructed for various flow rates and for various gas-liquid ratios. A comparison of these tests to Poettmann and Carpenter's correlation indicates that deviations occur for certain ranges of flow rates and gas-liquid ratios. Numerous curves are presented illustrating the comparison of this correlation with the field data. Poettmann and Carpenter's correlation deviates some for low flow rates and, in particular, for gas-liquid ratios in excess of 3,000 scf/bbl. These deviations me believed to be mainly due to the friction-factor correlation. However, Poettmann and Carpenter's correlation gives excellent agreement in those ranges of higher density. This was as expected and predicted by Poettmann. He pointed out that their method was not intended to be extended to those ranges of low densities whereby an extreme reversal in curvature occurs.l As a result of these experimental tests, correlations using Poettmann and Carpenter's method were established between the friction factors and mass flow rates which are applicable for all gas-liquid ratios and flow rates. Definite changing flow patterns do not allow any one correlation to be accurate for all ranges of flow. INTRODUCTION The ability to analytically predict the pressure at any point in a flow string is essential in determining optimum production string dimensions and in the design of gas-lift installations. This information is also invaluable in predicting bottom-hole pressures in flowing wells. Although this problem is not new to industry, it has by no means been solved completely for all types of flow conditions. Versluys,2 Uren, et al, 3 Gosline,4 May,5 and Moore, et al,6,7 were all early investigators of multiphase flow through vertical conduits. However, all of these investigations and proposed methods were very limited as to their range of application. Likewise, many are extremely complicated and therefore not very useful in the field. Only in the last decade have any significant methods been proposed which are generally applicable. The most widely accepted procedure in industry at the present time is a semi-empirical method developed from an energy balance, proposed by Poettmann and Carpenter8 in 1952. Their correlation is based on actual pressure measurements from field wells. Accurate predictions from this correlation are limited to high flow rates and low gas-liquid ratios. Although this method will be discussed in detail later, it should be pointed out that two important parameters, namely the gas-liquid ratio and the viscosity, were omitted in their correlation. The viscosity was justifiably omitted since their data was in the highly turbulent flow region for both phases, and most wells fall in this category. The gas-liquid ratio was incorporated to some extent in the gas-density term. In 1954, Gilbert9 presented numerous pressure gradient curves obtained from field data for various flow rates and gas-liquid ratios for the determination of optimum flow strings. However, no method is presented for predicting pressure gradients except by comparison to these curves. Baxendell, et al,10 proposed a method based on Poettmann and Carpenter's procedure

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