Search Documents

By G. M. Webster, G. E. Dawsongrove

The use of percussion-type side-wall cores as an aid in detection and evaluation of hydrocarbon shows and in examination of rock properties has become increasingly popular in recent years. Sidewall coring has an advantage over normal bottom-hole coring methods in that the zones may be surveyed with the electric log prior to selection of sample depths. This allows an efficient sampling of all zones of interest, a procedure not possible with bottom-hole coring techniques unless the costly procedure of continuous coring is utilized. Percussion sidewall coring presently provides the most rapid and lowest cost method of obtaining sidewall cores. It has been recognized for some time that certain discrepancies exist between core analysis data obtained from percussion cores and data obtained from other types of core Samples. This paper presents existing field and laboratory data that were available to the authors. Several examples of comparative core analysis data obtained from various fields are presented along with photomicrographs of thin sections of conventional and percussion cores. The results of this investigation are presented to throw more light on the possible disturbed condition of percussion sidewall cores. SAMPLE DISTURBANCE AND FRACTURING The principal cause of discrepancy between conventional core analysis data and percussion sidewall core analysis data is believed to be the shattering and re-adjustment of the grains due to the large impact shearing forces produced by mechanics of percussion coring. Figs. 1 through 5 are photomicrographs of thin section slides of formation samples obtained by means of diamond cores, conventional cores, and percussion sidewall cores. Figs. 1 and 2 show samples of the

By N. Mungan, J. L. Thompson

Laboratory displacement tests were performed to study oil recovery elficiency by gravity drainage in fractured systems under miscible conditions. The porous media used were cylindrical Berea and Blue Jacket sandstone cores containing a number of well-defined, artificially formed, vertical and subvertical fractures. Butane and Soltrol 130 were the two miscible fluids used. The purpose of the study was to examine the influences of the displacement rate, fracture density, fracture orientation, fracture permeability, matrix permeability, crossflow, core length and connate water on the oil recovery. It was found that displacement rate, matrix permeability and the subvertical fractures affected oil recovery most. The critical flow rate, based on the matrix permeability, was found to be a significant factor in the process. For displacement rates below the critical flow rate, the oil recovery efficiency appeared to be unaffected by the density of the subvertical fractures. At the high displacement rates, the fracture density becomes important, with the recovery being the most efficient in the core having the greatest number of the subvertical fractures. The magnitude of the fracture permeability, the fracture orientation, the core length and the connate water have little effect on the oil recovery efficiency. INTRODUCTION Gravity drainage under miscible conditions from relatively thick reservoirs can be a very efficient recovery process, especially at low flow rates where the gravity forces are dominant and, consequently, the adverse viscous fingers associated with the unfavorable viscosity ratio are minimized or eliminated. Such a process might involve injection of an LPG or some solvent bank at the crest and then driving the bank with dry gas. For some favorable combination of reservoir temperature, pressure and oil composition, there may be no need to inject any solvent bank, since enrichment of the dry gas by the light ends of the crude creates an in-situ solvent bank. Slobod and Howlett2 have studied the effects of gravity segregation in vertical unconsolidated porous media under miscible conditions in the laboratory. The main variables in their study were the viscosity ratio, the density differences and the rate of flow. In this study, the objective was to find the influence of fractures on gravity drainage under miscible conditions. A greater number of related variables were also studied. LABORATORY STUDY CORE DESCRIPTION Displacement tests were performed mostly in cylindrical Berea sandstone cores. In a few cases in which a much lower matrix permeability was desired, Blue Jacket sandstone cores were used. The vertical fractures were formed by cutting the cores with a thin circular blade or by parting them along a bedding plane. The vertical fracture plane always contained the axis of the cylindrical core. The subvertical fractures, cut with a saw, were inclined 45' from the horizontal. The line of intersection between the planes of the vertical and subvertical fractures was parallel to the circular faces of the cores in all cases except one, as depicted by VF3HF-2 in Fig. 1. Fracture geometry of the other cores is also given in Fig. 1. Tables 1 and 2 give the physical properties of the solid and fractured cores, respectively. The suffixes and prefixes shown for each core have been used throughout the paper so that the reader may discern the fracture geometry from the core numbers.

Jan 1, 1970

By Charles R. Stewart, William W. Owens

Reservoir performance predictions based on laboratory core test data assume that fluid flow is laminar for the laboratory test. A study has been made to determine the validity of this assumption for laboratory tests on various types of porosity found in producing limestone formation. Data are presented which show that turbulence and slippage can occur during laboratory tests on hetero-geneous-type porosity limestones, thus causing serious errors in measured single-phase permeabilities and two-phase relative permeability characteristics. In single-phase flow tests it is possible to eliminate turbulence and correct for slippage or to eliminate both factors by controlling test conditions. It is not always possible to control test conditions and thereby eliminate turbulence and slippage in two-phase .flow tests. A correction method is presented which can be used to calculate the true two-phase laminar flow relative permeability characteristic even though furbulence and slippage exist. .INTRODUCTION It is customary to make use of Darcy's law and modifications of this law, together with laboratory data on formation core samples to predict the performance of producing reservoirs. Such predictions are based on an assumption that fluid flow is in the laminar or streamline region for the laboratory test. It was the purpose of this inves- tigation to determine the extent to which turbulent flow may occur in laboratory fluid flow tests on hetero-geneous porosity limestones. Considering that turbulent flow conditions might exist in some laboratory fluid flow tests, additional emphasis was placed on the development of a method to correct for turbulence when laminar flow conditions could not be attained. FLUID FLOW CONCEPTS FOR POROUS MEDIA The Influence of Pore Geometry on Fluid Flow One of the more important factors influencing fluid flow in porous media is the geometry of the pore space which includes such characteristics of the pores as size, shape, distribution, roughness, uniformity, etc. In general, oil- and gas-producing formations can be divided into two broad types on the basis of pore geometry. One has been called sandstone-type porosity media, which is characterized by a small range in pore size, uniformity in shape of the pores, smooth pore surfaces and a regular and uniform distribution of pores. The other type has been called heterogeneous porosity media and is usually limited to the dolomites and limestones. This type is characterized by a wide variation in the size, shape, and distribution of the pores and rough, irregular pore surfaces. It is therefore apparent that conditions are much more favorable for turbulent flow* in heterogeneous-type porosity media than in sandstone-type porosity media. Interrelationship Between Turbulence and Gas Slippage In studying the problem of turbulent flow in laboratory tests on porous media, it is necessary to be aware of the interrelationship between slippage and turbulence for gas flow. As a result of slippage or the Klinken-berg effecta, apparent perrneabilities to gas are greater than the true value because there is no stationary layer of gas in contact with the walls of the flow channels. Gas slippage decreases as the mean free path of the gas molecules decreases. Since the mean free path of any gas decreases with increasing density, increases in static pressure result in lower apparent gas permeabilities. However, a reduction in gas permeability can also be due to turbulence. Therefore, in studying only turbulent flow in porous media, it is necessary to hold gas density, and slippage, constant or to reduce slippage to a negligible value by operating at high static pressures. Presentation of Laminar and Turbulent Flow Data A graphical relationship between permeability and a pseudo-Reynolds number, N,, will be used to show the two types of fluid flow, i.e., laminar and turbulent. The usual graphical method for such a description has been the use of friction factor-Reynolds number charts4. On such a logarithmic diagram, the laminar region appears as a straight line having a slope of 45 degrees. As the friction factor decreases and the Reynolds number increases, the turbulent region is reached and appears as a deviation from the 45-degree slope line. However, in petroleum engineering literature resistance of por-

