Search Documents

Jan 1, 1966

Jan 1, 1966

By C. R. Johnson, R. A. Greenkorn, R. E. Haring

Three confined five-spot miscible displacements at unity, favorable, and unfavorable mobility ratios were conducted in a shallow, water-saturated sandstone of Pennsylvan-ian age near Chandler, Okla. These studies, plus associated laboratory experiments, were designed to measure miscible displacement performance in a controlled natural system, using known scaling criteria to develop an approach to modeling the heterogeneous field system. We have concluded from these studies that: (I) displacement efficiency in the field is a pronounced function of mobility ratio, indicating that miscible fingering observed in simple laboratory models occurs in the field: (2) field displacements can be quantitatively predicted by scaled laboratory models if the degree and location of field permeability variations are preserved in the models: and (3) arbitrary simplifications of heterogeneity will not necessarily predict observed displacement efficiency, and the simpler the model, the more optimistic the prediction. INTRODUCTION Many processes to achieve miscible displacement of reservoir oil by injected fluids have been conceived and field tested by the oil industry. Among the better known are high-pressure gas, enriched gas, and LPG banks. The simplest form of miscible displacement—one fluid miscibly displacing another fluid of different viscosity but the same density—has been studied extensively in homogeneous laboratory models. Observations of unstable fingering have been made which explain the significant decrease in displacement efficiency as the mobility ratio (ratio of displacing fluid mobility to displaced fluid mobility) increases. General industry experience with field tests of miscible displacement projects, mostly LPG banks, has been premature solvent breakthrough and lower than predicted production rate increases. These results have been attributed to either unstable fingering, unusual or unexpected permeability stratification, or both. Miscible displacement data in a controlled natural system have never been reported, however. Also, it has not been shown that properly designed and constructed laboratory models quantitatively predict field-scale behavior. The purpose of the combined field and laboratory ex- periments reported in this paper is twofold. The first was to measure miscible displacement performance at different mobility ratios in natural rock approaching field size under precise, controlled conditions. The second purpose was to utilize known scaling criteria plus several approaches to heterogeneity to model the field. Comparison of model and actual field results should then determine whether or not the laboratory phenomena (manifested by miscible displacement efficiency) are exhibited in large, natural rock systems. We carried out our program by first locating a shallow, water-saturated reservoir whose rock properties were representative of oil-bearing reservoirs. Detailed reservoir description by core analysis and interference testing showed the field site to be heterogeneous. A sequence of controlled, aqueous-phase miscible displacements was conducted at unity, favorable and unfavorable mobility ratios. A central, confined pattern was used to obtain the displacement data. A laboratory program using sand-packed models was conducted to determine the modeling criteria necessary to simulate field behavior of miscible displacement in a heterogeneous system. SCALING THEORY The detailed derivations and descriptions of the scaling laws that apply to laboratory models of reservoirs are adequately described elsewhere, so the following discussion will be restricted to facets of importance in this study. For a displacement in which one liquid miscibly displaces another, the following dimensionless groups are required to have the same numerical value in the model as in the field: The model also must be geometrically similar to the field, be spatially oriented the same as the field (same dip angle), and have the same initial and boundary conditions as the field (same initial fluid saturations and same injection-production well arrangement). When these conditions are satisfied, the theory predicts that at dny dimensionless time (pore volumes of produced fluids) the dimen-sionless flow potential and dimensionless fluid concentrations will be identical at all dimensionless spatial locations within the model and the field. If this prediction is correct, then the local dimensionless velocities must be identical, thus the instantaneous fraction of displaced fluid produced and the cumulative recovery expressed as fraction of original fluids in place must be identical at all dimensionless times. The theory outlined above has been obtained by either dimensional analysis or inspectional analysis of differential

Jan 1, 1966

By J. Handin

