Search Documents

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 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 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 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 A. L. Barnes

Residual oil in watered-out reservoirs is a tremendous reserve which has been unrecoverable by established production methods. A study of the new recovery methods indicated that the forward combustion urocess might- recover oil from such reservoirs; however, no thermal recovery operating experience in a watered-out system was available. The Delaware-Childers pilot thermal test was undertaken to test the feasibility of thermal recovery in this watered-out reservoir. The pilot test consisted of a 2.22-acre, inverted five-spot in the 600-ft deep Bartlesville sand. The reservoir in the pilot area had a porosity of 20.6 per cent, a permeability of 118 md and an average sand thickness of 45 ft. The reservoir contains a 33' API, 6-cp oil. Combustion was started Nov. 22, 1960. The initial air injection capacity was 750 Mscf/D, but it was eventually increased to 2,000 Mscf/D. The test was surrounded by an active water flood; therefore, water production was initially high, but decreased as the heat wave moved toward the producing wells. Oil-bank arrival at an individual well was indicated by a drop in GOR and WOR, and an increase in oil production. Combustion-front arrival was evident at three of the pilot producers, and they were plugged. Cumulative oil production from pilot area wells war over 12,000 bbl. Operational difficulties were negligible and only conventional equipment was necessary. The combustion efficiency of this test averaged over 80 per cent. Results from coring showed that the leading edge of the combustion front tended to be wedge-shaped but a nearly complete sweep of the reservoir was eventually obtained. An isopach map based on evidence from 10 core holes and the existing wells showed that 126 acre-ft had been swept by the heat wave. Using this swept volume, an air requirement of 15.7 MMscf/acre-ft was calculated. It was calculated that 275 ,STB/acre-ft was consumed bv the heat wave. INTRODUCTION There is a large amount of oil remaining in reservoirs that have been water flooded. A study of ways to recover this oil showed that the forward combustion process might be applicable. Results from a number of forward combustion tests have been reported in the literature, but none of these tests were conducted in a watered-out system. The Delaware-Childers pilot thermal test was initiated in 1960 to define the operating characteristics of underground combustion in this watered-out Bartlesville sand reservoir, The purpose of this paper is to present the pilot test results in detail. FORWARD COMBUSTION PROCESS The forward combustion process consists of initiating combustion in the formation surrounding an injection well and driving this heat wave through the formation toward offset producing wells. As the combustion front progresses through the reservoir, oil and formation water are vaporized, driven forward in the gaseous phase, and recon-densed in the cooler part of the formation. These distilled liquids, water of combustion and gaseous combustion products, form a bank or three-phase region ahead of the burning front. This bank pushes mobile reservoir fluids toward the production wells. The rate of movement of the combustion front is controlled by the rate at which the nondistill able residue which serves as process fuel can be completely burned off the sand. The production performance from a heat wave conducted in a watered-out oil reservoir will differ from one conducted in a dissolved gas-depleted reservoir because of the difference in fluid saturations. The primary depleted reservoir contains connate water saturation, relatively high oil saturation and some gas saturation. The watered-out reservoir contains a highly mobile water saturation, residual oil saturation and little gas saturation. As the combustion front moves forward in either type of depleted reservoir, it drives the distillable oil, the water in place and the water of combustion ahead, and burns the nondistillable portion of the oil. The heat-wave process in a watered-out reservoir differs from the process in a

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 P. J. Root, J. E. Warren

An idealized model has been developed for the purpose of studying the characteristic behavior of a permeable medium which contains regions which contribute significantly to the pore volume of the system but contribute negligibly to the flow capacity; e.g., a naturally fractured or vugular reservoir. Unsteady-state flow in this model reservoir has been investigated analytically. The pressure build-up performance has been examined in some detail; and, a technique for analyzing the build-up data to evaluate the desired parameters has been suggested. The use of this approach in the interpretation of field data has been discussed. As a result of this study, the following general conclusions can be drawn: 1. Two parameters are sufficient to characterize the deviation of the behavior of a medium with "double porositym from that of a homogeneously porous medium. 2. These parameters can be evaluated by the proper analysis of pressure build-up data obtained from adequately designed tests. 3. Since the build-up curve associated with this type of porous system is similar to that obtained from a stratified reservoir, an unambiguous interpretation is not possible without additional information. 4. Differencing methods which utilize pressure data from the final stages of a build-up test should be used with extreme caution. INTRODUCTION In order to plan a sound exploitation program or a successful secondary-recovery project, sufficient reliable information concerning the nature of the reservoir-fluid system must be available. Since it is evident that an adequate description of the reservoir rock is necessary if this condition is to be fulfilled, the present investigation was undertaken for the purpose of improving the fluid-flow characterization, based on normally available data, of a particular porous medium. DISCUSSION OF THE PROBLEM For many years it was widely assumed that, for the purpose of making engineering studies, two param- eters were sufficient to describe the single-phase flow properties of a producing formation, i.e., the absolute permeability and the effective porosity. It later became evident that the concept of directional permeability was of more than academic interest; consequently, the degree of permeability anisotropy and the orientation of the principal axes of permeability were accepted as basic parameters governing reservoir performance. 1,2 More recently, 3-6 it was recognized that at least one additional parameter was required to depict the behavior of a porous system containing regions which contributed significantly to the pore volume but contributed negligibly to the flow capacity. Microscopically, these regions could be "dead-end" or "storage" pores or, macroscopi-cally, they could be discrete volumes of low-permeability matrix rock combined with natural fissures in a reservoir. It is obvious that some provision for the inclusion of all the indicated parameters, as well as their spatial variations, must be made if a truly useful, conceptual model of a reservoir is to be developed. A dichotomy of the internal voids of reservoir rocks has been suggested.7,8 These two classes of porosity can be described as follows: a. Primary porosity is intergranular and controlled by deposition and lithification. It is highly interconnected and usually can be correlated with permeability since it is largely dependent on the geometry, size distribution and spatial distribution of the grains. The void systems of sands, sandstones and oolitic limestones are typical of this type. b. Secondary porosity is foramenular and is controlled by fracturing, jointing and/or solution in circulating water although it may be modified by infilling as a result of precipitation. It is not highly interconnected and usually cannot be correlated with permeability. Solution channels or vugular voids developed during weathering or burial in sedimentary basins are indigenous to carbonate rocks such as limestones or dolomites. Joints or fissures which occur in massive, extensive formations composed of shale, siltstone, schist, limestone or dolomite are generally vertical, and they are ascribed to tensional failure during mechanical deformation (the permeability associated with this type of void system is often anisotropic). Shrinkage cracks are the result

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 V. J. Driscoll