By L. L. Handy

The most common method of identifying hydrocarbon-bearing strata in a well that penetrates many different formations involves measurement and interpretation of the electrical properties of the formations as determined by electrical logs. Even though this method is used extensively, and even though in a great many instances it is capable of indicating presence of oil or gas, situations arise for which it is extremely difficult, if not impossible, to deduce the presence of hydrocarbons. These situations may involve the following. 1. A thin formation, bounded by highly resistive formations, in which it is impossible to obtain the actual resistivity of the uninvaded zone with existing logging devices. 2. A formation in which invasion has been so extensive that a value for the uninvaded zone resistivity cannot be obtained. 3. A very shaly formation in which the resistivity index, I, is lower than that usually associated with productive formations. 4. Laminated formations comprised of thin productive sands separated by thin shale streaks in which the individual sand and shale streaks are too thin to permit measurement of uninvaded-zone resistivity with existing logging devices. 5. Productive formations in which the water saturation is high. To extend the utility of electric log interpretation to identification of hydrocarbons in all types of formations, there is strong incentive to find a method not subject to these limitations. Some time ago, in connection with research on the wettability of reservoir rock, an investigation was conducted in which the resistivities of cores were measured shortly after they were removed from a core barrel, and again after they had been extracted and restored to their original oil and brine saturation.' The resistivities after extraction were generally lower. Other tests made on the cores indicated that they were more nearly water wet after they were extracted; thus, it was assumed that the observed changes in resistivities were due to a change in wettability of the cores. Other experiments'," have shown that resistivities of rock samples are sharply dependent on wettability. These experiments have shown that oil-wet samples are more resistive than water-wet samples. To obtain an understanding of how the wetting properties of the surfaces of core material affect electrical resistivity, a series of experiments was conducted. Two groups of core samples were prepared for testing. One group contained brine, but no residual oil. The other group was saturated with brine, flooded with oil to a low water saturation, then flooded with brine to a final residual oil saturation. Resistivity measurements were made on each group. Both groups were then flooded with the original brine to which a chemical had been added that renders sand and clay surfaces preferentially oil wet, a so-called reverse-wetting agent. very little change in resistivity was observed in cores containing only water. The group containing residual oil, however, showed resistivity increases of 100 to 200 per cent. These experiments showed that the resistivity of a core containing oil could be altered by changing wet-tability of the core. Moreover, the possibility was introduced that reverse-wetting agents might be employed as the basis for a logging method for identification of oil-bearing strata. Since behavior of a porous rock containing gas and water might be expected to be similar to that of a rock containing oil and water, such a method should also be applicable to identification of gas-bearing zones. In principle the wettability of the invaded zone could be reversed without altering conductivity of the interstitial water or the hydrocarbon saturation therein. Those strata showing significantly increased invaded-zone resistivities would, therefore, contain hydrocarbons; those with no significant change would be filled only with water. Addition of a reverse-wetting agent to a hydro carbon -bearing zone which is, by nature, already preferentially oil wet would not result in an enhancement of its resistivity. It is generally believed, however, that most hydrocarbon-bearing strata are preferentially water wet.

By R. C. Hassinger, D. V. von Rosenberg

Transverse dispersion has received considerably less treatment in the literature than has longitudinal dispersion. Different methods for determining transverse dispersion coefficients have been used in different investigations, and the results obtained have not been consistent enough to permit accurate generalizations as to the effect of various physical parameters on the magnitude of these coefficients. A numerical solution to the differential equation describing transverse dispersion in the absence of longitudinal dispersion was obtained to enable one to calculate the dispersion coefficient from experimental results. The more general dispersion equation including longitudinal dispersion also was solved numerically to give quantitative limits of a dimensionless group within which the assumption of negligible longitudinal dispersion is justified. Possible experimental procedures were examined, and one utilizing a cylindrical packed column was chosen for the determination of transverse dispersion coefficients. Values of these coefficients were determined for a system of two miscible organic fluids of equal density and viscosity, for two sizes of packing material over a wide range of flow rates in the laminar regime. The dispersion coefficient was found to decrease, for a constant value of the product of packing size and interstitial velocity, as the size of the packing material particles increased. INTRODUCTION Longitudinal dispersion has received extensive treatment in the literature, and consequently is better understood than its orthogonal counterpart, transverse dispersion. Many mathematical models of dispersion processes assume that transverse dispersion is rapid enough to damp out any radial** concentration gradients and therefore may be neglected. Laboratory and production results, however, indicate that this is a poor assumption.l,2 Various experimental procedures for determining transverse dispersion coefficients have been used in previous investigations, but the results have generally been expressed by similar correlations. The transverse dispersion coefficients obtained, however, have often varied considerably for given values of the correlation parameters. We feel that further experimental determinations of transverse dispersion coefficients will help alleviate some of the inconsistencies in these empirical correlations. One assumption implicit in all previous investigations is that of negligible longitudinal dispersion in the experimental system. An attempt to justify this assumption often is made using intuitive reasoning, but it is apparent that this reasoning must break down as the condition of zero flow rate is approached. A mathematical examination of the equations describing the system yields physical limits outside of which the assumption of negligible longitudinal dispersion is invalid. BACKGROUND In a porous medium, the "effective molecular diffusivity" De is less than the molecular diffusivity D measured in the absence of a porous medium, due to the tortuous path which a diffusing molecule must travel. Various authors3-5 have reported values of the ratio De/D in the range of 0.6 to 0.7. When there is fluid flow within the porous medium, mass transfer occurs by convective dispersion as well as by molecular diffusion. These are separate phenomena and can be treated as such on a microscopic scale. However, the mathematical complexity is such that only extremely simple geometries could be considered, and the results hardly would be applicable to the complex geometries existent in actual porous media. A more

Jan 1, 1969

By L. H. Robinson