Triaxial compression tests with independently applied external confining pressures and internal pore pressures show that the ultimate compressive strengths of representative oil well cements are nearly linear functions of effective pressure — the difference between external and internal pressures on the jacketed cylindrical specimens (to 15,000 psi). The strengths are little affected by the test temperature to 350F (not to be confused with the curing temperature). At an effective pressure of 15,000 psi, strengths are in the rage of 30,000 to 50,000 psi, comparable to those of sedimentary rocks under similar conditions. The cements become very ductile even under low effective pressures; permanent shortenings of 30 per cent or more are attainable without rupture. INTRODUCTION Since the pioneering work of Richart, Brandtzaeg and Brown1 on the failure of cement under combined compressive stresses, it has been recognized that ultimate compressive strength is greatly enhanced by the application of confining pressure. More recently, Mc Henry2 showed that the strength of concrete was a linear function of the effective pressure (the difference between the external confining pressure on a jacketed specimen and the internal fluid pore pressure), at least for a range of 0 to 1,500 psi. The effect of temperature had not been investigated, and no previous systematic triaxial compression testing of naterials used for oilwell cementing seems to have been done. The present work was suggested by the late J. hi. Bugbee3 who stated that "consideration of the common application of high-pressure hydraulic fracturing to the initial completion or recompletion of wells, and the large pressure drawdowns in some producing wells, particularly those in abnormally high-pressure gas-condensate reservoirs, raises the question of what is a suitable cement strength for various completions. The intuitive belief exists that cement strength need be no greater than formation strength. Tests should, however, be conducted at downhole conditions." The ultimate compressive strengths of rocks penetrated by the borehole must rise several fold with increasing depth.4 This marked enhancement of strength is due to the influence of the effective pressure, the total weight per unit area of the overburden less the hydrostatic pore pressure. (The effect of temperature due to the geothermal gradients is relatively small for depths to 30,000 ft.) .A significant comparison of the strengths of rocks and cements at downhole conditions requires knowledge of the confined compressive strengths of cements as well. EXPERIMENTAL PROCEDURES The theory and technique of triaxial compression testing are fully discussed in earlier reports.4-6 Briefly, cylindrical specimens 1-in. long and 0.5-in. in diameter are jacketed in thin copper tubes of negligible strength, placed in the test chamber and subjected to an external confining pressure of kerosene and loaded axially by the piston at a strain rate of 1 per cent per minute. Pore pressures of water (or kerosene) are applied independently through the hollow piston and are maintained constant during the shortening of the specimen, Tests at sensibly 0 pore pressure are arranged so that any free water in the cement can escape to the atmosphere. (If egress of water were denied, pore pressure would rapidly attain the value of the external confining pressure because of reduction of pore space.) The test chamber can be heated for high-temperature experiments. Unless otherwise noted, the cement samples were air dried for about a week. Recorded during a test are pore and confining pressures, shortening and axial differential force (total force less the product of the confining pressure and the area of the piston). True stress-strain curves are constructed by dividing the measured differential force by the actual cross-sectional area of the specimen at a given

Jan 1, 1966

By K. P. Fournier

This report describes work on the problem of predicting oil recovery from a reservoir into which water is injected at a temperature higher than the reservoir temperature, taking into account effects of vis cosity-ratio reduction, heat loss and thermal expansion. It includes the derivation of the equations involved, the finite difference equations used to solve the partial differential equation which models the system, and the results obtained using the IBM 1620 and 7090-1401 computers. Figures and tables show present results of this study of recovery as a function of reservoir thickness and injection rate. For a possible reservoir hot water flood in which 1,000 BWPD at 250F are injected, an additional 5 per cent recovery of oil in place in a swept 1,000-ft-radius reservoir is predicted after injection of one pore volume of water. INTRODUCTION The problem of predicting oil recovery from the injection of hot water has been discussed by several researchers. l-6,19 In no case has the problem of predicting heat losses been rigorously incorporated into the recovery and displacement calculation problem. Willman et al. describe an approximate method of such treatment. 1 The calculation of heat losses in a reservoir and the corresponding temperature distribution while injecting a hot fluid has been attempted by several authors.7'8 In this report a method is presented to numerically predict the oil displacement by hot water in a radial system, taking into account the heat losses to adjacent strata, changes in viscosity ratio with temperature and the thermal-expansion effect for both oil and water. DERIVATION OF BASIC EQUATIONS We start with the familiar Buckley-Leverettg equation for a radial system:* This is sometimes referred to as the Lagrangian form of the displacement equation. If the injection water is at some temperature Ti which is above reservoir temperature TR, then lw, the flowing fraction which is water, is a function of both water saturation and fluid temperature at any given time and position in the reservoir. The temperature dependence of lw results from the temperature dependence of the viscosity ratio µo/µm of which lw is a function. Thus:

Jan 1, 1966

By J. van Lookeren

In the case of a reservoir where the oil underlies a large gas cap and overlies bottom water, production can be inzproved considerably if wells are perforated below the water-oil contact rather than opposite the oil-bearing section. The theoretical arguments for this improvement, which in some cases may be tenfold, are supported by experimental evidence. INTRODUCTION Reservoirs frequently consist of one sand body containing an oil layer underlying a thick, gas-bearing section and overlying bottom water (Fig. 1). Wells in such reservoirs are customarily opened up opposite the oil-bearing interval or, rather, opposite the lower part of this interval (Fig. 2). More often than not, it is inefficient or prohibited to blow off large quantities of gas with the oil. This im-psses a severe limitation on the drawdown, which results in low oil-production rates. Because of the small difference in specific gravity between oil and water, the latter will easily cone upwards, so that the production of water together with the oil is not prevented by this completion method. In such a situation wells should be opened deeper, in the originally water-bearing zone. The gas can be shut off by not perforating over the upper part of the interval (Fig. 3) or by a tubing hung close to the bottom of the interval (Fig. 4). Initially much water will be produced, but when steady-state conditions have been established (as depicted in Figs. 3 and 4) larger drawdowns are permissible than in a conventionally completed well, and correspondingly higher oil-production rates will result. Water cuts are only slightly higher because of the small effect of gravity on the water-oil flow distribution. CALCULATIONS The ratio of improvement has been estimated by means of approximate calculations. The completion depicted in Fig. 4 is the most convenient basis for these calculations. By postulating that vertical flow gradients can be disregarded, the Appendix shows that the well-known Dupuit formula can be obtained except for some complication due to the presence of two fluids rather than one. The Appendix makes the reasonable assumption that, for suitable boundary conditions, the oil layer occupies a constant fraction of the total fluid height throughout the reservoir. The oil- and water-production rates can thus be related to the liquid levels in the borehole as follows:

