Search Documents

By O. R. Wagner, R. O. Leach

Immiscible displacement tests were performed in a consolidated sandstone core over the interfacial tension range from less thdn 0.01 to 5 dynes/cm to better define how interfacial tension (IFTJ reduction can lead to increased oil recovery. Data obtained were displacement efficiency at breakthbrough vs IFT for both drainage and imbibition conditions. These tests simulate water flooding under oil-wet and water-wet conditions, respectively. Results of the study have shown that displacement efficiency under both oil-wet and water-wet conditions can be markedly improved by a sufficient reduction in IFT. In the particular porous media used and for the low pressure gradients employed, the IFT must be reduced to a value less than about 0.07 dynes/cm to achieve increased recovery at the time of breakthrough of the injected phase. Below 0.07 dynes /cm, further small reductions in IFT result in large increases in displacement efficiency. Observed increases in recovery were obtained at pressure gradients which are well below those which can exist in the interwell area of a reservoir under water flood. The effect of pressure gradient on recovery is discussed. INTRODUCTION A residual oil saturation remains in rock which has been water flooded because, under usual reservoir conditions, the driving force which can be generated is inadequate to expel oil trapped by capillary forces. Since these capillary forces can be reduced by reducing the IFT a frequently studied method of increasing oil recovery has been the use of surfactants to reduce the water-oil interfacial tension. Mungan l observed improved recovery in both water-wet and oil-wet systems at 1.1 dyne /cm IFT, finding that the amount of improvement was greater in oil-wet systems. Moore and Blum,2 working with visual flow cell micro-models, con- cluded that recovery cannot be improved in water-wet systems by injecting surfactant solutions. They calculated that, for a pressure gadient of 1 psi/ft in their model, the IFT must be reduced to 0.03 dynes/cm to release oil trapped under water-wet conditions. Berkeley et al.3 indicated that, at representative field flooding rates, the IFT must be reduced to 0.001 dyne/cm to improve recovery, Thus, there is a wide variation of opinion about the IFT levels needed to improve recovery under field conditions. These previous investigators all have used as their experimental method the addition of surfactants to the injected water. The use of surfactant solutions to reduce IFT creates two experimental problems: (1) the loss of surfactant through adsorption on reservoir rock can obscure the true IFT value which exists at the displacement front, and (2) at very low values of IFT emulsifica-tion of oil and water commonly occurs. It is difficult, therefore, to determine whether the increased recovery is caused by IFT reduction as such, or is instead caused by emulsification. Also the properties of the available surfactants limit the IFT range which can be studied. In the present study the purpose is to better define how interfacial tension reduction can lead to increased oil recovery. A matter of great interest is the amount of recovery improvement potentially achievable in this way. The study was made using very low pressure gradients which was well within the range achievable by water flooding in the interwell region of petroleum reservoirs. A unique experimental approach was chosen to avoid adsorption and emulsification problems, and to allow convenient control of IFT. The fluid phases used were the equilibrium vapor and liquid phases of the methane-n-pentane system. The interfacial tension level was varied by changing the equilibrium pressure of the methane-pentane system over the range from 1,200 psia to near the critical pressure (2,420 psia). All tests were performed at a controlled temperature of 100F. The IFT-vs-pressure relationship for the methane-pentane system was based on the data of Stegemeier and Hough,6 and on new data obtained in this study. With this experimental approach it has been possible to study displacement at lower values of IFT than have been previously investigated.

Jan 1, 1967

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

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

By D. L. Katz, M. J. Rzasa

The derivations of three methods of computing the static pressure gradients in natural gas wells have been presented to show the assumptions made. Charts were developed from which the pressure gradients may be read when the well-head pressure, the well fluid gravity, depth, and the average well temperature are given. A chart for estimating the well fluid gravity from the condensate content and separator gas gravity is included. The effect of the increased average well temperature after flow on the calculation of the static pressure gradient is discussed. Introduction Reservoir pressures have been calculated from well-head pressures for gas wells for many years. As the pressure measurements become more accurate, the need for a reliable calculation of the static pressure gradient often arises. This paper will develop the several methods for computing the pressure gradient in gas wells and make a comparison between them. The method of calculating static pressure gradients in common use is Eq. I (ref. I) or its counterpart which includes a factor for the deviation of the gas from the ideal gas law. P2 - P1 = P1(e0.0000347GX _ I) Eq. [I] in which P1 = pressure at well head, lb. per sq. in. abs. PP = pressure at bottom of well, Ib. per sq. in. abs. G = gas gravity X = depth of well, ft. An alternate method is to compute the average density of the gas in the well and multiply by the well depth in a manner similar to that used for liquid gradients. Derivation of Formulas A static pressure gradient is a special case of the general fluid-flow equation. Consider a pound of fluid flowing in a vertical column from point I to point 2, Fig. I. By an energy balance for fluid flowing from I to 2: U1 + P1V1 + u12/2g + X1 + q = U2 + P2V2+ u22/2g+X2 +W [2] in which U = internal energy, ft-lb. per lb. P = pressure, lb. per sq. ft. V = volume, cu. ft. per lb. u = average linear velocity, ft. per sec. X = height above datum, fl. q = heat absorbed by fluid, ft-Ib. pcr lb. W = work done by system, ft-lb. per lb. An energy balance on the fluid itself gives: U = ?TdS - ?PdV + etc. [3] in which T = temperature, deg- R- S = entropy, ft-lb. per deg. R. per lb.

Jan 1, 1945

By R. S. Sinha

A properly designed pre- and post-grouting program aided with coherent construction sequence and adoption of flexible design approach provide a practical combination to solve excavation problems in a hostile rock environment: intensely fractured and sheared. The Pacheco Pumping Plant demonstrates the utility of an engineered grouting program.

Jan 1, 1984

By S. L. GILLAN

OF the forty or more excursions provided for the delegates to the World Engineering Congress at Tokyo, the trip to the Ashio copper mine stands out as one of the most enjoyable. In every detail looking to our comfort and pleasure the Furukawa company, our hosts, were admirable. This company operates five other copper mines, also gold and coal mines. Ashio is the largest with an annual production of about 14,000 tons of copper. It has been continuously worked since 1620, and doubtless was one of the sources of copper used in the gigantic statues of Buddha of many hundred years earlier date.

Jan 1, 1930

By H. N. Hall

Reservoir and producing characteristics can govern the decision to use either a one-, two- or three-dimensional model for making predictions for gravity-drainage reservoirs. Examples of conditions requiring one-, two- and three-dimensional calculations are given. In 1961 the author presented a method for predicting one-dimensional gravity-drainage performance. This work has been exrendcd to obtain a three-dimensional model, which is described and for which a sample problem is presented. Introduction Many papers on gravity drainage or gas-cap drive'-3 have been presented for conditions where the reservoir is treated as a single areal segment so that, in effect, only one-dimen-sional flow is considered. Other authors'-Q ave given evidence showing that the one-dimensional concept would be unsatisfactory to predict reservoir performance for certain types of reservoirs. Experience has indicated that specific reservoir and producing conditions govern whether one-, two- or three-dimensional concepts must be used as a basis for reliable reservoir performance predictioils. Previously the author presented a paper outlining a method for predicting the performance of gravity-drainage reservoirs.' That paper treated the reservoir as a single areal segment so that it was applicable for only one-dimensional flow. The work has been extended to obtain a three-dimensional model that can be used to predict performance of gravity-drainage reservoirs. Although capable of being three-dimensional, it can be used equally well for one- and two-dimensional predictions. This paper discusses those reservoir and producing conditions governing the choice of model, and also describes a three-dimensional model suitable for predicting reservoir performance in the more complex situations. Discussion The Choice of Model To predict reservoir performance, it is necessary to represent the reservoir by a physical model. A mathematical model, based on the assumed physical model, then is used to predict the movement of oil and gas within the reservoir. The word "model" is used here to connote the combination of both a physical model and a mathematical model based on the physical model. Certain types of reservoir situations may be predicted suitably with a one-dimensional model, whereas two- and three-dimensional models may be required for other situations. This section illustrates reservoir conditions associated with the use of a specific type of model. Many types of reservoir structures are such that gravity drainage possibly could be an important factor in oil production. Production practices, as well as the nature of the structure. can influence the type of model that should be used. obviously, it is impossible to discuss all combinations of reservoir structure and production practices that might be encountered in gravity-drainage reservoirs. The following are typical examples of conditions where a choice must be made among a one-, two- or three-dimensional model to predict reservoir performance satisfactorily. Conditions Requiring a One-Dimensional Model Fig. la shows contours for an asymmetrical anticline in which gravity drainage might be expected to be important. Fig. Ib depicts a cross-section along Line A-A in Fig. la. Existence of a gas cap is a good clue that gravity drainage will be an important factor in producing oil from this type of reservoir. However, it is not essential since a gas cap will form as oil is produced from the reservoir. If the proper combination of reservoir permeability and withdrawal rates is present, the reservoir will behave like a

Jan 1, 1969

By M. A. Jones

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

Jan 1, 1967

By C. K. Eilerts

Finite difference equations were programmed and used to integrate the second-order, second-degree, partial differential equation with variable coefficients that represents the transient linear flow of gas-condensate fluids. Effect was given to the change with pressure of the compressibility factor, the viscosity, and the effective permeability and to change of the absolute permeability with distance. Integrations used as illustrations include recovery of fluid from a reservoir at a constant production rate followed by recovery at a declining rate calculated to maintain a constant pressure at the producing boundary. The time required to attain such a limiting pressure and the fraction of the reserve recovered in that time vary markedly with properties of the fluid represented by the coefficients. Fluid also is returned to the reservoir at a constant rate, until initial formation pressure is attained at the input boundary, and then at a calculated rate that will maintain but not exceed the limiting pressure. The computing programs were used to calculate the results that would be obtained in a series of back-pressure tests made at selected intervals of reservoir depletion. If effect is given to the variations in properties of the fluid that actually occur, then a series of back- pressure curves — one for each stage of reserve depletion — is required to indicate open-flow capacity and related flow characteristics dependably. The isochronal performnce method for determining flow characteristics of a well was simulated by computation. INTRODUCTION The back-pressure test procedurel is based on a derivation of the equation for steady-state radial flow of a gas, the properties of which are of necessity assumed to remain unchanged in applying the test results. The properties of most natural gases being recovered from reservoirs change as the reserve is depleted and pressures decline, and the results of an early back-pressure test may not be dependable for estimating the future delivery capacity of a well. The compressibility factor of a fluid under an initial pressure of 10,000 psia can change 45 per cent and the viscosity can change 70 per cent during the productive life of the reservoir. There are indications that the effective permeability to the flowing fluid can become 50 per cent of the original absolute permeability before enough liquid collects in the structure about a well as pressure declines to permit flow of liquid into the well. The advent of programmed electronic computing made practicable the integration of nonlinear, second-order, partial differentia1 equations2.3 pertaining to flow in reservoirs. Aronofsky and Porter4 represented the compressibility factor and the viscosity by a linear relationship, and integrated the equation for radial flow of gas for pressures up to 1,200 psi. Bruce, Peaceman, Rachford and Rice5 integrated the partial differential equations for both linear and radial unsteady-state flow of ideal gas in porous media. Their published results were a substantial guide in this study of integration of the partial differential equation of linear flow with coefficients of the equation variable. The computing program was developed to treat effective permeability as being both distance-dependent and pressure-dependent. In this study all examples of effective permeability are assumptions designed primarily to aid in developing programs for giving effect to this and other variable coefficients. The accumulation of data for expressing the pressure dependency of the effective permeability is the objective of a concurrent investigation. FUNDAMENTAL DIFFERENTIAL EQUATIONS DIFFERENTIAL EQUATION OF FLOW Linear transient fluid flow is represented by the partial differential equation,

Jan 1, 1965

By D. B. Kibler

THE Florida hard-rock field extends from Suwanee and Columbia Counties in northwest Florida to south of Croom, Florida, in Hernando County. This area is approximately 100 miles long and varies from 2 to 30 miles in width. Florida hard-rock phosphate was discovered at Hale Springs, about one mile northwest of the town of Dunnellon, Marion County, in the year 1888, by Albertus Voight, who lived within a few hundred yards of the spring. Soon after this discovery there was a great influx of land speculators and wildcat operators into the vicinity of Dunnellon. A few years later phosphate was found in Alachua, Levy, Citrus and Hernando Counties, and a great number of small crude plants were built from Newberry, Alachua County, all the way to Croom, Hernando County. Phosphate was found in many places as "outcrops." Looking for outcrops was the first method used by the early prospectors in searching for deposits. Several years passed before any systematic method of prospecting was used. EARLY MINING The early method of mining hard-rock phosphate was with picks and shovels, wheelbarrows, mules and wagons. The pick-and-shovel method was used altogether until about 1900, when some of the larger operators began to use railroad-type steam shovels of from ½ to 1-yd. capacity. Picks and shovels, with negro labor loading the phosphate ore (called "matrix") into wooden side-dump skip-cars, was continued until about 1907. From the discovery of hard-rock phosphate until approximately 1905, overburden was removed with mules and scrapers. At some places, the overburden was loaded into side-dump skip-cars by negro laborers and the cars were then pulled up an incline with a steam hoisting engine. The maximum depth of overburden that could be economically removed by hand labor, mules and scrapers was approximately 15 feet. About 1905 one of the companies erected a producer-gas powerhouse near Dunnellon and generated its own electricity and about the same time the larger companies began using hydraulic methods for the removal of overburden. From the year 1888 through the year 1939, approximately 75 companies were operating at one time or another in the hard-rock territory, on approximately 500 locations. The hard-rock deposits are more irregular and have higher mining costs than the pebble-rock deposits, therefore they have not been competitive in the domestic market, but because of the higher grade of the hard rock it is in demand for the export trade. Discovery of phosphate outside of the United States has made it increasingly difficult in recent years for the hard-rock companies to survive, so today only three companies are operating. These

Jan 1, 1941

By J. H. Henderson, J. Naar

The displacement of a wetting fluid from a porous medium by a non-wetting fluid (drainage) is now reasonably well understood. A complete explanation has yet to be found for the analogous case of a wetting fluid being spontaneously imbibed and the non-wetting phase displaced (imbibition). During the displacement of oil or gas by water in a water-wet sand, the porous medium ordinarily imbibes water. The amount of oil recovered, the cost of recovery and the production history seem then to be controlled mainly by pore geometry. The influence of pore geometry is reflected in drainage and imbibition capillary-pressure curves and relative permeability curves. Relative permeability curves for a particular consolidated sand show that at any given saturation the permeability to oil during imbibition is smaller than during drainage. ' Low imbibition permeabilities suggest that the non-wetting phase, oil or gas, is gradually trapped by the advancing water. This paper describes a mathematical image (model) of consolidated porous rock based on the concept of the trapping of the non-wetting phase during the imbibition process. The following items have been derived from the model. 1. A direct relation between the relative permeability characteristics during imbibition and those observed during drainage. 2. A theoretical limit for the fractional amount of oil or gas recoverable by imbibition. 3. An expression for the resistivity index which can be used in connection with the formula for wetting-phase relative permeability to check the consistency of the model. 4. The Iimits of flow performance for a given rock dictated by complete wetting by either oil or water. 5. The factors controlling oil recovery by imbibition in the presence of free gas. The complexity of a porous medium is such that drastic simplifications must be introduced to obtain a model amenable to mathematical treatment. Many parameters have been introduced by others in Iprogressing" from the parallel-capillary model to the randomly interconnected capillary models independently proposed by Wyllie and Gardner3 and Marshall4. To these a further complication must be added since an imbibition model must trap part of the non-wetting phase during imbibition of the wetting phase. Like so many of the previously introduced complications, this fluid-block was introduced to make the model performance fit the observed imbibition flow behavior. THE IMBIBITION MODEL ASSUMPTIONS The following assumptions are common to many two-phase models previously described. 1. The wetting characteristics of the medium are clearly defined. 2. When two fluids co-exist in the medium, the wetting phase will fill the smaller pores and the non-wetting fluid the bigger ones. 3. The pores have a circular cross section. 4. Apparent pore radii are related to the capillary pressure (drainage) by

By W. A. Klikoff, I. Fatt

The concept of fractional wet wattability is examined. Fractional water wettability of a reservoir rock is defined as the fraction of the internal surface urea that is in contact with water. Capillary pressure and relalive permeability of unconsolidated sand are shown to be functions of fractional wettability. INTRODUCTION The petroleum industry has long recognized that wettability of reservoir rock has an important effect on multiphase flow of oil, water and gas through reservoirs. As early as 1928 the American Petroleum Institute sponsored a study of wettability as part of API Project 27 at the U. of Michigan.' Despite 30 years of research, there is still little exact knowledge of the wettability of reservoir rocks. There are two parts to the wettability problem. After agreeing to a uniform nomenclature in regard to wettability,' the first question to be answered is, "What is the in situ wettability of a given reservoir rock?" If this can he answered the next question is, "What part does wettability play in determining the characteristics of multiphase fluid flow through the rock?" This paper represents an oblique attack on the problem of wettability. No attempt is made here to answer the basic question of wettability in situ. instead the consequences of the concept of fractional wettability are examined. Multiphase flow in sandpacks is shown to he highly influenced by fractional wettability. Jennings' has given 3 definition of wettability and the other terms used in discussing wettability. These terms must be applied to the physical situation existing In reservoir rock. A survey of the pertinent literature from 1928 to 1956 indicates that the concept of a contact angle was applied to reservoir rock in the same way it would be applied to a flat, homogeneous surface. Attempts were made to state quantitatively the wettability of a reservoir rock in terms of a contact angle which was presumably constant at all points on the very rough and heterogeneous interior surface of a porous rock. Calhoun, et al,3-6 prepared synthetic consolidated and un-consolidated porous media in which they claimed there was a known uniform contact angle. They then showed the effects on the capillary pressure and relative permeability characteristics of varying this angle. The API Project 47 at the U. of Texas' and others' have made extensive studies of an indirect approach to the contact angle through the use of heat of wetting data. Even if successful. however, this approach also states the wettability of porous rock in terms of a contact angle which is uniform over the entirc surface. If the angle varies from one part to another on the internal surface, there is no way of determining from the measurements the area distribution of contact angles. In 1956 Brown and Fatt5 suggested that the concept of a contact angle, as applied to reservoir rock, be abandoned. This suggestion was made because it is known that the internal surface of most reservoir rocks is composed of many different minerals, cach with a different surface chemistry and a different capacity to adsorh surface active materials from reservoir fluids. Furthermore, the operation of a contact angle in determining the form of a fluid-fluid interface is difficult to picture in the very complex geometry of a pore. Brown and Fatt proposed that the wettability of reservoir rock be stated in terms of the fractional internal surface area that is in contact with water or oil. All surfaces on which there is water are called water-wet; surfaces on which there is oil are called oil-wet. The fractional water wettability is then stated as a number which represents the fraction of the internal surface that is in contact with water. A symmetrical statement can be made for the fractional oil wettability. The concept of a fractional wettability as previously stated has in its favor the recognition of the heterogeneous mineral composition of most reservoir rocks. Another point in its favor is that fractional wettability can be measured quantitatively with relative ease. Hol-brook and Bernard10 se a simple dye adsorption test to obtain fractional wettability of reservoir rocks. Amott11 uses a combination of imbibition and displacement to arrive at a wettability index of reservoir rocks which seems related to fractional wettability in the range 0.25 to 0.75 fractional water wettability. Jennings" has shown the changes in relative permeability that take place when a porous material is changed from unity to zero fractional water wettability. He also shows that reservoir rocks in their natural state, but at room temperature and atmospheric pressure, have relative permeability characteristics which would indi-

By P. B. Crawford, W. L. Huskey

A study has been made to determine the effects of random vertical fractures existing in a reservoir matrix on the effective matrix properties and producing characteristics of a well. The study indicated that short vertical fractures may distort the isopotentials by 1 percent in the vicinity of the fracture. A correlation was developed which showed the effect of fracture density on the producing capacity of a well. The correlation indicated that fracture shape had little effect on producing capacity, but that total fracture length was closely correlated with the producing capacity. Effective values of reservoir permeability were found to correlate closely with fracture density when measured parallel to the direction of flow. INTRODUCTION In reviewing the petroleum literature, it is found that in most analytical solutions of hydrocarbon reservoir behavior the assumption is made that the rock is uniform throughout although core analysis and well logging surveys indicate that substantially all reservoirs are heterogeneous and some are fractured. Only in recent years have investigators studied the effects of reservoir hetrogeneity. Russell1 discussed the effect of fractures in reservoirs. The literature indicates that fractures greatly affect the exploitation of man reservoirs. This is described by Elkins and Skov2 and Little-field et al. 3 Several types of reservoir heterogeneity have been studied by unsteady-state methods. In 1959, Landrum et al. 4 showed that transient phenomena in hydrocarbon reservoirs could be studied by thermal models. pickering5 studied the flow of heat in stratified linear reservoirs consisting of plates of different metals separated by sheets of rubber which were suddenly exposed to a low-temperature source at one end. Cotman6 studied both laterally and vertically heterogeneous thermal model reservoirs with special reference to cross flow. Another type of reservoir heterogeneity involves the presence of highly permeable, circular lenses around or near the wellbore which was studied by Miesch7 using a thermal and a steady-state model. Many studies have been made of the effects of induced hydraulic fractures on well performance. 8-10 Work regarding vertical fractures that exist initially in the reservoir was reported by Givens, et al.11 Since there appears to be little published work on the importance and influence on reservoir performance of fractures that exist out in the rock matrix, this study was conducted to determine the effects of vertical fractures on pressure distribution, producing capacity and the effective permeability of the reservoir rock. The analogy between the flow of electrical current through a conducting body and the flow of fluid in a porous medium provides a convenient method of studying various aspects of petroleum reservoirs. For electrical radial steady-state flow,12 For steady-state fluid flow, With the above analogy, it is only necessary to establish the ratio of the values of the electrical units in terms of the corresponding flow units to make the analogy quantitative. The ratios are shown below.

By R. C. Mahon

During the early history of the Lake Superior district, tamarack timber was used almost exclusively in the wines. It had the strength and lasting qualities necessary for a satisfactory mining timber. But about fifteen years ago the tamarack trees were attacked by a blight which killed nearly all of them so today no tamarack timber is available. Hemlock, spruce, and poplar have neither the necessary strength nor the lasting qualities for mining timber, but a large supply of maple and birch is available. These have ample strength but, when placed underground, fail by decay in about three years. This means that, instead of replacing the timber sets once in fifteen years, it must be done five or more times in the same space of time which involves not only incraesed timber costs but more important, an increased labor cost.

Jan 1, 1949

By D. Havlena, A. S. Odeh

The use of the straight-line method of solving the material balance equation is illustrated by means of six field cases. Also, the application of statistical criteria to arrive at the most probable answer is shown. The theory underlying the straight-line method of solution and the applicability of the statistical criteria was presented in a previous paper.' The field cases include saturated and undersaturated oil reservoirs with and without water drive. The aquifers discussed are: limited radial, infinite radial, very small aquifer and infinite linear. The field cases also include a gas reservoir producing under water drive. INTRODUCTION In a previous paper,' the authors presented the theory underlying the solution of the material balance as an equation of a straight line. The appropriate equations for various material balance cases as well as the methods of analysis and interpretation with comments and discussion were also included. To illustrate the various theoretical cases treated previously, five field cases are analyzed in this paper by employing the straight-line method of solving the material balance equation (MBE) and one example previously published is referred to. The use of statistical criteria to arrive at the most probable answer is also shown. All the field examples presented, except Case 2, are excerpts from complete reservoir studies. To illustrate the method, only sections specifically dealing with the material balance principles are included. Additional geologic information and basic data are reported to better acquire an understanding of the cases and thus to better follow the reasoning that suggested the successful application of the straight-line method of solving the MBE. The six cases are: (1) saturated reservoir, small gas cap, limited aquifer; (2) saturated reservoir, very small gas cap, infinite aquifer; (3) undersaturated-saturated reservoir, very small aquifer; (4) highly undersaturated reservoir, no water drive; (5) high undersaturated one-well reservoir, limited aquifer; and (6) gas reservoir, infinite linear aquifer. THE D4 SAND, GUICO FIELD, VENEZUELA The D, sand, which was discovered in 1943, is presently in a depleted state. Since its discovery it has produced under water drive, gas-cap-gas expansion, and solution gas drive. In Nov., 1947, water injection was initiated to arrest further pressure decline. When discovered, the D, sand was a saturated reservoir with a gas cap /oil zone volume ratio m estimated volu-metrically at 0.0731, an average permeability of 500 md, a porosity value of 25 per cent, and an oil viscosity at reservoir conditions of 0.3 cp. The volumetrically determined stock-tank oil initially in place was 23.1 million bbl. The volumetrically weighted physical data and production data available until Nov., 1953 are reported in Table 1. In Ref. 1, the effects on the straight-line plot of various values of ro/r for a constant ?to, or of various dimen-sionless times for a constant r8/r10 were theorized and were illustrated in Fig. 3A of that reference. In this field case, the previously theoretically predicted effects are established. Thus, the MBE calculations using Eq. 3c of Ref. 1 were performed for various r6/r10 and dimension-less time values. Eq. 3c of Ref. 1 is:

Jan 1, 1965

By M. R. Tek, F. F. Craig, J. O. Wilkes, C. S. Goddin

The waterflood performance of a water-wet, stratified system with crossflow is computed by a finite difference procedure. The effects of five dimensionless parameters on tile oil displacement efficiency, water saturation con-tour.7 and crossflow rates are evaluated in the absence of gravity forces. Crossflow due to viscous and capillary forces is shown to exert a significant effect on oil recovery in a field-scale model of a two-layered. water-wet sandstone reservoir. The crossflow is at a maximum in the vicinity of the front advancing in the more permeable layer. Under favorable mobility ratio conditions, the comparted oil recovery with crossflow always is interrnerliate between that predicted for a uniform reservoir and that for a layered reservoir with no crossflow. INTRODUCTION The important erects of reservoir heterogeneity on waterflood performance are commanding increased attention in the technical literature. Much of this attention is centered on two categories of layered reservoirs: those in which layers are non-communicating and those in which crossflow of fluids occurs between the layers. In the first category, the reservoir is assumed to consist of discrete layers, each uniform within itself and differing from the others only in such properties as thickness, porosity and absolute permeability. The performance within each layer is calculated by one-dimensional flow theory, and the performance of the total reservoir is obtained by summing individual layer performances. Capillary and gravity effects usually are not considered. Representative publications dealing with thi5 type of reservoir are those of Stiles,' Dykstra and Parsons,' Hiatt,3 Warren and Cosgrove' and Higgins and Leighton." Prediction of performance for reservoirs in the second category is considerably more difficult since viscous, capillary and gravitational forces all play important roles in causing crossflow between layers. A number of authors have investigated the simpler problem of two-dimensional displacement flow in a stratified system with a mobility ratio of unity and negligible capillary and gravity effects.'; Others have considered two-dimensional, non-steady-state flow of a single, slightly compressible fluid in a stratified reservoir. A limited number of laboratory oil displacement tests in layered models with crossflow have been reported. Miscible floods (with resultant zero capillary forces) in layered five-spot models were conducted by Dyes and Braun," who studied the effect of mobility ratio with zero gravity forces, and by Craig et al. 12 who studied the effect of gravity forces at constant mobility ratio. Waterfloods in layered five-spot models (with cross-tlow due to capillary, viscous and gravity forces) were conducted by Gaucher and Lindley,"' who showed the effect of gravity forces in causing underrunning of the injected water and by Carpenter, Bail and Bobek, 14 who demonstrated the reliability of Rapoport's" dimension-less parameters for scaling layered systems. Waterfloods in rectangular layered models were conducted by Richardson and Perkins."' who investigated the effect of velocity at constant mobility ratio and with zero gravity forces, and by Hutchinson," who studied the effects of varying mobility, layer permeability and layer thickness ratios. The differential equations which rigorously describe waterflooding in a heterogenous porous medium are non-linear and do not facilitate analytical solution. By using finite difference approximations it is possible to obtain a solution to any desired degree of accuracy. Such a solution, using an alternating direction implicit procedure (ADJP), is described by Douglas, Peaceman and Rachford.18 In the present study, a computer program using ADJP explores systematically the effects of important parameters on waterflood performance of a two-dimensional, two-layered, field-scale model of a water-wet sandstone system. Particular attention is given to evaluation of the water saturation contours and crossflow rates at the interface between layers to gain improved understanding of the crossflow mechanism. PROCEDURE BASIC: FLOW EQUATIONS The basic flow equations for two-dimensional, two-phase, immiscible, incon~pressible flow in a porous medium are:

Jan 1, 1967

By J. E. Johnson

Discussion of the paper of Mark A. Amnion, presented at the Cleveland meeting October, 1913, and printed in Bulletin No. 70, October, 1913, pp. 1167 to 1 79. J. E. JOHNSON, JR., Ashland, Wis.:-Mr. Amnion's paper is of great value and it is not with any desire to reflect. on it that I wish to make some comments on the use of the scleroscope. We obtained one on trial some months ago and made very extensive tests with it, but without getting very satisfactory results. The readings on good car-wheels and bad ones were just about the same, but the factor in its use most likely to produce erroneous results is a variation in the thickness of the specimen being tested. In order to find out something on this subject we cut a thin disk, about -1,: in. thick, off a brass rod; then one 0.25 in. thick, and another about 1 in. thick, put the surfaces of all in about the same condition and tested them with the scleroscope. The readings varied widely; the resilience of the anvil evidently entered to a very high degree on thin pieces. By putting a soft cushion between such a thin specimen and the anvil, and forcing it up against the tube, the softness of the cushion cut off the resilience of the anvil and the results were more like those of the thick specimens. After a thickness of 1 in. was passed, the specimen seemed to have mass enough so that the resilience of the anvil no longer played an important part, and above that thickness results on different materials would be comparable, even though the thicknesses differed. But in order to get comparable results on thin specimens the thickness of those being tested should be the same. I should like to ask how the variations introduced in this way were overcome, so as to get strictly comparable results. Mr. AMMON :-.I have not had much experience in that, but the scleroscope is simply an indication of the surface harden-

Dec 1, 1912

By Asheim, Harald A.

Parameter estimates affect both the optimum field development decisions and the economic result. A general methodology for assessment of the effects of parameter uncertainty in offshore development decision situations is proposed, and used to analyze aspects of the initial development situation for the Statfjord field.

Jan 1, 1982

By R. O. Leach, O. W. Wagner

Because of unfavorable wetting conditions much residual oil is left when a porous material is Pushed by water. Methods suggested to change reservoir wetting to improve oil displncernrnt efficiency are generally expensitlr. The present 1aborator.y study was undertaken to gain an under.standinx of the factors which determine reservoir wettability, arid to find out if oil displacement efficiency might be improved by a wettahility change accomplished at low cost in on oil reservoir. Contact angle measurements were made on mineral surfaces using sevc.r~zl sets of reservoir oil and water samp1es. Results of the contact angle studies suggest that reservoir wetta-hility may he primarily determined by natural surface-active substances present in the reservoir fluids. The effect of changing sa1inity and pH of the water phase was studied. The re.suits suggest that gross changes in preferential wettability might be acc~o~npli.shed by injection of water containing simple chernicnls to alter pH or salinity in the reservoir. Such treatment could he much less expensive than injection of commercial surface-active agents. Waterflood tests have also been made using synthetic cores and oil and water having wening characteristics similar to those of reservoir fluids. Cores initially oil-wet were flooded in such a way that they were made prefermtial1y water-wct by the advancing flood water. This reversal in preferential wettability achieved greater oil displacement efficiency than when either oil-wet or water-wet conditions were maintained throughout the flood. For the systems studied, the higher the oil viscosity the greater the percentage improvement obtained over conventional waterflood recovery. This suggests that a flooding process making use of wettability-reversal may extend the oil viscosity range over which water flooding is attractive. Because a precise adjustment of reservoir wettability does not seem to be required, and because altering the pH or salinity in some reservoirs may be inexpensive, it appears that a waterflooding process employing wet-[ability-reversal could find .succesful field application. I NTRODUCTION The efficiency with which water will displace oil from a porous material is related to the nature of the capillary forces present. These in turn are controlled by the preferential wetting of the solid by the two fluids. Because of unfavorable wetting conditions, 30 per cent or more of the original oil in place may remain unrecovered in that portion of a reservoir flushed by water. This paper is concerned with the possibility of improving waterflood oil displacement efficiency by alterations in the wettability of the porous material. A laboratory study was made to gain a better understanding of the factors which control reservoir wettability, and to determine if the oil displacement efficiency could be improved by some inexpensive means of manipulating wettability of the porous medium. Contact angle measurements were made with several natural and synthetic oil, water and solid systems (1) to obtain a better understanding of how to duplicate reservoir wettability in the laboratory, and (2) to discover possible means for changing preferential wettability of natural reservoir systems. Flooding tests were also made in synthetic systems to determine if oil displacement efficiency could be improved by those wettability manipulations suggested by the contact angle measurements. Based on these studies a possible method for improving waterflood oil displacement efficiency is presented. This method involves causing an originally oil-wet porous material to become preferentially water-wet during the course of a water flood. The purpose of this paper is to present results of the laboratory studies. THEORY Rock surfaces in some oil reservoirs are believed to be covered with a firmly attached bituminous or other organic coating. Such surfaces would be preferentially oil-wet in the presence of both oil and water, regardless of composition of reservoir fluids. Other reservoirs are believed to contain rock surfaces not permanently coated with such materials, and which would be preferentially wet by water in the presence of water and oil free from surface-active substances. However, when the reservoir fluids do contain certain natural surface-active materials in sufficient quantity, rock surfaces acquire a degree of preferential oil wettability caused by adsorption of these natural surface-active materials on the solid. The equilibrium amount of these materials adsorbed per unit surface area is believed to depend upon their concentration in the bulk liquid phases.

By R. D. West

Miscible displacenzent recovers all oil in the area contacted by the injected .fluid, whereas water or immiscible gas drives usually leave substantial amounts of oil as residual. However, the Door mobility ratios associated with a gas-driven miscible displacement cause the sweep pattern efficiency to be much lower than that obtained with water flooding. One way in which the sweep eficiency in a miscible displacement process can be increased is by decreasing the mobility behind the flooding front. This can be achieved by injecting water along with the gas which drives the miscible slug. This water reduces the relative permeability to gas in this area and thus lowers the total mobility. The main operating conditions for the simultaneous injection -vrocess are that a zone of gas exists between the miscible slug and the leading edge of the water and that a su,@cient amount of gas be injected with the water to form the pas volume which is being left in the water zone. Laboratory model studies have shown that the ultimate sweep pattern efficiency can be as high as 90 per cent for a five-spot flooding system. If gar alone is used as the driving medium an ultimate sweep-out efficiency of about 60 per cent would be obtained in the same system. I INTRODUCTION The miscible displacement processes are a step towards total oil recovery. Conventional gas or water drives usually leave 25 to 5.0 per cent of the oil as residual in the swept portion of the reservoir. This residual can be eliminated if the oil is driven by a fluid with which it is miscible. At some reservoir conditions natural gas will become miscible with the oil. This is the "high pressure gas process".' More often, the oil does not contain enough light hydrocarbons to cause the gas to become miscible with the oil at reasonable pressures. In these cases a small band of fluid which is miscible both with the oil and gas must be kept between them2. Less than 2 per cent of the reservoir volume of the slug material is needed to keep the displacement miscible. Both processes work in the same manner, recovering all of the oil in the portion of the reservoir contacted by the injected fluids. The only difference is the manner in which the miscibility between the oil and the injected gas is obtained. Previous publications have contained detailed descriptions of these processes.1,2,3,4 However, total displacement of the oil in the swept region does not guarantee an efficient recovery process. The amount of oil to be recovered is also determined by the fraction of the reservoir contacted by the flood. This fraction is largely determined by the mobilities of the fluids. (The fluid mobility is the permeability of the rock to that fluid divided by the fluid's viscosity, k/p). This. dependence of the fraction swept on the mobility ratio has been shown in previous studies. Fig. 1 shows the ultimate fraction swept in a five-spot system as a function of the mobility ratio. The small drawings show the location of the areas left unswept for two different mobility ratios. The ultimate fraction of the reservoir swept is here considered to be attained when the producing stream contains less than 5 per cent oil at reservoir conditions. THE GAS-DRIVEN MISCIBLE DISPLACEMENT Since there is no oil left in the swept region after miscible displacement the mobility in this region is very high. It is often 50 times the mobility in the unswept regions. This means that the fraction of the reservoir contacted by the injected fluid will be less for a gas-driven miscible displacement than for a conventional water or gas drive. For a five-spot injection system, water would contact the entire reservoir volume, and the low pressure gas would contact about 90 per cent of this volume, while a gas-driven miscible displacement would only contact about 65 per cent of the reservoir. This poor sweep efficiency often offsets the benefits obtained through miscible displacement. Fig. 2 shows what the recovery curves for the three processes might look like for a five-spot system. The curves show the fraction of the in-place oil recovered as a function of reser-

By G. W. Tracy, R. D. Carter

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