Triaxial compression tests have been performed to determine the strength characteristics of limestone, sand-stone and shale rocks subjected to controlled stress conditions. This control tuns exercised by varying the liquid pressures within and around a plastic-encased rock specimen. The pressure in the pores of the rock was varied throughout the range from atmospheric to 15,000 psig; the external pressure was changed over the same range with various positive pressure differences between is and the intertnal pressure. The data show that the rock strength increased and the mode of failure changed as the pressuire. surrounding the rock became greater than the pressure in the pores of the rock. These observations and the resu1ts of microbit drilling experiments indicate that the increased rock strength under pressure may bc an important effect in reducing drilling rare, but that other factors are probably of even greater importance. INTRODUCTION A large part of the research conductcd to reducc thc cost of drilling is directed toward improving present drilling technology and dcveloping new drilling methods. In spite of the fact that drilling is essentiaIly the act of making rocks fail, many details of the mechanism of failure are not known. It seems probable, therefore, that continuing studies of the fundamentals of rock failure can lay a foundation for future improvements in drilling technology. One of the factors which affects rock failure is the stress applied to the rock. This effect has been known for many years.'.'.; As an example, stress studics have been published by engineers' concerned with the strength characteristics of concrete in large dams. More recently, a number of investigators4-12 have explored the effects of stress on the strength of geologic formations and the geologic implications of their findings. During the last few years some research'" on rock strengths has also been conductcd in an attempt to gain additional insight into the factors which affeet rock drillability. In detcrmining the strength characteristics of rocks under different stress conditions, investigators have conventionally used triaxial compression equipment. wherein a jacketed rock cylinder is uniformly loaded from all directions and then compressed longitudinally. The results of this work have demonstrated that the mode of failure (brittle or malleable) is dependent upon the loading pressure, and that the strength of the rock increases as the loading pressure increases. These results clearly manifest that the stress, which is determined by the loading pressure, affects rock strength, and this implies that rocks under stress should be stronger and hence harder to drill. During the drilling of a porous rock the stress conditions are determined not only by the pressure surrounding the portion of the rock undergoing failure (confining pressure) but also by the pressure of the interstitial fluids within the rock (pore pressure). The importancc of both these pressures on the drilling rate of rocks has already been demonstrated by a number of investigators. Murray and Cunningham" found in microbit drilling expcriments that at constant pore pressures the drilling ratc decreased as the hydrostatic pressure surrounding an unjacketed, impermeable rock increased. Later Eckel" found that for jacketed Iimestone specimens the differcnce between hydrostatic (wellbore) pressure and formation (pore) pressure had an important effect on drilling rate. These observations of the effect of pressure on drilling rate emphasize that for drilling studies triaxial compression tests should be conducted under conditions of confining and pore pressure which simulate conditions underground. In the past, triaxial tests have been made under conditions which do not reproduce those underground. In particular, most of the tests on porous materials have keen made at varying confining pressures but with atmospheric pore pressures. In only one group of tests' was any data on the effect of pore pressurc reported, and these few data were on materials not representiltive of those encountered in drilling earth formations. The research described herein was initiated to deter-mine the importance of both internal pore pressure and confining pressure on failure characteristics of limestone. sandstone and shales. Triaxial loading equipment was used to obtain strength data at pressures up to 15,000 psig. An additional objective of this research was to apply these results to drilling rate studies by comparing them with the previously published drilling data of Murray and Cunningham" and of Eckel12. The work reported here is but one step in determining and evaluating the fundamental mechanisms which control drilling rate. Additional research is needed to clarify other aspects of the drilling process. The better

By J. C. Melrose, W. R. Foster, J. G. Savins, E. R. Parish

Many differences can be imagined between gas-oil flow in which the gas is supplied at the face of the core and gas-oil flow in which the flowing gas was originally dissolved in the oil. If capillary pressure characteristics and flow requirements control gas saturation distribution, the gas would be expected to be located at preferred sites within the porous medium as determined by pore sizes. On the other hand, during solution gas drive the gas first appears as bubbles through a nuclea-tion process. Nothing in self-nucleation theory specifies at which sites the first bubbles should be formed. In all probability they will be randomly distributed throughout the porous medium. Furthermore, it is not at all certain that even at low rates of production the gas will redistribute itself after nucleation to the channels normally occupied by gas in simple gas flow. Stewart, et at, have shown that at least for some limestone samples, oil recoveries could not be predicted for all rates of production using any one set of relative gas and oil permeabilities. An important factor in controlling recoveries during solution gas drive was the rate of bubble formation, higher rates giving higher recoveries. Stewart, et al, attributed the increase in recovery to a better distribution of the gas phase in heterogeneous limestone samples than is obtained by simple external gas drive. Differences in recovery from these causes were not reported for sandstone cores. In the experiments to be reported here, oil recovery, pressure and producing GOR history were measured during solution gas drive for a 5-ft sandstone core. The results were compared with predictions from the Muskat method for computing solution gas-drive behavior using external gas-drive relative permeability. The effects of changing the rate of production and oil viscosity were studied. At high laboratory rates of average pressure decline, two observations were made which would not have been predicted by Muskat's depletion theory: (1) oil recovery increased with increasing rate of production for a given viscosity oil, and (2) oil recovery increased with increasing oil viscosity for a given high rate of production. Both of these observations are explained as consequences of diffusion control of gas saturations superimposed on the normal gas-oil flow requirements, Fur- thermore, discontinuous gas phase flow appears to be significant during solution gas drive. The laboratory tests were performed at rates of average pressure decline many times greater than the maximum possible rate of average pressure decline in an actual oil field. It is, therefore, not possible to draw any direct conclusions regarding the effect of rate on recovery for the solution gas-drive mechanism under actual field conditions. However, at the lower laboratory rates, recoveries were nearly independent of rate and could be predicted by the Muskat method, using external gas-drive relative permeability data. These results suggest that at normal oilfield rates the effect of rate on recovery for the solution gas-drive mechanism is negligible. EXPERIMENTAL PROCEDURES The core material for the pressure depletion studies was Bandera sandstone from an outcrop in the Mid-Continent. This sandstone was selected because of its low permeability (about 10 md), which would permit the development of substantial pressure gradients in the corn at moderate flow rates. The core was 5-ft long and 2-in. in diameter. Its properties are listed in Table 1. Relative gas-to-oil permeability ratios were measured by an external gas-drive method2. The results are shown for a short 2-in. core and for the 5-ft core in Fig. 8. Oils used in the pressure depletion experiments were kerosene and a highly refined white oil (standard white oil No. 3) with gas-free viscosities of 1.8 and 25 cp, respectively. The gas was a naturally occurring methane from Gough field, Inglewood, Calif. The oil viscosities, gas solubilities and formation volume factors are plotted as functions of pressure at 75°F in Figs. 1 and 2. Methane viscosities and compressibilities were obtained from the literature"'. A core mounting was required which could withstand up to 2,500 psi internal pressure. This was obtained by first encasing the core completely in a plastic resin (Scotch Cast, manufactured by Minnesota Mining & Manufacturing Co.) The plastic covered core was then inserted into a steel pipe equipped with screw caps so that the plastic coating could be pressured from

By P. van Meurs, C. van der Poel