Jan 1, 1966

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 T. K. Perkins, R. N. Hoffman, O. C. Johnston

A simplified theory of viscous fingering in miscible systems has been developed. It predicts the correct functional relationship between pertinent variables and permits the calculation of the order of magnitude of fingering behavior for simple systems, such as linear or radial. The theory is based on the following four observations: (1) cross-flow takes place only near the ends of fingers; (2) the relative finger width is about 0.5; (3) fingers can be suppressed by transverse dispersion, the suppression being quantitatively described by a critical ratio of dispersion times; and (4) the widths of extending fingers are close to the minimum size finger that can grow at any stage of displacement. Fingering is studied in two-dimensional linear and diverging radial systems, both theoretically and experimentally. For linear systems, the length of the fingered region is proportional to mean displacement; the finger width is proportional to square root of mean displacement; and there is a small initial region void of fingers because of suppression by longitudinal dispersion. For the radial system, two limiting conditions are recognized. If the mean displacement is small compared with the wellbore radius, the length of the fingered region is proportional to the mean displacement. The width is proportional to the square root of mean displacement. If the mean displacement is large compared with wellbore radius, length of the fingered region is proportional to mean displacement, but the number of fingers approaches a constant value. Also, in the radial case there is a small initial region void of fingers because of longitudinal dispersion. INTRODUCTION The behavior of viscous fingers in miscible systems has been of interest to the oil industry for many years. Previous studies have clearly shown the existence and nature of fingers in small models.1,2,3,4 Engineering techniques for extrapolating to reservoir situations have been proposed. 5,6 Still, because of the lack of a real understanding of the mechanics of fingering, there remains uncertainty and disagreement as to the best way to scale or calculate fingering behavior in the reservoir. In this paper we discuss a study of the mechanics of fingering in miscible systems (i.e., why do fingers behave as they do?). A simplified theory is developed which we believe will (1) predict the correct functional relationship between pertinent variables, (2) permit us to calculate the order of magnitude of fingering behavior for simple systems such as linear and radial, and (3) give further insight into the problem of scaling or otherwise extending the results to more complicated reservoir conditions. The paper includes the following sections: (1) brief summary of the mixing behavior of miscible fluids in linear and radial systems; (2) four fundamental observations of fluid flow (under fingering conditions); (3) based on these four observations, a simplified theory of fingering in linear and radial systems is developed; and (4) the theoretical equations are compared with experimental fingering data measured in laboratory models. A REVIEW OF DIFFUSION AND DISPERSION As will be shown later, the behavior of viscous fingers in miscible systems is controlled in large degree by the mixing between the two fluids. A quantitative description of fingers will first require a quantitative description of mixing. Fortunately, much work has been done to clarify this subject; a fairly comprehensive review has been presented in a previous paper. 7 For present purposes let us (1) briefly summarize the quantitative description of diffusion and dispersion (both longitudinal and transverse) within a differential element; (2) present a simplification of the "width of mixed zone"; and (3) describe the effect of geometry on width of mixed zone. If one fluid is displaced from a porous medium by another miscible fluid (in the absence of fingers)

Jan 1, 1966

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

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

Jan 1, 1966

By C. R. Kyle, R. L. Perrine