A method is given for checking pore volume reserve estimates and transient water inflzix factors utilizing early well performance. Principal applications are based on the observation that there is a convenient transient time-range in terms of well spacing in which observed pressure drop will give a good indication of fluid-in-place. Primary quantity determined is a porosity-feet value applicable in the vicinity of an observation well. Checks on porosity or net pay thickness and permeability values also are obtained. These are particularly valuable in improving quality of interpretive factors entering into pore-volume reserve determinations for larger fields. In smaller fields, direct determination of total field reserves and approximate location of boundaries may be made in some cases. Well interference and production data plus a build-up curve on the observation well are needed. The build-up curve may be taken at any time. Field applications are given. The method can be applied with the completion of the second well on a new discovery or to older reservoirs by major changes in withdrawal rates from a local area. Appropriate test length permits results to be obtained with equal accuracy on widely spaced wells. Criteria for optimizing test conditions are given. Provision is made for various boundary and well effects including incomplete observation well build-up. INTRODUCTION Due to low quality data, initial pore volume reserve and water influx estimates may often vary by a factor of two or more from results obtained at a later date. Unfortunately, by conventional methods, more reliable answers from performance usually are not available until after most significant economic decisions have been made. Furthermore, in older reservoirs, uncertainties may still exist as to original oil-in-place, recovery factors and residual oil saturations due to incomplete or inaccurate performance data. This paper outlines situations and appropriate arrangements of basic transient theory whereby well interference data can be utilized to improve the quality of estimates under certain conditions. In addition to production data on offset producers, key information required is a measured pressure drop at an observation well, corrected for any incomplete build-up due to previous production, plus a formation capacity value derived from a conventional pressure build-up plot on the observation well. The method, as outlined, is strictly applicable only to a horizontal single-layered system of uniform rock characteristics containing a homogeneous fluid of low and constant compressibility. However, under the set of test conditions devised, requirements can usually be relaxed so that it should be applicable to most low permeability gas, water and undersaturated oil reservoirs. This is an area where more reliable information is certainly needed. BASIC CONSIDERATIONS A, B, C and F as subscripts or in the text refer to producing Well A, observation Well B, a more remote producing Well C and field boundary F. Combinations such as AB refer to distances and dimensionless times in terms of distances between these points, while single letters refer to distances and dimensionless times in terms of the well-bore or field radius. Equations covering definitions of quantities and complex situations are given in the Appendix. Key equations in practical units are given in Table 2. Transient Pressure Behavior Fig. 1 illustrates pressure drop at Point B with time due to initiation of a constant rate of withdrawal from Point A. This is shown as since a quantitative under-

By C. R. Johnson, R. A. Greenkorn, L. K. Shallenbarger