From observations of the linear displacement of oil by water from a porous medium as visualized in transparent models, new insight into the mechanism of vircous fingering, as occurring in the case of unfavorable oil-water viscosity ratios, has been obtained. A mathematical description of the observed phenomena leads to formulas. expressing both oil production and pressure drop across the formation as a func-tion of cumulative water injection with the oil-water viscosity ratio as a parameter. Theoretical results are compared with those obtained from scaled laboratory experiments in sand-filled rubes in which water was injected at the one end and production was obtained from the other end. In the second part of the paper the same principles are used to calculate a more realistic production scheme in which production of oil and water is obtained from wells regularly distributed over the field. Results show satisfactory agreement with production data of an actual oil field. Application of the theory is extremely simple as only a limited number of readily available parameters have to be known and results are presented in an analytical form. INTRODUCTION When oil is displaced from a porous medium by water having a viscosity lower than the oil. the oil-water interface is essentially unstable and has the tendency to break up into what are called "fingers" or "streamers". This Viscous fingering has been noted before by Engelberts and Klinkenberg' and has recently been demonstrated in transparent models by van Meurs-. It is highly questionable whether it is at all possible to perform relative permeability measurements In the laboratory reproducing the same distribution of water and oil in the cores as Present in practice. Therefore, we based our descriptive theory on an idealized model which was "tailored" to the fingering mechanism as observed in transparent models under scaled conditions. This approach was the more attractive, as a highly simplified model had already yielded an adequate description of the process. Though we are aware our theoretical approach can be extended in several ways, it is given here in its simplest form to explain the fundamental principles and to avoid complicated formulas. SIMPLEST CASE OF A LINEAR DRIVE Observations and Theoretical Model In Fig. 1 a series of photographs shows how water penetrates into a formation saturated with an oil which is considerably more viscous than the water. The technique for obtaining these photographs, using transparent models, has already been described?. The fingers can be clearly observed and are seen to be mainly parallel to the flow direction. Note that the contours of these fingers are not smooth but show protrusions of water enclosing oil pockets.

By D. R. Parrish, W. E. Lamoreaux, W. W. Owens

Several methods are now being used by the industry for determining the gas-oil flow characteristics of reservoir rock samples. Most of the laboratory experirnerltal rlzethods can be classified either as steady-state, unsteady-state, or Hassler-type tests. Although the validity of the results obtained by any specific method is not questiorzed, in general, the complexity of these methods leaves much to be desired. This investigation dealt primarily with the development of a simplified unsteady-state (gas drive) laboratory test method for measuring kg/ko re1ationships of small core samples. A procedure for selecting uniform and representative sandstone rock samples was developetl since it was observed that the reliability of the test resu1ts was dependent upon having a sample which was uniform in porosity and permeability. The importance of factors known to influence the results of this type test, such as oil viscosity, gas expansion, and pressure differential were eva/uated. The method of calculating the test results was found to give valid data over a wide range of gas expansion. The limits of pressure differential and oil viscosity which would insure reliable results were defined. In cornparisorz with other presently used testing methods, the gar drive method requires considerably less apparatus, is less time consuming, and is more readily adaptable to a routine type core test. INTRODUCTION The detcrmination of the gas-oil flow characteristics of reservoir rocks has been given concentrated study by the oil industry for the past decade. Both analytical and experimental methods have been considered. The analytical methods have generally relied on a simplified characterization of porous materials in order that the fundamental laws of fluid flow and capillarity could be utilized. More direct experimental methods have been proposed by numerous production research laboratories with the primary goal being a simple and inexpensive, yet reliable method for measuring the gas-oil flow characteristics of reservoir rock samples. The experimental methods which have been reported in the literature generally can be grouped into one of the following classifications: 1. Steady state1,2,3,4 in which a fixed ratio of flowing fluids is forced through the test sample until equilibrium conditions of saturation and pressure gradient are established. 2. Hassler or stationary liquid method5,6,7 in which a uniform liquid saturation is attained, usually through capillary desaturation, then gas or oil is flowed through the sample at a pressure drop insufficient to disturb the pre-established saturation. :. Non-steady state method1,8,9,10 In which only the gas phase is injected into the sample, causing displacement of the oil phase and resulting in a system of continually changing average saturation and saturation gradient. The steady-state and Hassler methods have been most popular because, no doubt, of the simplicity of the actual relative permeability calculations. The problems involved in these methods were limited to apparatus design and measurement techniques. In the steady-state method, for instance, a constant flowing fluid ratio is established which, at equilibrium, results in an essentially uniform saturation over the entire length of the test sample. Therefore, the basic flow characteristics could be obtained by correlating the average saturation of the test sample with the kg/ko determined at various flowing gas-oil ratios.

Jan 1, 1957

By J. G. Richardson, J. W. Graham

The purpose of this paper is to present the results of theoretical and experimental studies of water imbibition. The imbibition processes are involved in recovery of oil from stratified and fractured-matrix formations in natural water drives and water flooding. An understanding of the role of inhibition in implementing the recovery of oil from such formations is deemed essential to proper control of these reservoirs to achieve maximum recovery. The theoretical studies involved development of the differential equations which describe the spontaneous imbibition of water by an oil-saturated rock. The dependence of the rate of water intake by the rock on the permeability, interfacial tension, contact angles, fluid viscosities and fluid saturatiorls is discussed. A few experiments were performed using core samples to determine the effects of core length and presence of a free gas suturation. The role of water imbibition in recovery of oil from a fractured-matrix reservoir by water flooding was investigated by use of a laboratory model. This model was scaled to represent one element of a frac-tured-matrix formation. Water floods were made at various rates with several fracture widths. Interpretations were made of the behavior expected in a system containing many matrix blocks. The presence of a free gas sntu.ration was found to reduce the rate of water imbibition. In the reservoir prototype of the fractured-matrix model, water imbibition rather than direct displacement by water was the dominant mechanism in the recovery of oil at low rates. INTRODUCTION Imbibition may be defined as the spontaneous taking up of a liquid by a porous solid. The spontaneous process of imbibition occurs when the fuid-filled solid is immersed or brought in contact with another fluid which preferentially wets the solid. In the process of wetting and flowing into the solid, the imbibing fluid displaces the non-wetting resident fluid. Common examples of this phenomenon are dry bricks soaking up water and expelling air, a blotter soaking up ink and expelling air and reservoir rock soaking up water and expelling oil. As increasingly better lithological descriptions have been made of the characteristics of petroleum-bearing formations, it has become obvious that imbibition phenomena which were once considered laboratory curiosities are of practical importance. For instance, in reservoirs composed of water-wet sand strata of different permeability in intimate contact, the tendency of water to channel through the more permeable stratum is offset by the tendency for water to imbibe into the tight sand and expel oil into the coarse sand. Also, in fractured-matrix formations the tendency of water to channel through the fractures is offset by water-wet matrix blocks. As some imbibition of the water into the of the largest fields in the world are fractured-matrix reservoirs, it has become increasingly important to understand all the factors involved in the imbibition process. Examples of fractured-matrix reservoirs are the Spraberry field in West Texas which produces from a fractured sandstone', the giant Kirkuk field in Iran', the Dukhan field in Qatar, Persian Gulf2, and the Masjid-I-Sula-main and the Haft-Kel fields in Southwestern Iran, which produce from fissured limestone3. Research into recovery of oil from fractured-matrix formations was stimulated by the rapid decline of oil productivity of wells in the Spraberry formation. One result of this research was the water imbibition process developed by the Atlantic Refining Co.4 Another idea was that much of the Spraberry oil could be recovered by conventional water-flooding procedures5. Subsequently, pilot floods were conducted in this field to test the feasibility of these ideas. It was felt that an understanding of the role played by imbibition processes in displacement of oil from a fractured-matrix reservoir could not be obtained from field data alone because of the many complicating factors and uncertainties involved. Therefore, theoretical and laboratory studies were undertaken to provide this understanding. Study of the equations which describe the linear, countercurrent imbibition process provided an insight into the role of various factors in the process, such as the permeability of rock and inter-facial tension. In addition to the theoretical studies, imbibition experiments were conducted with core samples to determine the effect on the rate of imbibition of such variables as core length and free gas saturation. The principal experimental studies were conducted by water flooding a scaled model of an clement of a frac-tu red-matrix reservoir to evaluate