A transparent model of a reservoir has been used to study some characteristics of instability in miscible displacement. The linear model dimensions are 1/4 in. x 91/2 in. x 10 ft. The model is packed with spherical glass beads. Thee groups of experiments have been performed. The objectives were to measure, respectively, the rate of growth of the transition zone, flow velocity variations in the transition zone, and the effect of a graded viscosity zone on stability. All experiments were compared, where possible, with theory. The transition zone grew almost linearly with distance traveled, regardless of flow rate, where adverse viscosity ratios were used. At very adverse viscosity ratios and large distances traveled, some reduction below linear growth was noted. This could be attributed to the effects of mixing. General similarity in flow pattern was commonly shown in different tests, although exceptions occurred at the higher viscosity ratios. Injection of tracer spots permitted direct experimental measurement of velocity variations. Statistical methods were used to compare observed velocity variations with recently published theories. Surprisingly good agreement with some theoretical predictions was obtained. In the experiments using an initially graded viscosity zone, macroscopic stability was reached before predicted by theory. Under conditions theoretically considered to be slightly unstable, any fingering tendency was too slight to become observable within a 10-ft long system. The general correlation of all experiments with published theory was good, and suggests that theoretical predictions can be extended to reservoir problems. INTRODUCTION Early laboratory studies of miscible displacement demonstrated the high recovery efficiency possible with this method. Linear laboratory models showed that, under some circumstances, all in-place hydrocarbons could be recovered using as little solvent as 5 per cent of the pore volume. These results led to a series of pilot field projects. Under realistic conditions, however, other laboratory studies showed that a small solvent bank may not recover much more oil than a conventional gas drive. The cause of this lack of efficiency is unstable viscous fingering. The size of the solvent bank required to give optimum recovery is, therefore, a question of great economic significance. The object of the experiments described in this paper was the quantitative measurement of unstable flow parameters, the verification of present theories on the miscible displacement process, and correlation of these theories with ex peri mental observation. Of the theories which have been advanced, those by Dougherty, Koval and Perrine show promise of predicting reservoir performance with reasonable accuracy.l-7 The present experiments were correlated, where possible, with the theories of Koval and Perrine. Dougherty's theory was not used because of the additional computational difficulty it presented and the belief that no further information would result. EXPERIMENTAL EQUIPMENT A transparent reservoir model, consisting of a vertical reinforced plexiglass column, was used for the experiments. The column was packed with glass beads, which served as the porous medium. To induce a known velocity disturbance in the flow, a double cosine wave restriction in the transverse cross-section was placed at the top of the column. To mark the flow field and provide a means of measuring the local fluid velocity, a dye-spot injection system was provided. The column was 10 ft high by 1/2 in. wide by 1/4 in. deep. The cosine restriction was placed 16 in. from the top of the column. The glass beads were smooth and spherical and had a mean diameter of 0.0102 cm. Liquids entered the top of the column through a manifold allowing three points of entry (left, center, and right, or any combination of these). A positive-displacement pump with variable-speed drive was used to control liquid flow rates. The pump transferred viscous oil from a reservoir to an accumulator, which displaced the test liquid into

Jan 1, 1966

Jan 1, 1966

By J. A. Jargon, H. K. van Poollen

Flow of slightly compressible fluids through porous media can be described by linear equations. Many mathematical solutions- are available for the boundary conditions of constant terminal pressures or constant flow rates, or a combination of both. Various techniques of formation evaluation use drastic rate changes, such as sudden shut-in. Other methods require that the well be produced at a constant rate. In practice, one is not always able m fulfill these rigid requtirements. This report shows a methor1 by which varying-rate data can be converted to a constant-rate pressure response, called the unit response function. After obtaining the unit response function, one may analyze it, using techniques described in the literature. The method uses tile principle of superposition. Necessary data is a time history of pressure and production. If pressure changes and/or rate changes are very drastic, the unit response function becomes erratic and oscillates The erratic and oscillating results are obvious, and under those circumstances, no answer can he obtained. It appears, however, that a response function without oscilla-tion will not contain any appreciable errors —- one gets either good results or no results. This report includes a mumber of examples of converting constant-rate data to constant terminal-pressure data and vice versa. In addition, a number of field-test applications are given. This method should be useful in converting mathematical solutions from one boundary condition to another. It already has proven useful in converting varying-rate data to constant-rate conditions. INTRODUCTION Pressure drawdown tests are performed on production wells to determine various reservoir characteristics, including permeability, distances to discontinuities, reservoir continuity and wellbore conditions. An ideal drawdown test is performed at a constant production rate with pressure observation during the drawdown period. A typical drawdown curve for a constant rate of production is shown in Fig. 1. In practice, a constant production rate cannot always be maintained. Fig. 2 shows a pressure drawdown curve for a varying rate. Usual methods of analysis are not applicable to tests such as this. It is, therefore, necessary to convert the varying-rate data somehow, before usual analytical techniques can be applied. This report describes a method of obtaining a response function. A response function can be defined as the relationship between pressure and time for a well producing at a constant unit rate. Such a response function can be analyzed by the usual analytical techniques. Hutchinson and Sikora' first approximated a resistance function (response function), and by applying the superposition principle, pressures were calculated using the varying rates. If no match was obtained, adjustments were made to the response function and the closer approximation was used again in an attempt to match calculated and observed pressure behavior. If a good match was obtained: the response function was then used to extrapolate the future aquifer behavior.

Jan 1, 1966

Jan 1, 1966

By Oren C. Baptist, Carlon S. Land