This paper describes a study, based on core data, of the directional permeability of a sandstone reservoir. Directional air permeabilities are explained and correlated with lithology by the tensor theory of permeability, which is extended to the more general case of heterogeneous anisotzopic porous media. In this case, the permeability tensor is made up of two components: (I) anisotropy (variation around a point) which correlates with bedding, and (21 heterogeneity (variation from point-to-point) which correlates with grain size. INTRODUCTION Previous studies of directional permeabilityl.5-8 have concerned themselves only with anisotropy at a point, rather than with both point-to-point differences and local anisotropies. In this paper, we present the existing tensor theory of permeability for anisotropic porous media and then extend it to the more general case of heterogeneous anisotropic porous media. The laboratory measurements of directional permeability are explained in view of this extended theory and the data are discussed in terms of the lithologic factors that correlate with it. The results show what the heterogeneity and anisotropy of the reservoir element are, and that the directional permeability correlates with lithology. Data used in this paper are measurements of air permeability in eight directions, spaced at 45° intervals, on 142 2-in. vertical plugs from 30 cores. The core data are meaningful in terms of directional permeability because the cores were oriented to within 45O during drilling and coring.3 Although we initially thought that only 60 per cent of the core material was reliable, subsequent study showed that almost all of it was reliably oriented. After determining the air permeabilities, the data were reduced to three independent variables for each plug: the major and minor permeability axes, and the direction of the major axis. These were obtained by converting the permeability data to the reciprocal square root of permeability and fitting the transformed data with ellipses according to the tensor theory of permeability.8 The point-to-point areal variation of the minimum permeability axis is related to grain size. The direction of the permeability axes, where the major axis exceeds the minor axis by at least 5 per cent (measured variation is about 4 per cent), correlates with the bedding. The permeability tensor used in this study must be considered as the sum of a scalar and a tensor, with the scalar being the minor axis permeability as a function of position, and the tensor the directional effect, which is additive permeability over the minor axis. In this case, the point-to-point variation or heterogeneity (minor axis) is substantially larger than the variation at a point or anisotropy. It may be important that this separation of "directional permeabilities" be recognized when considering migration of fluids due to permeability variation. Local migration may be due to anisotropy and point-to-point migration may be due to heterogeneity, but the direction and magnitude of these may not be the same. Furthermore, in truly heterogeneous systems, one would expect that anisotropy would be the smaller of the two effects. THEORY TENSOR THEORY OF PERMEABILITY FOR ANISOTROPIC POROUS MEDIA8 Darcy's law for flow in porous media in its usual form is where qis the flow rate vector, k is the permeability, is the viscosity, is the vector differential operator, and P is the Pressure. Consider flow, q, through a unit area in the direction of the unit vector n from elementary flow tubes parallel

Jan 1, 1965

By D. M. Waldorf

Steam permeability measurements have been made in the laboratory on several samples of natural reservoir materials. The steam temperatures and pressures were selected to simulate conditions which might exist in a reservoir during the injection of steam. For each sample tested, the experimental permeability to superheated steam was comparable to that measured with air and no evidence of plugging was detected. Some samples were exposed to water at various temperatures and plugging was found to occur in materials which contained significant quantities of monmorillonite clay. Temperature had little effect on the degree of plug-ning between 75 and 325 F. The measured pemeabilities tended to increase slightly with temperature, but the changes were small compared with the initial loss of per~neability on wetting. Sequential pemzeability measurements were made on two samples using air, water, steam, water and air, in that order. Both samples were water-sensitive and plugged extensively after the initial injection of water. Upon exposure to superheated steatm the samples dehydrated and their permenbilities to superheated steam were comparable to those initially measured with air. The remaining measuretnetzts with water and air confirmed that the water plugging was reversible and that the samples were not seriorrsly damaged during the tests. INTRODUCTION The swelling of water-sensitive clays during water floods has long been recognized as a potential source of reservoir damage. The recent extensive application of steam injection and stimulation has compounded this problem since both hot water and steam (as well as fresh water at reservoir temperatures) are, at sume time, in contact with the producing zone adjacent to the bore of a steam injection well. The purpose of this paper is to present data which compare the sensitivity of some natural sedimentary rock samples to water at various temperatures, and to super-heated steam. Some properties of montmorillonite clay are briefly reviewed, and comparisons are drawn between empirical data and the predicted behavior of the montmorillonite known to be present in the samples. PROPERTIES OF MONTMORILLONIT E CLAY Water initially adsorbs on dry N a -montmorillonite clay in discrete layers in the interlaminar space between clal platelets. The platelet spacing, which is 9.6 A (angstroms) for a dehydrated clay, has been observed to expand in discrete steps to 12.4, 15.5, 18.4 and 21.4 A spacings, indicating the formation of four discrete layers of regularly oriented water molecules.' The first two layers are easily formed by hydrating a dry sample to equilibrium in an atmosphere with carefully controlled humidity. The formation of the higher layers is more difficult. The usual X-ray diffraction patterns of the more highly hydrated samples indicate a gradual increase in the average spacing betwcen 15.5 and 19.2 A, followed by a discontinuous expansion to 31 A when the weight ratio of water to dry clay is between 0.5 and 1.2.' Platelet expansion above 31 A proceeds monotonically as the moisture is increased and no regular arrangement of the platelets ib observed. Water-sensitivity in sedimentary rocks is usually associated with Na-montmorillonite clay when it is in the noncrystal-line state. Mering3 found that the average lattice spacing of sodium montmorillonite hydrated at 68 F and 70 per cent relative humidity was 15.5 A, and that the spacing, at 92 per cent humidity was 16.5 A. The water adsorbed at the higher humidity has the same free energy as liquid water at 65.6 F. Kolaian and Low' used a tensiometer to measure the thermodynamic properties of water in diffuse suspensions of montmorillonite clays relative to pure water. They observed that water in suspensions as dilute as 6 per cent clay became partially oriented when left undisturbed. The bonding associated with this orientation was not extensive because the free energy difference between the water in suspension and pure water was only a few millicalories per mole. They also found that the measured free energy difference decreased rapidly with temperature and became negligible above 100 F. This evidence indicates that montmorillonites contained in sedimentary rocks would dehydrate to a crystalline structure when exposed to superheated steam, and that the rock permeability measured with steam would be equivalent to that measured with air. The effect of elevated temperatures on the swelline of montmorillonite clays in aqueous suspensions has not been investigated. The Gouy-Chapman diffuse-ion-layer theory has been used to predict the swelling pressure of clay suspensions in dilute salt solutions at room temperature with reasonable success. theory also correctly predicts the direction of the thermal response of Na-mont-morillonite swelling pressures in dilute salt suspensions, 9 Over the temperature range of 33 to 68 F, an increase in

Jan 1, 1966

By D. J. Pye

Certain high molecular weight synthetic polymers in very dilute solutions decrease water mobility in porous media 5 to 20 times more than would be expected from the solution viscosity. This indicates that increasing water viscosity to reduce adverse mobility ratios in water floods is economically feasible in many situations. Laboratory results are consistent with accepted theory with respect to improved areal sweep, permeability distribution and displacement efficiencies as would be predicted from the solution mobility measured in the core. The laboratory results have been confirmed by field pilot floods which recovered 80 per cent and 100 per cent more oil than comparative water-flood pilots in fields containing 16 cp and 130 cp oil respectively. Economic analyses based on pilot data indicate the method to be profitable. INTRODUCTION The desirability of improving the water-oil mobility ratio in waterflood-ing operations has long been appreciated, especially for high viscosity crudes. Methods of reducing the oil viscosity by application of heat, gas resaturation and miscible drives have met with some success. Conversely, it has been generally conceded that increasing the water viscosity would also be effective, but the amount of thickening agent required to effect an appreciable increase in viscosity is discouragingly high. The use of such materials as glycerin, sugar, or glycols is completely out of the question economically, so that considerable attention has been given to the use of the much more efficient synthetic water-soluble polymers as viscosity improvers; but even these materials are not economically attractive. If increased water viscosity could be economically realized, marked improvements in areal sweep efficiency, as discussed by Caudle and Wittel and others, would be realized in field operations. Similarly, the improved mobility ratio would bring about increased recoveries by correcting permeability distribution problems as discussed by Stiles2 and Dykstra and Parsons,' and improvements in displacement efficiency would be obtained as recognized by Buckley and Leverett4 and Welge.5 LABORATORY STUDIES THE RESISTANCE EFFECT In the course of research to develop a more efficient water-soluble polymer viscosity builder, we discovered that a very few types of water-soluble polymers exhibited a most unusual and very interesting property not previously observed. Among the earliest materials showing such activity were polymers containing acry-lamide. Fig. 1 shows a typical plot of the solution viscosity vs concentration for one of these polymer solutions as determined in an Ostwald viscometer. As usual, such a plot is an approximate straight line on semi-logarithmic paper. Now, an inspection of Darcy's law, shows that a viscosity value can also be determined in a formation sample if the pressure, flow rate and permeability are known. If the core permeability is previously determined with wa- ter or brine, then the viscosity of another aqueous fluid may be determined from the equation. The viscosities so determined for glycerin (glycol or polyvinyl alcohol solutions, for example) agree with the viscometer determinations. However, in the case of the water-soluble polymer solutions in which we are interested, the viscosities measured in the formation sample depart very markedly from the viscometer values. A typical relationship for a polymer solution is also shown in Fig. 1. This unusual departure from the expected response is referred to as a resistance property of the polymer and is quantified as "resistance factor". Since the resistance factor R can best be determined in a core under study, the term has been defined on the basis of the ratio of the brine mobility to the polymer solution mobility at residual oil saturation so that:

Jan 1, 1965

By R. L. Dalton, L. A. Rapoport, C. W. Carpenter

Procedures have been developed to study the performance characteristics of unconfined pilot water floods using scaled laboratory models. The effects of operating conditions on pilot performance for variozrs well patterns have been investigated. The limitations of the laboratory model studies are outlined. Relationships between the oil recoveries of pilots and large-scale floods have been established. The applications of model studies to the prediction of large-scale waterflood performance and pilot flood design are discussed. INTRODUCTION Pilot water floods are undertaken to evaluate operational procedures and to gain advance information about the performance of an extensive water flood. In general, however, the performance of a pilot flood is not directly indicative of what can be expected from a large-scale pattern water flood. This difference stems from the fact that the over-all flow configurations are not the same for the two systems. A water flood normally consists of an array of identical well patterns. In such an array, the perimeters of the well patterns are axes of symmetry and act as impermeable boundaries. Thus, an extensive, pattern water flood can be visualized as a repetition of "confined" floods. In contrast, a pilot water flood involving only one or a few well patterns is "unconfined". Here the well patterns are not balanced by other flood units; hence, the perimeter of the pilot area does not act as an effective boundary. Accordingly, only a portion of the fluids injected (and also of those originally contained within the pilot area) is captured by the producers of the pilot. The rest escapes into the surrounding reservoir. This usually causes the amounts of oil and water produced in the pilot to be different from those recovered by a producer in an extensive water flood. Moreover, the recovery observed in the pilot can be expected to be greatly influenced by operating conditions, i.e., by the bottom-hole pressures at the injection and producing wells. The above considerations show that full utilization of pilot test results and reliable prediction of large- scale water floods require a knowledge of the relationship between the recovery performance of confined and unconfined well patterns. To date, only limited basic information is available on the behavior of pilot water floods. Craig, et al,' have investigated the effects of initial gas saturation, oil-to-water mobility ratio, and well damage on the oil recovered from a single, "normal" five-spot pilot pattern (Fig. 1). Rosenbaum and Matthews2 studied the effect of initial gas saturation and mobility ratio on the ratio of production rate to injection rate for various groupings of five-spot patterns. Paulse113 investigated the effect of mobility ratio and distance from the pilot to the edge of the reservoir on the areal sweep in a single "inverted" five-spot (Fig. 1). However, none of these studies considered the effect of operating conditions on pilot flood behavior. Moreover, hardly any of the work to date dealt explicitly with the oil-recovery performance of the pilot floods. The purpose of this paper is to introduce the concept of interrelation between pilot flood behavior and operating conditions, to illustrate this interrelation by means of experimental data and to compare the recovery performances of dif-

By A. L. Jr. Pozzi, A. O. Garder, D. W. Peaceman

A new numerical method is proposed for the solution of multidimensional miscible displacement problems. Besides the usual stationary grid associated with numerical procedures, the method uses moving points for which positions and concentrations are computed each time step. Based on the method of characteristics for treating combined transport and dispersion, the method applies equally well to any number of dimensions and, in contrast with other numerical procedures, properly takes into account the physical dispersion, no matter how small. Extensive tests show that accurate solutions of one-dimensional problems can be obtained over a wide range of values of the dispersion coefficient, including zero. They also show that the moving points do not need to be uniformly spaced, and that increasing the number of moving points beyond 2/grid interval does not significantly improve the accuracy of the solution. Two-dimmsional calculations using the new method were carried out to simulate laboratory model displacements in which gravity caused overriding of solvent. Good agreement between observed and computed profiles and measured and computed concentrations of produced fluid was obtained for a viscosity ratio of 1.89; at higher viscosity ratios, fair agreement was obtained Computed results were found to be significantly affected by vertical dispersion, and negligibly affected by dispersion in the direction of flow, in agreement with conclusions based on experimental observations. INTRODUCTION Although miscible displacement processes are potentially of great economic significance, limited progress has been made in developing mathematical methods for predicting the outcome of a solvent flood. Because of the complexity of the partial differential equations which describe multidimensional miscible displacement, numerical methods are a natural approach to their solution. The task of finding suitable numerical approximations to one of the partial differential equations, which involves both dispersion and convective transport, has been particularly difficult. Dispersion in the absence of transport is described by the heat flow equation, a second-order equation of parabolic type. This equation has been very successfully treated by numerical methods. Transport in the absence of dispersion is described by a first-order equation of hyperbolic type. This equation has been treated numerically with some success in one dimension, both by Lagrangian and Eulerian techniques, 1.2 but extension to two or more dimensions has not been satisfactorily accomplished. Usually one of two things happens: either the numerical solution develops oscillations or it becomes smeared by artificial dispersion resulting from the numerical process. When transport and dispersion are considered simultaneously, this numerical dispersion may dominate the low physical dispersivicy that generally characterizes miscible displacement. Stone and Brian2 have developed a new difference equation for treating combined transport and dispersion in one dimension. By suitable choice of certain parameters, they succeeded in markedly reducing both oscillation and numerical dispersion. However, a satisfactory extension to higher dimensions has not yet been found. Peaceman and Rachford3 presented a centered-in-time and centered-in-distance difference equation combined with a 1,transfer of overshoot" procedure which was demonstrated to work well in one dimension. The technique was applied to the two-dimensional miscible displacement problem, but subsequent calculations with zero dispersivity showed little change in results compared with the calculations they presented. This indicates that their method involved a numerical dispersion of the same order as the physical dispersion. The present paper describes an alternate method for solving the combined transport-dispersion

Jan 1, 1965

By V. S. Gupt, W. H. Somerton

Rock may undergo great changes in physical properties when heated to high temperatures and then cooled, The temperature and intensity of reactions causing rock alterntiorl.s can he controlled by introducing certain chemicals during heat treament. Three typical outcrop sandsone samples were saturuted with common salt solutions, then heated to several maximum temperatures. After cooling, it was found that per-rrreuhilities had increased much more for salt-saturated samples than for samples not saturated with salt but heated to the same temperatures. This was only true, however. for samples heated above the melting points of the particrrlar salts. Potassium chloritle was particularly effective with Bandera sandstone. Samples saturated with potassium chlorirlr. solution and heated to 900C showed an 11-fold increa.se in perrrieahility. Samples without potasrirdrn chloritle but heated to the .same temperature. showed only a 2.7-jold incresre in permecahility. In application, it seems possible that injecting chemicalc into the formations from a wellbore followed by applying intensive borehole heating might promote reactions which would greatly improve permeability of the formalions. INTRODUCTION Great changes in the physical properties of rocks which have been heated to temperatures in the range of 600 to 800C have been reported earlier. Permeability increases of 50 per cent, and equivalent decreases in sonic velocity and breaking strength have been observed. Although there might be some suppression of certain reactions responsiblc for these changes when rock samples are heated under simulated reservoir pressure conditions. recent work has shown that these are more than offset by the increased importance of other reactions at high pressures. The reactions considered responsible for alteration of rock properties by heating include differential thermal expansion. dehydration, phase changes and dissociation of mineral constituents. Ceramists have long known that temperatures at which reactions occur can be raised or lowered by the presence of certain impurities in the sys. tcnl. Thus the possibility exists that by saturating rocks with appropriate solutions. desired reactions might be accelerated and unwanted reactions might be suppressec! when the rock is heated. Purpose of the present work was to investigate the effects of several common salts on the thermal alteration of sandstones. Types of reactions which might occur are reviewed and the findings of a number of thernlochemical alteration tests are presented. THEKMOCHEMICAL ALTERATION OF SANDSTONE MINERALS The most abundant constituent of most sandstones is quartz. Quartz can exist in at least eight different forms, but in the temperature range of immediate interest (to 900C) only two forms are of importance—alpha quartz below 573C and beta quartz above this temperature. Although the alpha-beta inversion of quartz is very important in thermal alteration of rocks. it is not particularly important here because the inversion temperature cannot be changed significantly by the presence of impurities or by the application of pressure. The transition temperature of quartz to tridynlite is 867C. but the transition is very sluggish. Tridymite is considered to be a metastable transition phase between two stable phases: quartz?cristobalite. However. the transition temperature for tridymite to cristobalite is 1,470C.4 The quartz-tridymite-cristobalite transition can be accelerated substantially by the presence of fluxing agents. Finely divided calcium, magnesium and titanium oxides accelerate the conversion while A12O may retard it. A Mixture of NaCl and CsCl reduces the temperature for cristobalite formation by as much as 420°C. Feldspars constitute the second most important mineral in sandstones. The melting points of feldspars range between 1,000 and 1,700°C depending upon the variety. No information available indicated any important role of fluxing agents in the thermal alteration of this group of minerals. Carbonate minerals are subject to dissociation at temperatures within the range of the present investigation. The dissociation of magnesite can start at a temperature of 373°C, but the reaction is sluggish and might not occur untll a temperature of 500°C or ,higher is reached. Dolomite dissociates in two stages at 500°C and 890°C, whereas calcite dissociation temperature is about 885C. Because CO2 is released in the dissociation of carbonates, ail such reactions are somewhat dependent upon CO2 partial pressure. In the absence of CO: in the surrounding atmosphere, the dissociation of calcite starts at 500°C.However. when 1 atm of CO surrounds the sample, the dissociation does not start until 900C. The other carbonates are apparently much less sensitive to a change in partial pressure of CO2. The aragonite calcite transformation can occur an!. where within the temperature range of 357 to 488°C. depending upon the presence of impurities such as barium. strontium. lead and perhaps zinc. The differential thermal analysis (DTA) curves of both magnesitc and dolomite vary with the presence of 'impuritics. The presence of iron in particular seems to affect the magnesitc DTA curve

Jan 1, 1966

By N. Mungan

Laboratory water floods were performed in oil-wet and waterwet alundum and Torpedo cores, displacing a refined oil with n-hexylamine or Triton X-100 solution. Also, some floods were performed in which a sucrose solution was displaced with n-butyl alcohol. The Purpose ot the tests was to see if the oil recotlery could be increased and to examine what role the inter-facial tension reduction and wettability change play in the recovery mechanism. it was found that the chromatographic transport of amine was influenced by core wettability. In oil-wet cores, the rate of advance of the amine band could be predicted from equilibrium chromato-graphy, while in water-wet cores the amine band moved faster than predicted, indicating a non-equilibrium process. By reducing the interfacial tension to 1.1 dyne/ cm, oil recovery was increased. More oil was recovered from Dri-filmed cores than from water-wet cores. Reversing the wettability of the porous media from oil-wet condition also resulted in some additional oil recovery. Neutral wettability floods did not increase oil recovery. INTRODUCTION Interfacial forces in petroleum reservoirs are responsible for retention of large quantities of residual oil. Increasing exploration costs have created an incentive to attempt recovery of this residual oil by new and novel processes. One novel process involves changing the interfacial forces by introducing a chemical into the reservoir during water flooding. Recently, several investigators have carried out laboratory displacement tests using amines. lt3.4 The amines adsorbed onto initially water-wet core surfaces, changing them to oil wet. Additional oil was recovered when the cores were flooded in a manner reversing the. wettability from oil-wet conditions. In all these studies, the wettability reversal was considered responsible for increased oil recoveries although it was noted that amines also reduced oil-water interfacial tension. In the present study, displacement experiments were performed in water-wet and Dri-filmed cores. The purpose was to separate the effects of lowering the interfacial tension from changing the wettability. LABORATORY STUDY All experimental data given in this paper have been obtained at a temperature of 70F + 1 and prevailing atmospheric pressure. FLUIDS AND MATERIALS The oil used in these experiments* was washed with H2S04 and passed through two 100-200 mesh silica-gel beds to remove unsaturates and surface active impurities. Normal hexylamine** of practical purity was used, and the alundum cores*** were composed of nearly pure A1203. Fused quartz (SiOi) and Lucalox (A12 )t plates were used in the contact angle measurements. DISPLACEMENT TESTS Displacement tests were performed in alundum and Torpedo cores. Before displacement runs, cores were acidized with 0.1 normal HCI, flushed with distilled water and fired at 1,200F to make them water wet. This procedure also reduced water sensitivity of the Torpedo cores. Oil-wet cores were prepared from water-wet cores by flowing several pore volumes of 1 per cent by weight Dri-film SC-87 solution in heptane through the cores, flushing the cores with nitrogen, evacuating, and finally heating the cores at 300F for 4 or 5 hours. Between floods, cores were cleaned by conventional procedures and treated again as above.

Jan 1, 1965

By W. E. Howlett, R. L. Slobod

Jan 1, 1965

By F. W. Jessen, R. W. Durie

The dissolution of salt in the development of salt cavities is controlled by free-convection boundary layer flow along the salt surfaces. It is the purpose of this paper to expand upon results published previously, and primarily to illustrate the influence of surface features on the boundary layer and hence the dissolution process in the development of underground salt cavities. General dissolution equations are presented that apply to a surface inclined to the vertical. A salt surface inclined so as to be overlain by water is dissolved by boundary layer flow, while an inverted surface is dissolved by a combination of boundary layer and cellular flow. Detailed investigation of the boundary layer was limited to laminar free-convection flow. Truly laminar free-convection boundary layer flow may be encountered only in laboratory solution cavities, but the boundary layer will control dissolution even where cavity dimensions are such that free-convection turbulence will occur. The factors found to significantly affect the laminar boundary layer will affect dissolution with a turbulent boundary layer. INTRODUCTION Underground salt cavities have been used extensively for the storage of LPG products, and the use of these cavities is expanding each year. The increasing requirements for such cavities suggest the need for a better understanding of the mechanisms involved, for the purpose of the control of cavity geometry and the most effective use of the injected water. A previous paper outlined a quantitative approach to the determination of salt dissolution rates in the laminar free-convection dissolution of salt. It is the purpose of this paper to expand upon the results previously published and to point out significant factors in the application of these results to the dissolution of field cavities. The dissolution of salt in laboratory scale cavities is controlled by the free-convection boundary layer, which is the region of flow along the salt surface caused by the difference in buoyancy forces resulting from the salt dissolution. The transfer of salt into solution at the salt face causes an increase in salinity, and hence a greater density, in the water at the salt surface. The region of higher density tends to fall through the lighter, less saline water and a flowing boundary layer is the result. The water remains stationary at the contact with the salt surface,l so the resulting velocity profile is such that the influence of gravity is opposed by a viscous shear force at the salt surface. For salt to be transferred into solution from a solid surface to a boundary layer of fluid, it must diffuse through a stationary film of water at the salt contact. According to Fick's first law, the rate of mass transfer by diffusion is the product of the diffusion coefficient and the concentration gradient, with diffusion transfer in the direction of decreasing salinity. Once the salt has passed into the boundary layer, it will be transferred by a combination of molecular diffusion and convection. The convection transfer will in turn affect the rate of dissolution by influencing the gradient of the salt concentration at the surface. In the previous investigation by the authors to determine the basic mechanism of the salt dissolution process, an equation was developed that agreed very well with the experimental rates of free-convection salt dissolution from a smooth vertical surface. The technique was taken from the work of Eckert2 in the analogous heat transfer system and will be referred to as the integral method. The rate of salt transfer from a vertical section of height H was found to be:

Jan 1, 1965

By M. M. Mebta, G. W. Dean, W. H. Somerton

With the advent of underground heating operations, interest has developed in the alteration of rock properties by high-temperature treatment. In the present work a number of sandstones were heated to temperatures in the range of 400°C to 800°C under both atmospheric and simulated reservoir pressures. Pertneabilities increased by at least 50 per cent and sonic velocities and breaking strerlgths decreased by an equivalent amount. Differential thermal expansion and other reactions of constituent min-era1 grains are the causes of these alterations. INTRODUCTION In the underground combustion of petroleum reservoirs, temperatures of the order of 600C are reported to have been reached in the combustion zone.' At this temperature rocks are subject to extensive thermal alteration. Temperatures of this magnitude and higher may also occur in subsurface formations when subjected to bottom-hole heating, thermal drilling operations, and underground nuclear explosions. Temperatures of this magnitude might also be generated by conventional rock drilling methods at points of bit-tooth contact. In earlier work, the permanent deformation of rocks resulting from heating was reported. Major structural damage of rocks occurs due to differential thermal expansion of mineral constituents. A number of mineral alterations including crystal inversions, loss of water of crystallization and dissociation, may also contribute to changes in physical structure and properties of rock. In the present work, samples of three typical sandstones were heated to several temperatures up to a maximum of 800C and then allowed to cool to room temperature. Heating was done under both atmospheric pressure and simulated reservoir pressure conditions. Physical properties of the samples were measured before and after heating and comparisons made. Measured properties included permeability, sonic velocity, breaking strength and fracture index. Changes in physical properties were compared to changes in mineralogical characteristics as determined by thin-section. X-ray diffraction and chemical tests. EQUIPMENT AND PROCEDURE Two outcrop sandstones (Bandera and Berea) and one sub-surface sandstone (St. Peter) were selected for the tests. These samples have a wide ranee in composition and physical properties as shown in Table 1. The first series of tests was made on 2-in. diameter by 5-in. long test specimens. Test specimens used in all later work were 3/4-in. diameter by 1 1/8-in. long, this being the specimen size required for heating at simulated reservoir pressures. After careful washing, the cores were oven dried at 100 ± 5C for a minimum of 24 hours before the tests were run. Test specimens were heated in an electric furnace at a constant rate of temperature increase of 3C per minute. When maximum temperature of the run was reached, the sample was allowed to soak for one hour. The furnace was then cooled to room temperature at the same rate of 3C per minute. The entire heating operation was designed for reproducibility without subjecting the test specimens to excessive thermal shock. For samples heated under simulated reservoir pressures, a pressure cell designed by Dean was used (Fig. 1).3 The core sample was inserted into a thin-walled (0.006 to 0.01-in.) copper cup which was then mounted in a high-pressure cell. Provisions were made for the application of internal pore pressure as well as confining pressure. Tests showed that the thin-walled copper cup closed tightly around the core and satisfactorily transmitted confining pressure to the core. The core was heated by placing the entire cell into the electric furnace. The heating program was the same as that used in the atmospheric pressure runs: tempera-ture rise of 3C per minute to maximum temperature of the test, soaking at maximum temperature for one hour, and cooling at a rate of approximately 3C per minute. The cell was designed to withstand 5,000 psi at 1,000C. However, since it was considered likely that repeated heating and cooling would in time weaken the steel, 2,000 psi at 850C was set as a working limit. In the present series of tests, the pore pressure was held constant at 750 psi and the confining pressure at 1.500 psi. The pressure source was a high-pressure nitrogen tank. The two pressures were controlled manually but are accurate well within ± 50 psi.

Jan 1, 1966

By J. D. Rice, P. L. Wearden, W. L. Hensou

Solutions to the unsteady-state partial water-drive reservoir performance problem can be obtained through the use of analogue computers or high-speed electronic digital computers. The solutions that have previously been resolved for use on digital computers, however, demand a knowledge of the aquifer parameters. Generally, the analogue computers now in use do not require this knowledge. A solution is presented herein where the aquifer performance, expressed in terms of difference equations, is related to the reservoir performance as expressed by the modified Scbiltbuis material-balance equation. A numerical procedure for a medium-sized digital computer also is presented in which a solution to the set of equations defining the aquifer and reservoir performance is obtained and the aquifer parameters (permeability and sand thickness) are automatically optimized while simultaneously matching the known pressure behavior. Predicted pressure behavior can be calculated using rates from any assumed future production practice. The procedure provides an output format which presents the cumulative, incremental and average rate of natural water influx, the per cent gas-cap expansion, the calculated reservoir pressure, the measured reservoir pressure, and the difference between the measured and calculated pressures. Results of a test problem are presented in comparison with results obtained by the Bruce Analyzer, Ohio Oil Co. 's Pace General Purpose Analogue Computer, and Sun Oil Co.'s Single-Pool Electronic Reservoir Analyzer. These results indicate that unsteady-state reservoir perlormance for a single-pool system can be adequately simulated by a numerical method employing a digital computer and that the special-purpose analogue computer can be supplanted by this method. INTRODUCTION The numerical approach to the unsteady-state partial water-drive reservoir performance problem heretofore has been limited in that solutions have been resolved only for certain fixed parameters. Also, prior to the advent of high-speed electronic digital computers in the last few years, calculations were so tedious and time-consuming that they ordinarily were neglected, or some other means of obtaining the solutions were used. The analogue approach had been used for obtaining solutions to the transmission-of-heat problem and, consequently, means were developed for applying the electrical analogue to the unsteady-state reservoir performance problem. Bruce submitted a paper for AIME publication in Sept. 1942, in which he described an electronic device (now called the Bruce Analyzer) for analyzing the water-drive performance of any reservoir, based on either steady or unsteady-state flow characteristics. 1 Sun Oil Co. adopted the analogue approach in 1949 and now has two special-purpose high-speed electronic reservoir analyzers in use.2,3 Two main features of Sun's analyzers are (1) the high repetitive rate of reproducing the entire history of reservoir performance cyclically a number of times per second, and (2) the elimination of the necessity of knowing the properties of the associated aquifer. Since the advent of larger and faster digital computers, the possibility of incorporating this second feature into a numerical approach appeared feasible. Previous work has been done and material published covering the application of numerical methods to the behavior of the pressure with respect to time at any point in an associated aquifer? However, the aquifer parameters must be known from some source. The material presented in Ref. 4 has been used to a great degree in developing the numerical method described herein. DIFFERENTIAL EQUATION AND BOUNDARY CONDITIONS The problem of simulating the performance of a reservoir and the associated aquifer is to determine the characteristics of the aquifer such that the