By C. S. Brooks

Several of the commonly used methods for deter-mination of effective grain densities of minerals have been evaluated for accuracy and precision, and their applicability to various types of petroleum reservoir minerals tested. Consideration was given to the effect of consolidation, swelling clay, and wettability on these grain-density determinations. Procedures with helium displacement, a modified low-pressure Boyle's law porosimeter, and liquid displacement (water, n-decane, and n-octanol) were tested. The helium-displacement procedure was used primarily as an independent reference for which the effects of fluid-solid interaction, e.g., adsorption and wettability, were negligible. The procedures with the Boyle's law porosimeter and with liquid displacement are the methods in common use. Minerals investigated included calcite, West Texas limestone core samples, California sandstone core samples, and Wyoming bentonite. The number of samples tested, while not large, was adequate to permit application of statistical methods. Grain-density measurements on calcite (Iceland spar) provided a basis of comparison for the accuracy and precision of the procedures under investigation for a material with a high degree of chemical and physical homogeneity and for which fluid-solid interactions were at a minimum. Liquid-displacement procedures (water or n-decane) provided the greatest reliability and precision for measurement of the grain densities of the limestone tested. Gas-displacement procedures (helium or Boyle's law porosimeter) provided the greatest reliability and experimental convenience, if not the greatest precision, for the measurement of the grain densities of the swell- ing clay, Wyoming bentonite, and the shah sandstones which -contained appreciable amounts of swelling, in-terstitial clays. These experiments were designed to evaluate the applicability of the grain-density procedures in common we to some of the less tractable minerals—tight limestones, shaly sandstones, and swelling clays. The results should provide a useful basis for a selection of the method most suitable for application to the various petroleum reservoir minerals encountered and especially for use in the estimation of formation porosity. INTRODUCTION These experiments were undertaken to evaluate the accuracy and precision of several commonly used methods for determination of effective grain densities of minerals, and to test their applicability to various types of petroleum reservoir minerals. Consideration was also given to the effect of consolidation, clay, and wettability on these grain-density determinations. Procedures used included measurements with helium displacement, with a modified low-pressure Boyle's law porosimeter, and with liquid displacement (water, n-decane, and n-octanol). The minerals investigated included calcite, West Texas limestone core samples, California sandstone core samples, and Wyoming bentonite. The accurate determination of the porosity of petroleum reservoir formations, in particular those of low porosity, has received much attention. Both bulk-volume and grain-volume measurements necessary in the evaluation of porosity have been examined by many investigators with the objectives of improving both the accuracy and the precision of these determinations. In this investigation statistical criteria for selection of the most reliable and suitable procedure for use with a given type of mineral are demonstrated. The grain-volume measurement procedures appear to be the most susceptible to improvement.

Jan 1, 1958

By P. E. Baker

An experimental study of heat flow in steam flooding was carried out with steam displacing water in a plane-radial fluid-flow model. Temperature distributions in the model reservoir, overburden, and substratum were measured at frequent intervals during flooding by means of a fixed array of thermocouples. Heat losses to overburden and substratum and the division of heat between the steam and hot-water zones in the reservoir were calculated from the temperature data by performing numerical integrations over the appropriate volumes. Temperature profiles and heat distribution results are presented. Some results were (1) the fraction of injected heat that is lost to overburden and substratum, when expressed as a function of time, did not depend on injection rate; (2) a significant portion of the injected heat was contained in the hot-water zone ahead of and underlying the steam front in the reservoir. The division of heat between the steam and hot-water zones did depend on injection rate; and (3) considerable gravity override was observed. The experimental results are compared with those of several published theoretical studies of heat transfer in hot-fluid injection; approximate agreement on heat losses is found with some of the most idealized theories. A valid comparison cannot be made in regard to steam us hot-water zone; basic interpretation of the hot-water zone, in terms of a vertical steam front, is clouded by gravity override in the experiments. INTRODUCTION The experimental study reported here primarily concerns two problems of heat transfer in steam flooding that have economic and engineering importance. These are (1) the fraction of injected heat lost to the overburden and substratum, and (2) the division of the remaining heat between a zone containing steam and a zone of hot water at lower-than-boiling temperature. The volume of the steamed region, which has a direct effect on oil recovery, can be calculated readily from these results. Several valuable theoretical studies have been published on this problem, but no experimental tests of the theories have previously appeared. Tests were needed because the effects of the necessary simplifications in the theories could not be determined otherwise. The basic, untested theoretical problem is the effect of "decoupling" heat and mass transfer, and the resulting inability to include the effects of heat transfer and condensation on fluid flow. The theoretical problem solved so far has been a heat transfer problem in which fluid is present as a heat carrier and flows in a stipulated manner, either linear or plane-radial. The simplifications that have sometimes been made in regard to thermal conductivity by, for example, Marx and Langenheim,l Lauwerier2 and Thomas3 apparently do not have an important effect on heat loss calculations, although the results of the different methods disagree to a certain extent, as shown by Spillette4 and by Ramey.5 A recent paper by the author6 compared the results of experiments in hot-water flooding with several theoretical results and showed a rather poor agreement. It will be seen that the experimental heat loss results for steam flooding, however, show acceptable agreement. These comparisons are discussed in a later section. EXPERIMENTAL APPARATUS AND PROCEDURE The system studied was one in which a sand pack, initially water saturated and at a uniform temperature, was flooded with steam at a constant rate and temperature. The temperature distribution in the "reservoir" and adjacent "strata" was measured with an array of thermocouples, and heat losses were determined from these measurements. The sand pack contained no oil because the only concern of the experiments was heat flow, and it was important to have the flow pattern as simple

Jan 1, 1970

By N. L. Hancox

Previous researchers have used the cutting and rejoining model to represent a porous body when attempting to place Archie's law on a theoretical basis and calculate rock permeabilities (F = Ø-m, where F = formation factor, = - porosity and m = a constant). In using the model, a porous solid is assumed to consist of sections with parallel sides containing straight. continuous pores of total porosity 4 (pore to total volume or, for a homogeneous body, pore to total area). The joining of pores at an interface between two sections is assumed to be random which results in a common pore area Ø2. This reduced area Ø2 is then taken as that area which limits electrical conduction or fluid flow in the body. Agreement between experimental results and results calculated using a model is not always good Because of this, various modifications of the model are presented — for instance, the random nature of the joining is questioned. Formation factor measurements on specially constructed porous specimens are described in this paper. it is concluded that the common pore area limits flow completely only if sections containing straight pores are very thin, and that, in most natural materials, the common pore area at an interface must be greater than that given by random joining. However, the model is still useful in permeability calculations, provided that the fractional area limiting flow is determined from formation factor measurements. INTRODUCTION Two important properties of porous media are permeability k, which describes the ease with which a fluid can flow through a porous medium, and formation factor F, defined as the ratio of the resistivity of a saturated specimen to the resistivity of the saturating liquid which describes the electrical properties of the saturated material. k is defined by Darcy's law,l k= µqL/P Ae . . . . . . . . . . . . . . . [1] where µ = liquid viscosity, q = flow rate, L = length of specimen, p = pressure across the specimen and = - the effective area over which flow occurs, F is given experimentally by Archie's law,2 F = Ø-m ,..............(2) where 4 is porosity defined as the pore to total volume or, for a homogeneous body, the pore to total area; and m lies between 1.3 for loose sand and 2 for cemented sand. Consider a porous body of total cross-sectional area A and length L, saturated with a liquid of resistivity R and with an effective area for conduction of Ae. Then, from the definition of F given previously, F=A/Ae , ................(3) or, from Eqs. 2 and 3, Ae = A Øm ..............(4) Thus, in any attempt to calculate k or F, some assumptions must be made concerning the pore structure of the body and a value derived for A,. To determine Ae, Childs and George3 and Wyllie and Gardner4 suggested an elementary statistical representation often referred to as the cutting and rejoining model of a porous body. The body is assumed to be composed of two or more parallel-sided sections of porosity 4, containing straight, continuous pores with a common pore area of Ø2 at an interface between two sections. The value is given by a random joining of the pores in either section. Two such sections,

Jan 1, 1969

By E. F. Johnson, V. O. Naumann, D. P. Bossler

A method is presented for calculating individual gas and oil or water and oil relative permeabilities from data obtained during a gas drive or a waterflood experiment performed on a linear porous body. The method has been tested and found both rapid and reliable for normal-sized core samples. INTRODUCTION Individual oil and gas or oil and water relative permeabilities are required for a number of reservoir engineering applications. Chief among these is the evaluation of oil displacement under conditions where gravitational effects are significant, such as a water drive or crestal gas injection in a steeply dipping oil reservoir. Numerous proposed methods of obtaining relative permeability data on reservoir core samples have been too tedious and time consuming for practical use, or have yielded questionable and sometimes inconsistent results. A method has been developed by which the individual relative permeability curves can be calculated from data collected during a displacement test. The method is based on sound theoretical considerations. Using this method, with a properly designed experimental procedure, relative permeability curves can be obtained using core samples of normal size (i.e., 2 to 3 in. in length and 1 to 2 in. in diameter) within a few days after receipt of the core. In a recent publication D. A. Efros1 describes an approach to the calculation of individual relative permeabilities that is based on the same theoretical considerations. We believe the approach described in the present paper is more adaptable to practical application than the method implied by Efros. In addition, comparisons with independently determined relative permeabilities are furnished to substantiate the reliability of the new method. DERIVATION Previously the theory of Buckley and Leverett as extended by Welge2 as been used to calculate the ratio of relative permeabilities. In the derivation which follows, this theory is further extended to permit calculation of the individual relative permeabilities. The theory assumes two conditions which must be achieved before the method is applicable. They are that the flow velocity be high enough to achieve what has been termed stabilized displacement,3 and that the flow velocity is constant at all cross sections of the linear porous body. In stabilized displacement the flowing pressure gradient is high compared with the capillary pressure difference between the flowing phases. The high pressure drop insures that the portion of the core in which capillary effects predominate will be compressed to a negligibly small fraction of the total pore space. The assumption of constant flow velocity at all

By R. E. Gilchrist, R. F. Nielsen

When a gas is displaced by mother in a porous medium, and there is a relative immobile liquid present. there is a transition zone in which the gas composition varies from essentially that of the original gas to that of the injected gas. The change in extent of this zone with time may be governed by a mixing process in the direction of flow and by the rate of establishment of gas-liquid equilibrium, if the gases are appreciably soluble in the liquid. The relative importance of these two factors in deterri~ining the growth of the transition zone is discussed, with appropriate equations and experimental data. INTRODUCTION The various solution methods of oil recovery involve a transfer of the components of a mixture, whether within a single phase or between a liquid and gas. The similarity of a natural porous medium to systems used in separation and chromatographic processes makes it seem reasonable that certain concepts used by chemical engineers in "mass transfer" might be applicable to petroleum reservoirs. One approach, the HETP (height equivalent to a theoretical plate), involves a distance or length of time over which equilibrium might be con. sidered as established. Another approach involves direct use of a rate of mixing or rate of establishment of equilibrium. The latter approach has been preferred in relation to production problems. DISCUSSION AND PROCEDURE When a gas or liquid in a porous medium is displaced by one of different composition but miscible with it there will be set up a transition zone in which the composition varies from cssentially that of the original fluid to essentially that of the injected fluid. If the fluid being displaced is a gas, and there is a relatively immobile liquid present, the transition zone may be

By R. I. Lahring, C. G. Inks

A surfactant waterflood, using a block copolymer of ethylene and propylene oxide (PLURONIC L64*), wa., started in 1956 in Kansas to determine the commercial feasibility of such a process. The production formation was the Bartlesville sand at a depth of 1,400 ft. Leares in the test area had attained both primary and secondary production peaks, and were on a normal decline when the tests were initiated. The injection wells were taking water at a reasonable pressure and none of the operating conditions were particularly unusual. The leases were selected primarily because of normal operating conditions, and not because of any serious problems demanding remedie>,. Oil production and water injection data were recorded during the 10-year test period, and their compilation indicates that the surfactant resulted in an increase of about 9 percent in secondary oil production. Other benefits derived from the surfactant injection were lower injection pressure that resulted in reduced power consumption and reduced erosion of the pumps; and lower calcium scale buildup. These fringe benefits more than offset the cost of the surfactant used. Introduction During the past several decades, considerable research has been carried out on secondary recovery of crude oil that remains in the producing formation despite the most efficient, current primary production practices. Methods are constantly sought of increasing recovery of oil, while at the same time improving the economic efficiency of operations. One method that has received intensive study over the years is the use of a surfactant as a waterflood additive. Laboratory tests show conclusively that dilute solutions of surfactants remove considerably more oil from uncon-solidated sand cores than does water alone.'- The loss of surfactant by adsorption at solid-liquid interfaces has caused serious doubt that surfactants can be used economically as waterflooding additives. This is especially true with anionic and cationic surfactants, and to a lesser extent with nonionic surfactants. This work was undertaken to determine the effectiveness in increasing oil production by injecting water containing small quantities of a nonionic surfactant. The surfactant selected was PLURONIC L64*, a block copolymer based on ethylene and propylene oxide, having the following structure The selection of the test site, a typical shoe-string sandstone type of field located at Batesville pool, Toronto District of Woodson County, Kans., was made because of its normal primary and secondary recovery operating conditions, and not because of any serious operating problems. When the tests were begun, the field was experiencing the customary post-peak production decline of a waterflood. Using for the tests a field that was on a normal production decline, while maintaining constant all variables other than the addition of the surfactant. was of varamount importance if success or failure was to be attributed to a particular cause. It was reasoned that if a surfactant was added to the injection waters and nothing else was done to influence the oil production or injection of water, the effects of the surfactant on the flood could be determined precisely. Production History The first well in the tested field was drilled in 1935. Gas repressurization began immediately, but it was not until Aug., 1938, that the major portion of the field was developed. During the period of repressurization, 1,700 MMcf of gas was injected. As can be seen from the curves of Fig. 1, the production peak from gas injection was attained in 1937 and a normal production decline followed. Early in 1942, water injection wells were drilled and injection of water Was initiated. The amount of water injected increased as more injection wells were drilled until a peak injection rate of 2,700 BID was reached in

Jan 1, 1969

By O. C. Holbrook, George G. Bernard

A new theoretical treatment has been obtained for the behavior of pattern waterflood injection wells when closed in. Two cases are treated: Case I where oil and water are assumed to have the same properties, and Case 2 where they arc different. In applying the method, one plots log (p — p,) vs closed-in time, where p is well-bore pressure at any tims and p, is static pressure. The value of p. is determined by trial and error as that value which makes the plot linear at large time. A value for the permeability-thickness product can be determined from the intercept of this linear part, and a value of the skin factor from the injection pressure at time of closing in. Application of the method to data from water floods at three fields seems to give reasonable results. For the case of unit mobility ratio, it is proved that this new method should give the same value for permeability-thickness product as the conventional pressure build-up method. In addition, the new method gives correct values for static pressure, whereas the conventional method does not, often indicating negative static pressures. The new method may be used in cases where the surface pressure persists after closing in as well as in cases where it does not. INTRODUCTION It is of considerable interest and importance to be able to determine the characteristics of the reservoir in an area surrounding a water injection well. Thus, if we can determine early in the life of an injection well that there is a considerable "skin effect", remedial measures can be started before a full-scale pattern flood begins. Similarly, if it can be shown that a gradual buildup of skin effect is occurring with time, measures to free the water of plugging material can be taken. Determination of static pressure in the water-injection well may show that the water is entering a thief zone and not the desired reservoir. Finally, determination of the permeability of the sand around the injection well will allow estimation of the future relation between injection pressure and rate. It should be possible to determine average reservoir permeability, skin effect and static pressure from pressure fall-off data. However, at the time we began work on this subject, it was thought that no adequate theory on which to base such determinations' was available. According to the conventional method which considers the reservoir to be filled with one fluid of small compressibility (see Van Everdingen, Joers2, and Nowak2), shut-in pressure is plotted vs log where is injection time, and At is closed-in time. The physical significance of injection time, may well be questioned in this case, since in a reservoir completely filled with a single fluid (as required by this theory) and with input and output rates equal, the pressure behavior after an initial transient is independent of t,. Attempts by our Tulsa area to use this theory led to negative values of static pressure in most cases. Because of these limitations of the method discussed above, it was decided to attempt to develop a new theory of pressure behavior in water injection wells, one which would apply when there is a gas saturation, as is so often the case in water floods. In the following treatment the assumptions and basic equations are given first, then the method of application of the equations. A complete example is given to clarify details of application. All difficult mathematics has been placed in the appendices so that the reader can follow the text without difficulty. However, if he wishes only to apply the results without knowing the basis for them, he can learn how to do this from reading only the sections entitled "Plotting of Experimental Results" and "Example." ASSUMPTIONS AND BASIC EQUATIONS Statement of Problem It will be assumed that a horizontal layer of constant thickness contains in its pore system a mixture of oil, gas and water. While water is being injected into this pore- system through a well at constant rate, an oil bank is built up, gas being expelled from the space taken by the oil as shown in Fig. 1. The saturations within each

By S. C. Jones, J. A. Davis

This paper describes displacement mechanisms of micel-lur solution slugs, displuced by a thickened water "mobility buffer", in a glass micromodel and in consolidated Berea sandstone cores. Colored motion pictures, taken through a microscope, recorded interactions on a pore-to-pore bash in the micromodel. Production and differential pressure data from consolidated cores previously waterflooded to residua1 oil saturation revealed displacement mechani.sm.s macroscopically. Oil is displaced miscibly, while interstitia1 water is displaced immiscibly by the injected micellar. solution slug. Part of this water is dispersed and carried near the lead ing edge of the .slug, and some may remain as a residual saturation. Thickened water displaces any residual water. It also emulsifies with the. slug material, thereby displacing it. The displacement sequence forms banks of oil and water with relatively constant fractional flows, saturations, and mobilities. The mobility of oil and water flowing ahead of the .slug defines .slug viscosity required for stable displacement. Equations are presented that describe. several aspecis of the displacement flow behavior: oil breakthrough, oil cut during production of the stabilized bank, breakthrough of banked water and the resulting decreased fractional flow of oil, fluid mobilities, etc. Breakthroughs, calculated from the equations, are compared with experimental data. Slugs of micellar sollution, ranging from 1 to 5 percent PV, produced unusually high oil recoveries (up to 100 percent) in previously waterflooded consolidated cores. Investigation of the effect of mobility buffer size on oil recovery indicated that recovery is reduced in 3-in. X 4-ft cores for mobility buffer sizes .smaller than 50 percent of the pore volume. Introduction A new oil recovery process,* conceived in 1962, utilizes a slug of micellar solution' that has sufficient viscosity to avoid an unfavorable mobility ratio with respect to the oil and water it displaces. These micellar solutions are composed of hydrocarbons, water, and surfactants, and may contain electrolytes and co-surfactants. The micellar solution is displaced by a "mobility buffer". The mobility buffer is a bank of reduced-mobility fluid, preferably an aqueous polymer solution. The entire system is displaced by unthickened water. This paper describes some of the microscopic and macroscopic displacement mechanisms observed during laboratory development of the fluids for a potential field test. Microscopic phenomena were observed photographically in a specially constructed glass flow model and macroscopic data were taken from constant rate, tertiary displacements conducted in Berea sandstone cores. Qualitative Description of Micellar Solution Flooding This section qualitatively describes displacement mechanisms that occur when a slug of micellar solution is injected into a previously waterflooded core. The slug completely removes oil from the portion of core it contacts, and banks oil in the process. The oil saturation in this bank and the flowing WOR reach definite, fixed values that are determined by water and oil mobilities and slug characteristics. This zone of constant WOR and oil saturation will be called "stabilized oil bank" in the following discussion. When this stabilized bank reaches the outlet end of the core, the oil cut jumps from zero to a constant value equal to the fractional flow of oil. Saline water has little or no solubility in the slug. Because of this. the slug displaces in-place brine immiscibly. During slug injection, most of the brine is displaced ahead of the slug, some is carried along in the leading portion of the slug, and part of the brine may become immobile. The immobile, residual brine is later miscibly displaced by thickened water. Displacement of all the oil and only part of the water during slug injection results in a higher fractional flow of oil ahead of the slug than occurs when all the water is displaced after injection of thickened water is begun. The slug spontaneously emulsifies fresh (low-salt content) water. By emulsification, the slug can become infinitely diluted, and hence totally displaced, by thickened water. Fluid saturations in the core are shown schematically in Fig. 1. Parts A and B show conditions before and after the start of the thickened-water injection, respectively. Since thickened water displaces any residual brine by-

Jan 1, 1969

By A. R. Gregory, A. L. Podio, K. E. Gray

Dynamic elastic properties of dry and water-saturated Green River shale samples were computed from compressional- and shear-wave velocity measurements. P- and 5-wave velocity measurements were made in three mutually perpendicular directions with respect to the bedding planes. Measurements were also made in several different directions by varying the angle between the bedding planes and the direction of propagation of the wave for angles of 0, 30, 45, 60 and 909 The oriented samples were subjected to both confining pressure and axial loads, in excess of the confining stress, in the direction of propagation. In general, P- and S-wave velocities increased with increasing stress levels, with a corresponding increase in Young's modulus. Water saturation caused the P-wave velocity to increase and the S-wave velocity to decrease. Elastic moduli decreased upon saturation, except for Pois son's ratio, which increased, indicating some degree of weakening of the material. The samples showed a moderate degree of anisotropy; this was to be expected from the laminated nature and shallow occurrence of Green River shale. INTRODUCTION This paper presents some results of an experimental determination of the elastic coefficients of anisotropic materials (in particular, finely layered rocks and minerals such as Green River shale) from measurements of dilatational and shear-ultrasonic-wave velocities. Ultrasonic techniques have been used extensively in nondestructive testing. Several methods have been proposed by Mcskimmin, 1 and some of these have been used successfully to measure ultrasonic velocities in rocks. Hughes and Cross, 2 Wyllie et al., 3 and Birch 4,5 developed pulse first-arrival techniques for the measurement of dilatational and shear velocities. Williams and Lamb 6 used the method of cancellation of a traveling wave, which was later modified by Myers et al.' and perfected by McSkimmin. 8 Although this method is highly accurate, it has not been used as widely as the pulse-transmission methods recently reported by Jamieson and Hoskins 9 King,10 and Mattaboni and Schreiber. 11 It has been common practice to use some form of crystal transducers, either quartz or ceramic, that has been cut or polarized in different directions in order to generate either compressional or shear waves. However, accurate determination of shear wave velocities has been difficult due to problems that arise in obtaining a pure shear wave from cross-polarized crystals, which usually also generate a small amount of compressional energy. As reported by Gregory, 12 this energy can be seen as a long precursor preceding the sharp break of the shear first arrival. The need for generating pure shear waves led to interest in mode-conversion techniques, which are based upon conversion of the mode of vibration through wave reflection or refraction at a discontinuity. Arenberg13 showed that for certain materials and for certain angles of incidence it is possible to generate pure shear modes by reflection at a boundary. Jamieson and Hoskins 9 used a pyrex glass-air interface for generating pure shear waves, and King 10 used this method successfully for measuring shear-wave velocities in rocks. Gregory l4 arrived at a similar result by refraction of a wave at an aluminum-oil interface. A plane compressional wave, traveling in the oil phase, is incident on the aluminum at an angle larger than the critical angle for compressional waves, and thereby generates a purely transverse, plane-polarized wave in the aluminum. During the last few years methods have been developed that allow the simultaneous determination of shear and compressional velocities in solids.

Jan 1, 1969

By D. O. Wyble

Five water-driver fields, for which pressure predictions had been made within a few years after discovery, were re-studied when four to five years of additional pressure and production data became available. It was fortnd that forecasts of the pressure in water-dri~fe reservoirs may be subject to considerable unccrtainty. Predictions car be erroneous when there is insufficient presswe history to define the cllaracteristics of the aquifer, and when there are errors in the pressure, production or geologic data. INTRODUCTION The efficient operation of an oil or gas reservoir requires that the important characteristics of the reservoir be determined soon after its discovery. It is especially desirable to evaluate the pressure-production relationship quickly, for this relationship provides a basis for deducing the type of production mechanism that is operative, and for making pressure predictions. The pressure predictions, in turn, provide a basis for calculating efficient rates of production, determining whether water or gas should be injected, etc. The importance of pressure prediction is such that forecasts should be made as accurate as possible, and the limits of error well defined. Historically, pressure forecasts for some fields have proved very accurate, while for others the forecasts were so erroneous that they resulted in the adoption of wrong, or at best, unnecessary production practices. Pressure predictions made for water-drive fields are especially subject to error. Accurate forecasting for these fields requires that much be known not only about the reservoir proper but about the adjacent water-filled rock as well. The water-containing rock, or aquifer, may extend for many miles, and often there is only meager data regarding its thickness, permeability and porosity. Much of the information must be inferred from the previous pressure behavior of the reservoir. A brief pressure history may provide little information about the characteristics of the aquifer. Hence, pressure predictions made early in the life of water-drive fields can be incorrect. The work described in this paper was done to determine the degree of reliability which can be expected of pressure predictions made within a few years after water-drive fields are discovered. Five fields, for which the original pressure predictions had been based on less than seven years' history, were re-examined with the aid of four to five years of additional pressure data. Where the predictions were inaccurate, the reasons for failure were ascertained. The studies were accomplished by use of the electric reservoir analyzer.' The analyzer was first adjusted to repeat the original pressure-production match using the same values for the various parameters as were used in the original study. Then, the actual production rates experienced during the time elapsed since the initial study were set into the analyzer, and the pressure behavior indicated by the machine for this period was compared with that which actually OCcurred. Where the pressure profile traced out by the analyzer deviated from the observed field behavior, the analyzer was readjusted to obtain a better over-all match. The old and the new solutions were compared to determine the reasons for any differences. COMPARISON OF PREDICTED RESERVOIR PRESSURE BEHAVIOR WITH ACTUAL PERFORMANCE Field A The original pressure-production match is entirely adequate as evidenced by the pressure match of the field's history subsequent to Jan., 1952 (see Fig. 1)'.