Laboratory research has been conducted to evaluute the effect of clay hydration on the permeability to gas of water-sensitive reservoir sands. Samples of a .sandstone containing trace amounts of montmorillonite and a sample of montmorillonite were .studied in the laboratory to detertnine whether swelling or dispersion was the cause of permeability reduction in these samples. Heliuin, containing various amounts of water vapor, was used to hydrate the clay minerals and to determine the gas permeability at various stages of clay hydration. The amount of water adsorbed by the samples using this method is small. The nonwetting-phase permeability at higher water saturations war investigated by saturating the with water and measuring the permeability to humid helium while decreasing the water saturation, Relative-permeability curves obtained from results of these procedures were used to estimate the effect of the swelling of trace amounts of mont/tlorillonite on the permeability of the .samples. Most of the damage to the permeability when reservoir sands containing trace amounts of montmorillonite are exposed to fresh water is due to dispersion and movement of clays. Blockage of pores by the increased volume of expanded montmorillonite is believed to result in permeability damage that is small in comparison to the observed damage to the samples tested. INTRODUCTION Studies have shown that permeability is severely damaged when sands containing only small amounts of montmorillonite are contacted by fresh water.15 When samples of sands containing large amounts of montmorillonite are placed in fresh water in the laboratory, these samples may completely disintegrate, forming an unconsolidated mass of larger volume than that occupied by the dry sample." In this case, it is apparent that the swelling of montmoril-lonite has destroyed the pore structure of the sand. If only a trace of montmorillonite is present in a sand. samples may remain intact when saturated with water, although the permeability to water is a small fraction of the gas permeability of the dry sample. Many workers in the field of water sensitivity have attributed this reduction in permeability to the blocking of pores and reduction of pore size by the increased volume occupied by expanded mont- niorillonite. if the sand contains a detectable amount of montmorill'onite or mixed-layer clay containing rnontmorillonite. Logically3 the smaller amount of montmorillonite present in a sand, the smaller should he the effect of montnlorillonite swelling on permeability; however, the quantity of montmorillonite sufficient to cause severe damage by swelling is not known. Although hundreds of samples have been tested in our laboratory, no correlation has been established between the amount of montmorillonite in samples and the permeability reduction caused by fresh water. To many petroleum engineers, the phrase "clay swelling" is synonymous with "water sensitivity", or "permeability reduction" implying that any formation damage due to the hydration of clays is caused by swelling. Although all clays adsorb water on their surfaces, montmorillonite is the only clay mineral commonly found in reservoir rocks which adsorbs water between intercrystalline layers, resulting in expansion of the clay particle. As montmoril-lonite swells, the first few layers of water adsorbed between platelets are strongly held and well oriented, and the montmorillonite retains its crystalline structure, although expanded. As swelling of sodium montmorillonite continues, the platelets become farther apart and the forces orienting the platelets in the crystalline structure become weaker, resulting in a less orderly orientation of platelets. In an abundance of water, small groups of platelets may become detached from the original monl-rnorillonite particle and may be dispersed throughout the water phase. Because of its swelling properties, sodium montmorillonite is very easily dispersed in water. Particles of other clay minerals. such as illite and kaolinite may also be dispersed in water. causing water sensitivity of sands not containing montmorillonite. The presence of an immobile layer of water adsorbed on the surface of clays has been considered a possible cause of the low permeability to water of dirty sands. Grim states that the thickness of the layer of immobile water held by sodium montrnorillonite is three nlolecular layers or 7.5 A (angstroms), with some orientation of water extending to 100 A. Assuming a very thick, immobile water layer adsorbed on the surface of a pore represented by a capillary tube, the maximum effect of the water layer on permeability can be calculated. Using a pore radius of 10 ' cm and an immobile water layer of 50 A. the calculation shows the permeability to be reduced only 2 per cent. Similar calculations can be used to show that the effect of electro-osmotic counterflow is of the same order of magnitude as that of bound water. The reduction of the permeability to water by either an immobile water layer

Jan 1, 1966

By L. H. Robinson

A proposed method of drilling utilizes sequential detonation of two types of explosive charges delivered to the hole by a conventional drilling fluid through pipe. A shaped charge first produces a long thin hole. A second charge, called a gauging charge, is pumped into the thin hole, explodes, and enlarges it to full gauge. In one suggested configuration, the explosives are contained in cylindrical capsules 4- to 5-ft long. A drilling head is positioned 2- to 3-ft off bottom. When a capsule is pumped to the bottom of the hole and stops, pressure inside the drill pipe increases and differential pressure detonates the explosive. Primacord along the length of the capsule disintegrates the capsule. The debris circulates out of the hole with the drilling fluid. The gauging-charge action was examined using 2-ft cubes of Berea sandstone. The sandstone (18 per cent porosity) was vacuum-saturated with water and surrounded by concrete. To prevent reflected tensile pulses, the acoustical impedance of the concrete matched that of the sandstone. A 750-psig pressure was applied to a drilling fluid within a shaped-charge-formed hole before detonating the gauging charge. In these tests, 190 gm of rock were removed for each gram of explosive detonated. A similar series of experiments on low-porosity (0.40 per cent) metamorphic limestone provided a more rigorous test. Extensive fracturing from the gauging-charge action prevented an estimate of the volume of the expanded hole. The diameter of the hole was irregular, ranging from about 2 to 4 in, or about fivefold larger than the hole created by the shaped charge. INTRODUCTION In the past, explosives have found a varied but somewhat limited use in the oil field. An early well stimulation technique, still in use today, in vo lves shooting ni troglycerin to en large the borehole and to increase matrix permeability near the wellbore. Bullets and jets are now the conventional method of perforating casing and tubing. Detonations are frequently used to "back off" or "shoot off" stuck pipe and, of course, seismo-graphic work has long depended on explosive-generated pulses to find attractive earth structures. Efforts to use explosives in drilling operations have received limited examinations1-3 which, in one respect, is surprising because explosives offer one way to transmit large amounts of energy directly to the bottom of the hole. Recently, Ostrovskii l described a Russian method of drilling with explosives that utilizes small capsules containing explosives which are pumped into the mud stream to the bottom of a borehole. As the capsules leave the end of the drill pipe, a detonation mechanism is armed. Impact detonates the explosive. In another configuration, two liquids (an oxidizing agent and fuel) in the capsule are separated by a frangible membrane. Separately, the liquids are nonexplosive. Near the bottom of the drill pipe, the capsule passes through a narrow constriction which ruptures the membrane. The mixing of the liquids produces an explosive. Impact is again used to detonate the explosive. About the time Ostrovskii was developing his method, Humble Oil & Refining Co. was independently doing research on another explosive drilling concept based on the sequential detonation of two types -of explosive charges. The proposed technique consisted of alternately pumping a shaped charge and then a "gauging charge" down a well. The shaped charge should drill a long, tapered hole, and the gauging charge should expand it to full gauge. This paper presents the visualized drilling concept and gives results of laboratory work which showed that sequential firing of charges can form hole at an attractive charge-to-rock removal ratio. An important finding was that increased environmental pressure, which is to be expected in wells, increases rock removal efficiency of the gauging charge.

Jan 1, 1966

By R. F. Burdyn

The inadequate use of centrifugation to economically recover solids from weighted drilling fluids reflects the need for better equipment and techniques for this putpose. Laboratory studies in the development of an improved separator are described in which an operating equation is derived and tested. Results show that the concentric cylinder geometry employed effectively separates bnrite from a suspension and that the operating equation provides a good approxinmtion for scale-up. INTRODUCTION Our current drilling technology frequently requires a high-density drilling fluid obtained by addition of barite. In the course of drilling, formation solids which are too fine to be removed either by screening or settling become suspended in the drilling fluid and gradually the volume of solids in the mud increases. The volume fraction of solids must be limited (if a satisfactory set of rheological parameters are to be maintained). A centrifugal separatorprovidesan economical way of accomplishing this. The barite recovery process can be considered as a separation of two solids. One, the light solids, composed of formation and added solids, has a specific gravity of 2.6 to 2.7; the other, barite, has a specific gravity of 4.2 to 4.3. This density difference, plus the fact that the average light-solids particle size is much smaller than the average barite particle size, permits separation by a centrifuge. In drilling fluids some of the coarse particles of the light-solids-range will settle faster than fine particles of the barite-particle range. As a result a complete separation of the two species is not possible. since the object of the process is not merely recovering the maximum amount of barite but includes as well removing the maximum amount of light solids, an optimum barite recovery efficiency exists. From a practical standpoint this optimum cannot be determined in the field for each drilling fluid system, and in practice the separation is less than optimum, with some sacrifice of barite. Drilling technology has included centrifugal separators for barite recovery for more than a decade. Results have been reported by a number of investigators indicating that the process is practical and economica1. l-4 The decision to use a centrifuge is based on economics in which direct cost savings and the indirect benefits in rig time derived from improvement of the drilling fluid are important factors. One would expect that centrifugal separation of barite from drilling fluids would significantly affect barite consumption; however, this is not the case. The Minerals Yearbook shows an annual domestic barite consumption in the drilling industry of nearly 1 million tons.5 By rough estimate there are perhaps 80 separators presently in field use. Assuming half of these in use at any one time, operating an average of four hours per day, at recovery rates averaging 3,000 Ib of barite per hour, total annual recovery is about 90,000 tons. This is less than 10 per cent of he total barite used. I conservatively estimate that barite consumption in drilling operations can be reduced by 30 per cent through greater utilization of centrifugal separators. To encourage more widespread acceptance of centrifugal separators in the drilling industry, improved equipment and techniques would be very desirable. The present paper, covering theory and results obtained from a laboratory model, is the first in a series on he development of an improved mud separator for field use. THE CONCENTRIC CYLINDER GEOMETRY AS A SEPARATING DEVICE Consider the geometry shown in Fig. 1, consisting of two concentric cylinders separated by an annular space. These are arranged so that the outer cylinder is fixed and the inner one can be rotated about its axis on shafts sealed against the ends of the outer cylinder. If the outer cylinder is fitted with openings at either end, the inner

Jan 1, 1966

By N. Mungan

Formation damage, i.e.. reduclion in permeability, has been generally attribuled to clay minerals which expand or disperse upon contact with water that is less saline than the connate water. Luboratory, studies show that penneahility reduction can also occur in formalions containing only nonexpandable clays such as illite or kaolinite, and can be caused also by changes in pH. Furthermore, pH changes can damage even formations that are essentially free of clays. It is suggested that permeability reduction is due to the small passages being blocked by particles, which may be dispersed clays, cemenlion material or other fine parricles. These particles are dislodged by dispersion of clays due to changes in salinity or by dissolution of calcareous cement by acids, or of silicaceous cement by alkaline solutions. In working with reservoir cores, it was found that extracted cores damaged more easily and extensively than nonextracted cores. The extent of damage depended also on tenlperatltre. INTRODUCTION Permeability is an important property of porous media and has been the subject of many studies by engineers and geologists. Many of these studies are conccrned with formation damage, i.e., reduction in permeability, resulting from exposure of oil-producing formations to water substantially less saline than the connate water. This effect causes understandable concern since during drilling, completion and production phases formations are often exposed to fresh water. The damage resulting from contact with relatively fresh water has been attributed to expansion and dispersion oF clay minerals. During laboratory investigation of the use of NaOH as a wettability reversal agent to increase oil recovery from oil-wet reservoirs, several cores used in the displacement studies suffered loss in permeability. Despite the traditional usage of NaOH for conditioning aqueous mud systems, the role of the caustic filtrate in wellbore damage seems to have been overlooked. Browning2 as recently reported on the effects of NaOH in dispersing clay minerals but he was concerned only with complications that may arise in drilling massive shale beds. The following study was made to examine the role of pH and salinity changes in core damage. Where cores from reservoirs were used, tests were performed with extracted and nonextracted cores both at room and reser- voir temperatures, since it was felt that the test environment and core condition may affect the results. Because of its limited coverage and exploratory nature, this study is not intended to provide answers to field formation damage problems. It is hoped that it will encourage research into new aspects of the permeability reduction problems, particularly those allied to new recovery and production processes. PROCEDURE In all permeability tests, fluids were pumped through the cores at a constant volumetric rate. Only deaerated fluids and reagent grade chemicals were used. The fluids were passed through two ultrafine filters before injection to remove any entrained particles. The cores, with the exception of the unconsolidated cores, were mounted in Hassler holders. Water was used to transmit pressure to the sleeve. The inlet endpiece had two entry ports which permitted scavenging one fluid with another to avoid any mixing in the small holdup volume. The cores were flushed with CO2 gas, evacuated for 5 to 6 hours and saturated with the first liquid at a pressure of 1,000 psi for 24 hours to eliminate any free gas from the cores. Pressure differences up to 20 psi were measured by transducers, calibrated in inches of water and continuously recorded. For greater pressure drops, gauges were used. All reservoir cores were cleaned with a light refined mineral oil, then with heptane, and finally dried with CO2. Compatibility tests showed that no precipitates formed when mineral oil and the crudes were mixed. Some cores were extracted in Dean-Stark-type solvent extractors using xylene and trichloroethane and dried in a vacuum oven at 450F. Each test consisted of a sequence of water, test solution, and again water flow. RESULTS AND DISCUSSION STUDIES IN BEREA CORES Salinity Contrast Berea cores 2-in. in diameter and 12-in. long were cut from sandstone quarried in Cleveland, Ohio. The clay minerals were identified by X-ray diffraction to be chlorite, kaolinite, illite and incerlayered illite. Flow of fresh water or 30,000 ppm brinc does not cause any permeability reduction (Fig. 1). However, after injection of brine the core is readily damaged by fresh water. Damage starts almost instantly as the fresh water injection is begun, and at a cumulative injection of 1.2 PV fresh water, the permeability has dropped from 190 to 0.9 md. Upon continued injection, the effluent contains clay minerals dislodged from the core. The final core per-

Jan 1, 1966

By A. Satter

Studies of wellbore heat transtnission during the injection of a hot fluid, as either gas or liquid, have appeared in he literature. The present investigation takes into account the effect of condensation, which is of practical sig-nificance when considering steam injection operations. Calculation procedures are given for superheated and for saturated and undersaturated steam. Effects of injection rate, time, pressure, temperature and depth of a well on heat losses have been analyzed. The benefits of using a packer are discussed. Also presented are correlations for estimating heat losses involved during injection of saturated steam. INTRODUCTION Studies of wellbore heat losses during the injection of a hot fluid have appeared in the literature. Ramey' has presented a useful solution for calculating temperature of a heated gas as a function of depth and injection time. Ramey and others have also presented solutions for the hot liquid case.'. ' No method is available to take into account the effect of condensation, which is of practical significance when considering steam injection operations. This paper presents a method of estimating the quality of a condensing fluid as a function of depth and time. The approach is basically the same as that used by Ramey. The following injected fluids and conditions are considered in this study: 1. Injection of superheated steam A. Decrease in temperature until cooled to the saturated state B. Condensation until condensed completely into hot (saturated) water 2. Injection of saturated or undersaturated steam; condensation until condensed completely into hot (saturated) water. Although this investigation does not include cooling of water after steam is completely condensed, Ramey's solution for the hot liquid case can be used for this purpose. The effects of various process variables on heat losses During the injection of steam into a well are investigated. THEORY QUALITY OF CONDENSING STEAM Predicting the behavior of condensing steam flowing down the injection wellbore requires an estimate of the quality of steam, i.e., the mass fraction of vapor in the mixture. Let us divide the total depth of the injection well into several segments which might be, but do not have to be, of equal magnitude. Consider condensation within a given depth interval AZ, the bottom of which is located at a depth Z from the surface (Fig 1). The following equation, derived in the Appendix- relates the quality of steam at the bottom of an interval to that at the top of the interval:

Jan 1, 1966

By A. M. Khan, I. Fatt

A new type of steady-state thermal conductivity apparatus is presented and described. Results of measurements on Berea sandstone, Solehofen limestone, Bandera sandstone and Bakelite are presented. The measurements were made in a two-hour period at temperatures up to 80C and axial stress to 125 bars. Results show that the thermoelectric heat pump and tbermistor-controlled electric heater can be used in a steady-state apparatus for measuring thermal conductivity of mck. Steady-state is reached within two hours. With Urethane foam insulation, not more than 3 per cent of heat is lost by radial flaw. The apparatus can be used as a comparator or for absolute thermal conductivity determinations. INTRODUCTION Calculations of heat flow through the earth, both for scientific purposes and engineering design, require thermal-conductivity data on rock. However, the high variability of rock often makes handbook thermal-conductivity data of little value. For engineering purposes in particular, it is desirable to have measured data on samples of the rock mass under study. Laboratory methods for measuring thermal conductivity of rock use either time-dependent1-3 or steady - state temperatures .4,5 The time -dependent (transient) techniques give thermal diffusivity from which thermal conductivity can be calculated if the specific heat and density of the sample is known. Most transient techniques require lengthy and involved sample preparation. Also, the accuracy tends to be low. Steady-state techniques are more exact and are the most commonly used. The steady-state methods, however, also have their difficulties. A long time is often required to reach steady state, troubles have been encountered in establishing good thermal contact between the heating and cooling elements and the rock sample, and unaccounted-for heat losses are difficult to prevent. In this paper we present a description of a new type of steady-state thermal conductivity apparatus. With it, measurements have been made in two hours at temperatures up to 80C and axial stress to 125 bars. Sample preparation is simple and no thermal conductivity standard is needed.6 APPARATUS Our steady-state apparatus is an improvement . over the conduction pile used by Zierfuss and van der vliet4 and Somerton.5 A linear heat flux is established through a conduction pile consisting of the rock sample and a standard material of known thermal conductivity in series. The thermal conductivity of the sample is calculated from the ratio of the temperature drop across the sample

Jan 1, 1966

By M. A. Torcaso, P. Raimondi

To study mass transport in systems simulating oil recovery processes, different porous media were saturated with a mobile (carrier phase) and a stationary phase. Slugs of carrier phase containing a small amount of solute were displaced with pure carrier phase. By analogy to the chromato-graphic processes, the velocity of the solute can be predicted from a knowledge of the partition coefficient and the saturation provided that equilibrium between the two phases exists. Equilibrium was found to exist for different porous media, solutes and rates. The conditions were varied over the range normally encountered in the laboratory and in the field. The longitudinal dispersion of a solute undergoing interphase mass transfer was also investigated. INTRODUCTION The production of hydrocarbons by gas cycling, enriched gas drive and C02 or alcohol displacement involves, among other factors, relative motion between two phases and compounds, hereafter called solute, which are soluble in both phases. The solute is carried forward by the faster flowing phase at a lower velocity than the average velocity of that phase. Retardation of the solute is caused by chromatographic absorption and desorption in the slower flowing phase and by the degree of departure from equilibrium. At equilibrium the concentration of solute in the two phases can be related by the equation* CSw = K Cso ,............(1) where C, and C are the concentration of solute in the aqueous and oleic phases respectively and K is the equilibrium ratio, or partition coefficient. Displacement theories must contain an explicit assumption with regard to equilibrium, i.e., whether the compositions can be related by Eq. 1. The existance of equilibrium depends, in general, on the relative velocity between the phases. Unfortunately, other factors such as gravity segregation and viscous fingering, also depend on velocity. For this reason, whenever effects of rate on displacement were observed,2-4 it was practically impossible to discern what caused them — lack of equilibrium or the factors mentioned above. Equilibrium between phases has been the subject of extensive studies in fields such as extraction or chroma to graphy. It has received only small attention in flow through the type of porous media encountered in oil production.5 For this reason a method was developed which makes it possible to study the movement of a solute as it is affected by rate, type of porous media, partition coefficient and carrier phase, but in the absence of segregation or fingering. The information obtained enables one to determine when the assumption of equilibrium can be made. Briefly, the method consists of (1) saturating the core with a mobile and an immobile phase, (2) injecting a slug made up of the same fluid as the mobile phase and a small concentration of mutually soluble solute, (3) measuring the lag and the peak height of the slug at arrival and (4) correlating these variables with fluid properties such as partition coefficient and mixing constants of the medium. PROPOSED MECHANISM The principles of chromatography are combined with the equation of longitudinal mixing to predict the velocity of a solute slug compared to the bulk velocity and the peak height of a slug.6-9 The equation so obtained is valid under equilibrium conditions only. Therefore, a comparison between experimental and predicted results will give a measure of departure from equilibrium. This work was done with either the oleic or the aqueous phase being immobile. For simplicity, the following development is based on the case where the oleic phase is immobile. However, the treatment is the same in either case.

Jan 1, 1966