Search Documents

By Fred H. Poettmann

The vaporization characteristics of carbon dioxide in a League City natural gas - Billings crude oil system were studied at three temperatures, 38°. 120°, and 202°F and for pressures ranging from 600 to 8,500 psi. Variation of carbon dioxide concentration up to 12 mole per cent in the composite showed no effect on the equilibrium vaporization ratios (K values) of the hydrocarbon constituents or on the K value of carbon dioxide itself. It was shown that carbon dioxide is more soluble in crudes than in distillates which is contrary to the behavior of methane. A working chart of carbon dioxide K values is presented. INTRODUCTION The study of the equilibrium vaporization ratios of mixtures of paraffin hydrocarbons has been rather thorough.2,6,7,8,9 In the past few years considerable attention has been paid to the vaporization characteristics of the so-called noncondensable gases such as nitrogen, carbon dioxide, and hydrogen sulfide in mixtures of hydrocarbons. since they usually occur to some extent in most crude oils and natural gases.1,3,4,5 Knowledge of this behavior is useful to both the production and refining phases of the petroleum industry. This paper reports the equilibrium vaporization ratios (K's) of carbon dioxide in a mixture of League City natural gas and Billings crude oil, and compares them to those obtained in a natural gas-distillate system. The equilibrium vaporization ratios for the hydrocarbon components in this system had previously been studied by Roland.' In addition to the determination of the K values for carbon dioxide, the K values for methane and ethane were also determined in order to observe what effect, if any, the presence of carbon dioxide had on these K values. The concentration of carbon dioxide was also varied in order to observe the effect of this variable on the carbon dioxide K values. EXPERIMENTAL PROCEDURE The apparatus used in this study cotlsisted of a stainless steel equilibrium cell of about 2 liters capacity. The cell was mounted on trunions permitting rocking in a thermostatically controlled oil bath. Two high pressure valves fitted with steel tubing were mounted on the top of the cell. one was used for sampling the equilibrium gas phase and the other for sampling the equilibrium liquid phase by means of an induction tube within the cell. Stainless steel tubing from the bottom of the cell led to a mercury reservoir and manifold which was connected to a free-piston type pressure gauge manufac- lured by the American Instrunlent Ctr. and to a volumetric. putrip. The temperature of the oil bath was measured by means of a ralibrated mercury-in-glass thermometer. The recorded temperatures are believed to be accurate to ±0.5 °F. The pressures are correct to 22 psi. The crude oil used in this study was stock tank oil obtained from the Wilcox formation in the Billings Field, Noble County. Okla. The natural gas was obtained from the League City Field. Galveston County, Tex. The oil was treated with anhydrous calcium sulfate in order to remove the last traces of water. To insure a supply of constant composition gas at room temperature the cylinders of League City gas were cooled to about 30°F, inverted, and the condensed liquid was allowed to drain from the cylinders. The analysis of the gas and crude are tabulated in Table I. The carbon dioxide came from Pure Carbonic, Inc., and was .stated to have a purity of 99.5 per cent or better. The procedure used to obtain samples of the equilibrium liquid and vapor was similar to that employed by others making use of the rocking type equilibrium cell.6,7,8 The equilibrium cell was evacuated and calculated quantities of carbon dioxide, natural gas, and crude oil were charged to the cell to the desired pressure. In charging the equilibrium cell an attempt was made to maintain the ratio of the natural gas to crude oil as close as possible to that employed by Roland. After the cell was charged, samples of

Jan 1, 1951

By Fred H. Poettmann

The vaporization characteristics of carbon dioxide in a League City natural gas - Billings crude oil system were studied at three temperatures, 38°. 120°, and 202°F and for pressures ranging from 600 to 8,500 psi. Variation of carbon dioxide concentration up to 12 mole per cent in the composite showed no effect on the equilibrium vaporization ratios (K values) of the hydrocarbon constituents or on the K value of carbon dioxide itself. It was shown that carbon dioxide is more soluble in crudes than in distillates which is contrary to the behavior of methane. A working chart of carbon dioxide K values is presented. INTRODUCTION The study of the equilibrium vaporization ratios of mixtures of paraffin hydrocarbons has been rather thorough.2,6,7,8,9 In the past few years considerable attention has been paid to the vaporization characteristics of the so-called noncondensable gases such as nitrogen, carbon dioxide, and hydrogen sulfide in mixtures of hydrocarbons. since they usually occur to some extent in most crude oils and natural gases.1,3,4,5 Knowledge of this behavior is useful to both the production and refining phases of the petroleum industry. This paper reports the equilibrium vaporization ratios (K's) of carbon dioxide in a mixture of League City natural gas and Billings crude oil, and compares them to those obtained in a natural gas-distillate system. The equilibrium vaporization ratios for the hydrocarbon components in this system had previously been studied by Roland.' In addition to the determination of the K values for carbon dioxide, the K values for methane and ethane were also determined in order to observe what effect, if any, the presence of carbon dioxide had on these K values. The concentration of carbon dioxide was also varied in order to observe the effect of this variable on the carbon dioxide K values. EXPERIMENTAL PROCEDURE The apparatus used in this study cotlsisted of a stainless steel equilibrium cell of about 2 liters capacity. The cell was mounted on trunions permitting rocking in a thermostatically controlled oil bath. Two high pressure valves fitted with steel tubing were mounted on the top of the cell. one was used for sampling the equilibrium gas phase and the other for sampling the equilibrium liquid phase by means of an induction tube within the cell. Stainless steel tubing from the bottom of the cell led to a mercury reservoir and manifold which was connected to a free-piston type pressure gauge manufac- lured by the American Instrunlent Ctr. and to a volumetric. putrip. The temperature of the oil bath was measured by means of a ralibrated mercury-in-glass thermometer. The recorded temperatures are believed to be accurate to ±0.5 °F. The pressures are correct to 22 psi. The crude oil used in this study was stock tank oil obtained from the Wilcox formation in the Billings Field, Noble County. Okla. The natural gas was obtained from the League City Field. Galveston County, Tex. The oil was treated with anhydrous calcium sulfate in order to remove the last traces of water. To insure a supply of constant composition gas at room temperature the cylinders of League City gas were cooled to about 30°F, inverted, and the condensed liquid was allowed to drain from the cylinders. The analysis of the gas and crude are tabulated in Table I. The carbon dioxide came from Pure Carbonic, Inc., and was .stated to have a purity of 99.5 per cent or better. The procedure used to obtain samples of the equilibrium liquid and vapor was similar to that employed by others making use of the rocking type equilibrium cell.6,7,8 The equilibrium cell was evacuated and calculated quantities of carbon dioxide, natural gas, and crude oil were charged to the cell to the desired pressure. In charging the equilibrium cell an attempt was made to maintain the ratio of the natural gas to crude oil as close as possible to that employed by Roland. After the cell was charged, samples of

Jan 1, 1951

By Joseph Fry, Ralph H. Espach

In the spring of 1943, when it was evident that the Tensleep bandstone in the Elk Basin Field, Wyoming and Montana, held a large reserve of petroleum, Bureau of Mines engineers obtained samples of oil from the bottom of nine wells and analyzed them for such physical characteristics as the volumes. of gas in solution. saturation pressures or bubble points, shrinkage in volume caused by the release of gas from solution, expansion of the oil with decrease in pressure, and other related properties. The composition of the gas in solution in the oil was studied. The pressures and temperatures existing in the reservoir and the productivity characteristics of the oil wells were determined. The data obtained indicate that the oil in the Tensleep Reservoir of the Elk Basin Field has unusually varying physiral characteristics, such as a saturation pressure of 1,250 psia and 490 cu ft of gas in solultion in a barrel of oil at the crest of the structure and a saturation pressure of 530 psia and 134 cu ft of gas in solution in a barrel of oil low on the flanks. The hydrogen sulfide content of the gas in solution in the oil varies from 18 per cent for oil on the crest to 5 per cent for oil low on the flanks of the structure. Of even greater significance is the fact that these and other variable characteristics of the reservoir oil are related to the position of the oil in the structure. Many geologists and petroleum engineers have considered that all the oil in a petroleum reservoir has rather uniform physical characteristics and that equilibrium conditions prevailed in all underground accumulations of oil and gas; that such is not always so is borne out by the results of the study by the writers. INTRODUCTION The Rocky Mountain region is one in which may be found striking examples of the unusual in oil and gas accumulations, as is evident from the following: The high helium content (7.6 per cent) of the gas in the Ouray-Leadville limestone sequence in the Rattlesnake Field, New Mexico, and gases of similar helium content in other fields; 50" to 55' API gravity distillate in solution in carbon dioxide gas and recoverable through retrograde condensation, in the North McCallum Field, Colorado; the occurrence of gas, oil, or both in closely related structures contrary to the usual concepts of gravimetric segregation: the accumulation of gas and/or oil in structures closely related to other structures that apparently are more favorable but do not contain oil or gas accumulations; the high hydrogen sulfide content (as high as 42 per cent) of the gas associated with oil in some fields in the Big Horn Basin, Wyoming; and the wide range of fluid chararteristics found in the Elk Basin reservoir. Elk Basin, an interesting old oil field that has been producing oil from the Frontier formation since 1915, is situated in a highly eroded basin resulting from the erosion of the crest of an anticline and some of the underlying softer shales. The field came back into national prominence during 1943 when it became known that it was the largest single reserve of new oil discovered in the United States that year. The Tensleep sandstone was found to contain oil in November. 1942, when a well drilled to a depth of 4,538 ft (44 ft into the Tensleep sandstone) flowed oil at the rate of 2,500 B/D. By the end of 1949, 137 oil-producing wells and five dry holes had been drilled, and approximately 32 million bbl of oil had been produced. Approximately 6,000 acres may be considered productive of oil in the Tensleep Reservoir, and estimates of the oil that will be produced average 200 million bbl. The Tensleep Reservoir has further interest because it ha-greater closure than any oil field in the Rocky Mountain region; the closure of the Elk Basin anticline is variously estimated at 5.000 to 10,000 ft. of which the top 2.00 ft of the structure contained oil. SUBSURFACE OIL SAMPLING Fig. 1 is a structural map of the Elk Basin Tensleep Reservoir, on which the nine wells used in this study and the numbers correvponding to the well designations hereafter referred to are shown. Wells 1. 2, 3, 4, and 8 were tested and sampled during October and November. 1943. and Wells 5, 6. 7, and 9 during June and July, 1944. An electromagnetic type sampler developed by the Bureau of Mines and described by Grandone and Cook' was used in obtaining the subsurface oil samples. As the wells were tubed nearly to bottom, the sampler was run as far as possible in the tubing hut never below the top perforations. The following procedure was used in testing and sampling the wells: A well was shut in for at least three days, after

Jan 1, 1951

By Joseph Fry, Ralph H. Espach

In the spring of 1943, when it was evident that the Tensleep bandstone in the Elk Basin Field, Wyoming and Montana, held a large reserve of petroleum, Bureau of Mines engineers obtained samples of oil from the bottom of nine wells and analyzed them for such physical characteristics as the volumes. of gas in solution. saturation pressures or bubble points, shrinkage in volume caused by the release of gas from solution, expansion of the oil with decrease in pressure, and other related properties. The composition of the gas in solution in the oil was studied. The pressures and temperatures existing in the reservoir and the productivity characteristics of the oil wells were determined. The data obtained indicate that the oil in the Tensleep Reservoir of the Elk Basin Field has unusually varying physiral characteristics, such as a saturation pressure of 1,250 psia and 490 cu ft of gas in solultion in a barrel of oil at the crest of the structure and a saturation pressure of 530 psia and 134 cu ft of gas in solution in a barrel of oil low on the flanks. The hydrogen sulfide content of the gas in solution in the oil varies from 18 per cent for oil on the crest to 5 per cent for oil low on the flanks of the structure. Of even greater significance is the fact that these and other variable characteristics of the reservoir oil are related to the position of the oil in the structure. Many geologists and petroleum engineers have considered that all the oil in a petroleum reservoir has rather uniform physical characteristics and that equilibrium conditions prevailed in all underground accumulations of oil and gas; that such is not always so is borne out by the results of the study by the writers. INTRODUCTION The Rocky Mountain region is one in which may be found striking examples of the unusual in oil and gas accumulations, as is evident from the following: The high helium content (7.6 per cent) of the gas in the Ouray-Leadville limestone sequence in the Rattlesnake Field, New Mexico, and gases of similar helium content in other fields; 50" to 55' API gravity distillate in solution in carbon dioxide gas and recoverable through retrograde condensation, in the North McCallum Field, Colorado; the occurrence of gas, oil, or both in closely related structures contrary to the usual concepts of gravimetric segregation: the accumulation of gas and/or oil in structures closely related to other structures that apparently are more favorable but do not contain oil or gas accumulations; the high hydrogen sulfide content (as high as 42 per cent) of the gas associated with oil in some fields in the Big Horn Basin, Wyoming; and the wide range of fluid chararteristics found in the Elk Basin reservoir. Elk Basin, an interesting old oil field that has been producing oil from the Frontier formation since 1915, is situated in a highly eroded basin resulting from the erosion of the crest of an anticline and some of the underlying softer shales. The field came back into national prominence during 1943 when it became known that it was the largest single reserve of new oil discovered in the United States that year. The Tensleep sandstone was found to contain oil in November. 1942, when a well drilled to a depth of 4,538 ft (44 ft into the Tensleep sandstone) flowed oil at the rate of 2,500 B/D. By the end of 1949, 137 oil-producing wells and five dry holes had been drilled, and approximately 32 million bbl of oil had been produced. Approximately 6,000 acres may be considered productive of oil in the Tensleep Reservoir, and estimates of the oil that will be produced average 200 million bbl. The Tensleep Reservoir has further interest because it ha-greater closure than any oil field in the Rocky Mountain region; the closure of the Elk Basin anticline is variously estimated at 5.000 to 10,000 ft. of which the top 2.00 ft of the structure contained oil. SUBSURFACE OIL SAMPLING Fig. 1 is a structural map of the Elk Basin Tensleep Reservoir, on which the nine wells used in this study and the numbers correvponding to the well designations hereafter referred to are shown. Wells 1. 2, 3, 4, and 8 were tested and sampled during October and November. 1943. and Wells 5, 6. 7, and 9 during June and July, 1944. An electromagnetic type sampler developed by the Bureau of Mines and described by Grandone and Cook' was used in obtaining the subsurface oil samples. As the wells were tubed nearly to bottom, the sampler was run as far as possible in the tubing hut never below the top perforations. The following procedure was used in testing and sampling the wells: A well was shut in for at least three days, after

Jan 1, 1951

By C. Özgen Karacan

Degasification and methane control in underground coal mining is an important area augmenting conventional ventilation in order to prevent possible fires and explosions due to excessive methane emissions. Coal bed reservoir engineering has been applied for many years to coalbed methane production. In recent years, modeling the flow and draining methane, not only from coal bed itself but also from the overlying strata, have been recognized as important areas of expertise for supplementing underground ventilation. This paper demonstrates the applications of coal bed reservoir engineering and modeling techniques for optimizing and controlling methane emissions in longwall and continuous miner sections. Practical applications and some considerations are presented for reservoir modeling, borehole performances, and face emission predictions for various coal bed and borehole configurations. A section on field reservoir studies and well bore log analyses for methane control are also presented as a complementary way of controlling methane and optimizing ventilation.

Jan 1, 2009

By J. S. Aronofsky, A. S. Lee

Results and experimental procedures are presented covering a preliminary laboratory investigation of the compaction of reservoir rocks and its effect on porosity and permeability. "Egective" compressibility was measured by subjecting core material to an overburden load and reducing pore pressure. Preliminary results showed relatively high effective compressibilities which varied from 1.8 X 10-' pore volume/pore volume /lb / in.~ for limestones to 7.7 X 10-6 pore volume/pore volume/lb/in.' and even higher for sandstone reservoirs conmining large amounts of fine clay material. It appeared from the initiol results that the effective compressibility of the material could not be correlated directly with porosity and that other factors, such as the amount of clay material present might have to be taken into account. The work to date has indicated that effective compressibility factors will have to be measured for the particular reservoirs heing studied. Oil permeability of core materid was measured under various effective overburden loads from 0 to 8,000 psi. It has been found than reduction of permeability in clean sands was relatively small and the results agreed with published data. For sandstones containing large amounts of clay, preliminary tests indicated very large reductions in permeability with increasing effective overburden pressures. It is planned to continue laboratory work on this phase because of its importance in low-permeability reservoirs. INTRODUCTION During the past few years increasing interest has been shown in the differences obtained between conventional laboratory measurements on cores at atmospheric pressure and those when reservoir conditions are simulated. Since these data are used as the basis for calculations and the prediction of reservoir behavior, their evaluation under reservoir conditions is most important. The effect of the compressibility of the rock upon the reservoir behavior, as fluid pressure declines, is an important subject which has been insufficiently studied. The total pressure on any plane through a reservoir is the result of the grain-to-grain rock pressure due to the overburden and the pressure of the interstitial fluid on the plane. As reservoir fluid pressure declines, the net load on the rock phase must increase, and the result is the compaction of the rock structure which produces a decrease in pore volume. Simultaneously, the rock grains expand into the pore space as the fluid pressure declines. The result of these two effects is defined as the effective compressibility of the rock and is expressed in units of pore volume/original pore volume/psi. Work was started by Imperial in 1954 on the effects of pressure on the porosity of reservoir rock. These studies established effective compressibility factors for some limestones and dolomites. Since that time effective compressibilities have been obtained for other Canadian and some U. S. reservoir rocks. In conjunction with the study of porosity, an investigation has also been initiated to determine the effects of pressure on permeability. One outcrop and several reservoir sandstone cores have been tested to date. ROCK COMPRESSIBILITY Apparatus anD Procedure A diagram of the apparatus used in these studies is shown in Fig. 1. The equipment was constructed so that a uniform pressure, equivalent to the reservoir overburden, could be applied by the oil pump and transmitted on all sides of the core through the Lucite sheath. The sample was completely saturated by circulating brine through the mounted core. Pressure in

By B. G. Hurd

A liquid-Freon permeameter suitable for making routine permeability determinations on small plug samples is described. The instrument is characterized by simplicity of design and ease of operation, and can be assembled inexpensively from stock items of laboratory equipment. It combines the unique advantages of both liquid and gas permeam-eters while eliminating many undesirable features of both general classes of instruments. Precision and accuracy of the specific Freon permeability measurements compare favorably with results of conventional liquid and gas permeameters. INTRODUCTION The liquid-Freon permeameter described in this paper was designed and built for special investigations of the properties of petroleum reservoir rocks. However, its unique advantages over conventional gas and liquid permeameters make it eminently suitable for routine use in the core analysis of small plug samples. In general, gas permeameters are very popular because of their simplicity of design and ease of operation. However, apparent gas permeabilities are normally higher than true specific permeabilities unless special corrections for gas-slippage are applied. Since these corrections require two or more measurements at different internal gas pressures, much of the time advantage of the specific gas permeability measurements is lost. Liquid permeability measurements, on the other hand, require no slip correction, but liquid saturation of the sample is time-consuming and difficult to insure. Often, too, the saturating liquid must be extracted or dried from the plug before it can be used for another experiment. The liquid-Freon permeameter described herein retains the inherent accuracy of the liquid permeameter, while affording the speed and ease of operation of the gas permeameter. Freon-12 (dichlorodifluoromethane) is a haloge-nated hydrocarbon, immiscible in the liquid phase with water. It may be considered as an oil phase in laboratory fluid-flow studies. Like butane and propane, it is a gas at normal room temperatures and atmospheric pressure but can be liquefied at moderate pressures. Complete liquid saturation of a plug sample with any of these fluids can be easily achieved by condensing vapors in the sample and operating the permeameter at a pressure above the vapor pressure of the liquid. Freon-12 is superior to butane or propane for permeability measurements only in that it is relatively nontoxic and completely incombustible, and thus presents no fire or explosion hazard. Freon-12 is believed to be completely nonreactive with most petroleum-reservoir rock samples. It volatilizes immediately from a plug on depressurizing to leave a clean, dry sample suitable for subsequent experiments. Thus, the inherent risk of changing physical properties of rock samples by cleaning and drying operations between permeability measurement and subsequent experiments is eliminated. DESCRIPTION OF THE PERMEAMETER The design of the liquid Freon-permeameter is shown schematically in Fig. 1. Regulated nitrogen pressure caps are used both as a driving pressure

By J. W. Whalen

The design, operation and evaluation of an instrumental method for the determination of fluid saturation in porous materials during multiphase flow studies is described. The presence of a magnetic tracer comprising approximately 20 per cent by weight of the aqueous phase is necessary in the utilization of this method. The results of two-phase relative permeability studies on two samples of porous media are presented. The realizable accuracy of the magnetic susceptibility apparatus in fluid saturation determination is estimated to be 2 saturation per cent. The response characteristics of the instrument, together with its simplicity of operation and maintenance are sufficient to recommend its use in routine laboratory studies. INTRODUCTION The problem of fluid saturation determination in laboratory multiphase flow studies characterizing oil recovery processes is one which has been given exhaustive study. Briefly the requirements which are considered essential for a method sewing this purpose are: (1) that the measurement be made external to the porous rock sample under study with no interruption of the fluid flow pattern, (2) the saturation indication must be independent of fluid distribution in the volume scanned by the measuring device, (3) the sensing element should have a small field of definition enabling several independent measurements to be made along the length of a core sample, and (4) the properties of the fluid phases should not be too greatly influenced by required tracer substances. Several methods of saturation determination have found acceptance as meeting most of the above requirements. These include resistivity,1 X-Ray absorption,2,3,4 radiotracers,5,6,7 and neutron8 diffraction The merits of these techniques have been discussed elsewhere in detail and will not be considered here. It may, however. be pointed out that the elaborate instrumentation required by the methods mentioned above, with the exception of the resistivity method (which is undesirable from other standpoints) encourages the development of methods of saturation determination which, while meeting the requirements listed, are simple to operate and maintain. It is the purpose of this paper to describe the instrumentation and evaluation of a technique for the determination of fluid saturation which combines relative simplicity with an accuracy and applicability comparable with methods now in practice. DESCRIPTION OF THE METHOD If the magnetic susceptibility of one of the fluid phases contained within a porous material is sufficiently greater than that of the remaining phases, the fluid content of the material with respect to that phase may be determined by an induction measurement. This may be accomplished by measuring the induced voltage in the secondary of a transformer which contains the porous sample and its interstitial fluids as a portion of the magnetic flux path. Since the fluids used in multiphase flow studies involving petroleum reservoir samples possess uniformly low magnetic susceptibilities, the use of a magnetic tracer in one of the fluid phases is required. In view of the relatively high paramagnetic susceptibility of salts of the transition elements, together with their high water solubility. these compounds become logical choices as tracer materials. Cobaltous chloride, having a specific magnetic susceptibility of 17.5 x 106 cgs electromagnetic units, was used as the tracer material throughout this work, in a concentration of approximately 20 per cent by weight in the aqueous phase. MAGNETIC SUSCEPTIBILITY APPARATUS Basically, the magnetic susceptibility apparatus consists of a transformer serving as the detecting head, a

Jan 1, 1955

By C. H. Stewart, R. A. Fagin

A two-dimensional, three-phasereservoir simulator was programmed for a large memory digital computer. It was designed to provide a practical solution to describing the complex physical relation between the natural forces and the physical properties of a heterogeneous reservoir when subjected to a specific set of conditions. A reservoir study is briefly described to illustrate application of the model. A full volumetric account of three phases (oil, gas and water) is performed simultaneously throughout an integration net representing the reservoir. Absence of one or two of the phases is treated as a special case of the more general situation. Expansion (or contraction) of all phases, including rock expansion, is performed so that the pressure calculation is the general unsteady-state case. To account for the large variations of subsurface elevation encountered in some reseruoirs, and to allow for segregation of the various phases, a gravity head term is included in the basic drive potential. Appropriate fluid and rock properties are used in polynominal surface form (functions of pressure and/or depth) or they can be entered as space variables at each position of the integration net. An unsteady-state water influx calculation, based on the method of van Everdingen and Hurst, l was connected to the boundary of the matrix to simulate aquifers of various sizes. In addition to reservoir calculations, three-phase fluid flow from the producing depth to the wellhead, including provisions for gas lift, was incorporated in the simulator. A workover routine was also built which can automatically switch to a different set of production relations when a gas-oil ratio or water fraction reaches a limit; or it can shut-in the well if prescribed. INTRODUCTION This paper describes a reservoir engineering mathematical simulator used to represent the complex interaction of natural forces and physical properties of a reservoir during natural depletion or with various injection schemes. The simulator, which was programmed for a large memory digital computer, is a two-dimensional calculation which handles three mobile fluid phases simultaneously (oil, gas and water). Basic requisites for the method are individual well production and pressure data, hydrocarbon fluid properties, geological data (producing depth and net sand), capillary pressure data, relative permeability data and permeability and porosity information. Matching the past performance of a combination drive reservoir often has yielded information concerning continuity and the validity of basic data. Detailed predictions of future performance can be made for continuation of current depletion methods (natural depletion) as well as for various types of recovery by gas or water injection. Combination injection cases and pattern studies can also be performed. Workover programs, gas lift and different types of artificial lift programs have been investigated using a technique similar to that described by Kern and Nicholson2 except that conditions of pressure and saturation at the block within which the well is located are used rather than average reservoir conditions. Drilling additional wells to optimize profit was explored, both as to number and location, by placing wells at different spots within the reservoir matrix. Special depletion processes can be examined, such as upstructure drainage and lateral (or strike) waterfloods in thin oil columns. In one case the mathematics of the simulator were modified to calculate the displacement in the vertical plane rather than in the horizontal plane. In this manner specific reservoir problems can be studied, such as coning of gas and/or water around production points, fingering along permeable stringers or, more generally, frontal advances in a heterogeneous section.

Jan 1, 1967

By B. P. Nutter, M. Lebourg, J. A. McAlister

The Multi-Flow Evaluator (MFE) is a new system of tools providing an original approach in drill-stem testing. It improves control during the testing operation and gives a more accurate evaluation of the fluid recovered while providing additional pressure information for reservoir analysis. The tools are operated entirely by up and down motion of the drill stem. This up and down operation provides a positive means of control and offers easily observed surface indications of tool operating position. An unlimited number of shut-in and flow periods may be taken with this tool while in the hole. The key to the success of the MFE system of tools is the safety seal packer. Until this development, operation of the tools by up and down motion only had proven unreliable. The success ratio of the MFE now exceeds the ratio achieved by conventional tools. The equipment includes a 2,750-cc chamber in which a representative sample of the flowing formation fluid is trapped at the end of the last flow period and brought to the surface under pressure. The sample can be evaluated at the wellsite or transferred under pressure for laboratory analysis. The sampling feature of the tool allows a sample to be obtained from a reservoir which has suffered minimum influence from production. A representative sample can be obtained for laboratory or empirical analysis by employing a testing technique to minimize drawdown during flow periods. Interpretation methods to take advantage of this additional information are presented and supported by actual field examples. INTRODUCTlON The Multi-Flow Evaluator (MFE) represents the most recent technical advance in formation testing. Drill-stem testing techniques prior to the MFE have been limited in three areas—(1) operaiion, (2) recovery analysis, and (3) reservoir analysis. The tools in the MFE system allow maximum control during testing and provide test information superior to that obtainable with conventional tools. The system of tools and the operation procedures are described to illustrate the mechanical advantages over present tools. Methods of pressure build-up and fluid

Jan 1, 1966

By Mohamed Mortada

A practical analytical method is developed in this paper which provides the practicing reservoir engineer with a handy method for analyzing oilfield interference problems. The procedure employed entails the principle of superposition of individual eflects. A set of charts has been developed and can be used to evaluate the individual effects. The charts represent non-steady state pressure distribution in an extensive aquifer due to a constant rate of production from an oilfield which can be approximated by circular geometry; they are constructed from solutions to the diffusivity equation considering a circular sink as distinguished from a point sink. The solutions are obtained partially from mathenlatical considerations and partially from analog and digital computations. By virtue of the method of images, the charts can be used also to treat the case of production near a fault. In addition, provisions are made for treating oil fields, which are produced by partial water-drive. Treating an actual field case is deferred to a subyequent writing; however, illustrative example is presented to demonstrate the method of application. The pressure distribution in linear and limited radial aquifers is currently under study. INTRODUCTION With the increasing difficulty of discovering new oil reserves, attention is necessarily focused on the efficient development and production of existing reservoirs. Efficient production depends on the availability of reliable means for predicting the pressure-production behavior of oil fields. The prediction is usually performed by analytical methods or computers. Analytical methods are often limited to systems which can be approximated by simple geometry. However, analytical methods have the distinct advantage of being readily performed without special equipment. Analog and digital computers lend themselves to a variety of problems, but are not always locally available and consequently their use is limited. One particular situation in which analytical techniques are needed involves studies of multiple, water-drive oil fields sharing a common aquifer. Such oil fields are in hydrodynamic communication. Therefore, production from any of the oil fields is accompanied by a pressure decline which is transmitted through the aquifer to the other oil fields and manifested as pressure interference. The rate of propagation of the pressure decline is such that the pressure may be significantly reduced many miles away from a producing pool. An example of pressure interference is found in the Woodbine formation in East Texas. There. a dozen or more pools, including several of major size, draw on the Woodbine formation for water drive. The pressure in the water-bearing portion of the formation has been significantly reduced over an area of several thousand square miles. This reduction of pressure has been such that when the Hawkins pool was discovered in 1940, the pressure in it had already been reduced 280 psi below the initial value in the aquifer. This resulted from prior production from other Woodbine sand pools, including, in particular, the East Texas field some 20 miles to the east of Hawkins.' The above example indicates that the interference effect can be significant and should be taken into con-

Jan 1, 1956

By R. N. Hoffman

A new technique for determining capillary pressure curves has been developed and tested. The technique differs from previously reported centrifuge techniques in that the centrifuge is slowly accelerated from zero to the maximum desired speed rather than being held constant at particular, progressively higher speeds. An important advantage of this technique over other methods of determining capillary pressure curves is the short time required to obtain the desired amount of data over the chosen pressure range. For example, to obtain 30 data points between 1.2 and 104 psig with a 1.55-in. long, 3/4-in. diameter core using a brine-air system, 6.6 hours were required with this technique. An equally important development of this paper is an analytic method for the conversion of the data from the centrifuge experiment to capillary pressure curve data. Previously there has been only an approximate conversion available. Although the capillary pressure curves determined by this technique appear to be as accurate as those determined by other techniques, the accuracy could be improved if certain variables, not treated in this experiment, were investigated. Among these are the dynamic distortion of the centrifuge equipment and imperfect initial saturation of the cores. INTRODUCTION Pirson 1 defines capillary pressure in a porous medium as "the differential pressure that exists between two fluid phases at their interfaces when they are distributed under static equilibrium within a porous material". Capillary pressure in rocks is known to be a function of fluid saturation, among other things, and a capillary pressure curve is defined for the purposes of this paper as a plot of the capillary pressure-wetting phase saturation relationship for a particular rock sample. Several methods are used for determining capillary pressure curves for small rock cores. Prominent among these are the semi-permeable barrier, mercury injection, and, to a lesser extent, centrifuge methods. The semi-permeable barrier method is currently the most popular. It features simplicity in both execution and the mathematical conversion of the experimental data into a capillary pressure curve. The main disadvantages of the semi-permeable barrier method are the time required — as long as two months — to obtain several points for the curve, and a fairly low maximum pressure before breakthrough of the non-wetting phase into the barrier occurs, for example, about 32 psig for a brine-air system. It is for these reasons that other methods such as the centrifuge method have been introduced. High accelerations and the absence of a barrier result in quicker attainment of saturation equilibrium at a given pressure. However, the centrifuge method involves much more expensive equipment and more difficult procedures and calculations than does the barrier method. The purpose of this investigation has been to improve the equipment and procedures of the centrifuge method and to develop an analytic method for the conversion of the experimental data into a capillary pressure curve. Hassler and Brunner did the original work in the determination of capillary pressure using a centrifuge. In their work the centrifuge speed was increased in a step-wise manner, each speed being held constant until saturation equilibrium was reached in each core. Saturation equilibrium was indicated when the volume of liquid collected in the graduated pipette of the core holder remained constant. According to Hassler and Brunner, equilibrium was reached in "a few minutes to one-half hour or more". In the centrifuge method, as opposed to the barrier method, the fluid saturation of the core is not a constant throughout the length of the core, but varies with the radius of centrifugation. Also, the capillary pressure cannot be read directly but must be calculated from a knowledge of the centrifuge speed and other parameters. For these calculations, Hassler and Brunner derived the following equations:

By D. L. Luffel

This paper presents two main developments: (1) a short, constantly accelerated centrifuge method of measuring drainage capillary pressure properties of cores, and (2) a new analytic method of treatment of centrifuge measurements to obtain saturation-capillary pressure relations. Toward the first objective, a clever experimental method is presented which is certainly more rapid than previous centrifuge methods. Toward the second objective, a new analytic method is derived to eliminate the approximate method of Hassler and runner. The derivation of this new analytic method, however, appears to be in error. The new analytic solution is developed as Eq. 5 in Appendix A. At the beginning of this derivation, the average wetting phase saturation of a core is expressed as, where is the centrifuge angular velocity. Eq. D1.2 makes explicit the dependence of S and S upon The second equation shown in Appendix A expresses the relation between the three variables Since for a given porous medium the saturation is uniquely related to capillary pressure, the functional dependence of S upon P is trivial, i.e., Pc.ri does not actually appear. Thus, the function 5 (r, ) may be expressed as a new function, S*, of PC, , alone The author transforms the variable of integration in Eq. Dl through use of Eq. D2 and writes

Jan 1, 1965

By A. M. Khan, I. Fatt

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

Jan 1, 1966

By R. F. Krueger, L. C. Vogel

A design is presented for an electrical analog computer which can solve non-steady state heat transfer problems in an extensive radial forrnation containing a moving cylindrical source. The computer is used to simulate a .simplified thermal oil recovery process in which heat transfer from a moving, constant temperature source takes place radially by conduction only. Temperature distribution curves are .shown for several different assumed modes of travel for the heat source. The data are used to estimate the residual fuel requirements necessary to maintain a self-propagating isotherma1 front for the particular system being studied. Although the computer is designed to represent an unique system. conversion factors can he adjusted to show the effects of changing the assumed values of system constants, such as thermal conductivity of the formation and temperature of the source isotherm. Examples are given to illustrate the relative influence the assumed values of these variab1es will have upon the quantity of fuel required to maintain the source. The maximum effect of heat transfer by the flowing air stream is estimated mathematically, and full requirements for this system are shown to he markedly reduced. Inasmuch as the examples discussrd refer to a simplified process, the results are not directly applicable to a practical field operation. It is believed, however, that the data il1ustr-ate some general trends which are important in thermal recovery processes. Data from laboratory and field experimentation can he used to modify the computer to take into account the influence of heat transfer by various mechanisms, for example, by injected air and by fluids ahead of the hot zone. With .such refinements it would be possible to estimate more accurately the limiting values for the rate of travel of the high temperature front and for the required air- injection rates. Use of this computer should aid in evaluating the economic feasibility of oil recovery by various thermal processes. It is hoped this paper will stimulate further work by others to help accomplish this objective. As the costs of exploration for and development of new oil fields rise. increasing the recovery of oil from established fields becomes more and more important. In fields containing low gravity crude oil, ultimate recoveries by primary means are sometimes less than 10 per cent and often total only 10 to 20 per cent of the oil in place. Conventional secondary recovery methods are not always effective in increasing this total. Methods which have been proposed for increasing the recovery of viscous oils include thermal processes for increasing the mobility of the oil. One method is to provide a moving heat source generated within the formation to heat the oil and a flowing gaseous medium (for example, air) to sweep the oil to the producing wells. The possibility of such a process has intrigued engineers for many years.'. '." but little or no information has been available in the literature relative to the technical or economic feasibility of this method. Recently, however, interest has been revived by publication of laboratory and field test results indicating the feasibility of maintaining

Jan 1, 1956

By C. Scala, F. Bernstein

General expressions are derived relating the streaming potential to the electrochemical potential of a permeable junction separating two electrolytic solutions. By the methods of irreversible thermodynamics, these potentials are related to the flux of mobile species arising from differences in pressure and chemical composition. The flux is conveniently expressed as a transference number, the equivalents of a mobile species transferred through the junction per faraday of electricity. At constant pressure, the flux equation gives the electrochemical potential across the permeable junction. The streaming potential across the junction is obtained by imposing the condition of constant chemical composition. The familiar expression for the borehole SP is shown to be a particular solution of the general flux equation. The derived exprersions were verified in experiments with seven shale samples. Both the theory and experimental data provide additional proof of the presence and magniturle of the streaming potential in shales. INTRODUCTION In the interpretation of electric logs, there has been much concern over the correction to be made to the recorded SP for an electrokinetic or streaming potential originating from the difference in pressure between the mud column and the formation. It was long recognized by M. R. J. Wyllie that such an electrokinetic potential source mas in the mud cake opposite permeable sands.' More recently, M. Gondouin and C. Scala2 showed an additional but opposing electrokinetic source of potential to be in shale beds. Their work was supported by the very extensive and excellent work of Hill and Anderson.3 n the present contribution, the general relationships between the electrochemical potential and the corresponding streaming or electrokinetic potential are examined and clarified. A notable advance towards understanding the relation between the various sources of electrical potential in a borehole has been the application of the now common methods of irreversible thermodynamics to describe the passage of ions and water through charged electrolyte junctions such as in shales.' The advantages of these methods are twofold. Unlike classical thermodynamics, they can be applied to non-equilibrium sys- tems, i.e., systems in which flow takes place, provided the flow processes have reached a steady state. They also permit conventional thermodynamic properties to be related to the steady-state dynamic properties. These methods arc used here specifically to derive expressions relating the potential behavior of shales in terms of ion and water transport. Experimental verification of the expression for the shale streaming potential is given which in turn offers theoretical support to the existence of a streaming potential in shales.' Importantly, the relative contribution of transport parameters to the SP is clarified. The theoretical treatment is essentially that of Lori-mer, Boterenbrood and Hermans3 ut has been extended to derive the general expression for the SP.

By Harold B. Janzen, M. A. Nobles

Campbell and Martin1 described a new exploratory logging tool which was developed as the result of a new concept of the phenomenon that takes place when a porous formation containing mobile hydrocarbons and irreducible interstitial water is invaded by mud filtrate. They postulate that invasion will create a bank of interstitial water between the flushed zone and the undisturbed formation. This postulate was confirmed by Wyllie.' Thus, if R,, is greater than R, and mobile hydrocarbons are present, the resistivity profile with increasing radii from the borehole is high:low:high, as shown in Fig. 1. If mobile hydrocarbons are absent, the bank of interstitial water is indistinguishable from the undisturbed formation and the resistivity profile, when R is greater than R,, is high:low, as shown in Fig. 2. The zone containing the bank of interstitial water and residual oil was named the "low zone". Campbell and Martin conducted extensive electrolytic wedge and resistivity analog model studies to select a logging system that would reflect the presence or absence of the low zone. The system they selected was that of a family of lateral curves. The model studies indicated that the presence of a low zone could be detected in approximately 60 to 70 per cent of the producing sand reservoirs with a family of lateral curves provided that the sonde effectively occupies 90 per cent of the diameter of the hole. These studies were conducted with D, = 5d. The presence or absence of the low zone is best determined by a plot of the apparent resistivity vs spacing as shown in Figs. 3 and 4. A dip or inflection in the profile that repeats at successive depth intervals is considered diagnostic. To reduce the shunting effect of the hole and mud, a sonde of insulating material that is surrounded by a series of flexible rubber bumpers was designed. The bumper efficiency is considered to be the ratio of the bumper diameter to the hole diameter. Field tests conducted by Campbell and Martin on 10 wells with a tool that provided for station readings indicated that dips or inflections could be detected in 19 out of 25 known producing zones. INSTRUMENT DEVELOPMENT Halliburton developed a logging system to simultaneously and continuously record eight lateral curves. The instrumentation is described in the literature.' The original logs made

By C. D. Stahl, S. M. Farouq Ali

This investigation attempts to describe and simulate the alcohol displacement process by means of a cell model, as employed in chemical engineering practice. The proposed model is more simple than previously proposed models. and utilizes parameters chosen on a theoretical basis. The model successfully reproduced the formation of the stabilized bank and the breakthrough of alcohol, the latter depending on one of the model parameters, which may be correlated with the length of the porous medium. Moreover, the effects of the phase behavior of the liquid system involved, as observed in experimental studies, were reproduced. Several variations of the basic model were devised and tested on a digital computer. These included the cases in which: (1) the actual value of fractional flow was used in cell-to-cell computations; (2) the number of cells was varied within the same run; and (3) incomplete rather than complete phase equilibrium was assumed within each cell. The proposed cell model clarifies the basic mechanism of the process. Detailed concentration profiles obtained for each cell, for instance, showed the mechanism of bank formation in relation to the phase behavior characteristics. The results obtained indicated a varying degree of phase equilibrium concommitant with changes in the velocities of the phases in an actual alcohol displacement. This condition was approximated by changing the number of cells during the simulation. Interesting information was obtained on the influence of path length on the efficiency of alcohol displacement, which has been the subject of some controversy. Certain limitations preclude the use of the proposed model as a substitute for experimental studies. The results obtained were, nevertheless, of value in interpreting the experimental results. INTRODUCTION During recent years considerable effort has been directed toward an understanding of alcohol displacement, the process whereby oil and water are recovered from a porous medium by the continuous injection of a solvent. The complex nature of the physical process involved has so far defied a complete mathematical treatment. Other methods of approach, amounting to an overall material balance, have been proposed, yielding useful information on certain aspects of the process.1-3 Taber et al, in particular, defined the displacement mechanism in terms of the phase behavior of the alcohol-oil-brine system invoIved.2 Wachmann reported a mathematical treatment of alcohol displacement subject to certain simplifying assumptions? Donohue proposed the use of a "cell model" for simulating alcohol disPlacement.5 The nature of the assumptions involved limited the utility of the model. The present work attempts to examine the variables involved in the simulation of alcohol displacement and discusses several possible versions of the basic cell model. Under certain conditions the model results are similar to the experimental results. In particular, the spontaneous formation of of the stabilized bank and the effects of the system phase behavior were successfully reproduced. PREVIOUS WORK ON CELL MODELS Cell models and the theoretical plate concept are often used in solving chemical engineering problems in which an explicit mathematical solution may be difficult or impossible to obtain. Examples of such applications occur in distillation, gas-liquid chroma tography,6 reactor technology, absorption, etc. In petroleum engineering, such a model was used by Attra7 to describe non-equilibrium gas drive, and by Higgins and Leigh ton8 to calculate sweep efficiency in water flooding. Aris and Amundsen pointed out the equivalence between the diffusion model and perfectly mixed cells connected in series.9 Deans proposed a three-parameter cell model to simulate two-component

Jan 1, 1966

By A. H. Heim

A method and apparatus for measuring gas porosities of rocks are described. The apparatus can be assembled from commercially available components. In principle, measurements are made by volume substitution at constant pressure. The maximum error is not more than 0.3 porosity per cent. Typical results are given. INTRODUCTION Determining the porosity of rock samples is one of the most important and yet most varied types of measurement in core analysis. Among the many techniques devised are the so-called "gas porosity" methods. An old and well known example is the Washburn-Bunting method.' The U. S. Bureau of Mines2-' described and later improved the apparatus for a now widely used method generally known as the "Boyle's law" method. In the present form of the Washburn-Bunting method,' the volume of air in the pores of a rock sample at atmospheric pressure is extracted and then collected in a graduated burette at atmospheric pressure. The volume of air is read directly as the pore volume of the sample. The absolute error in reading the collected volume of gas is independent of the total volume; thus, the relative error is larger when the volume is small, as it is for rocks of low porosity. In addition, the sample after measurement contains mercury, which limits its use for other analyses. The Bureau of Mines (or Boyle's law) method measures directly the solids volume of a sample from which the pore volume and porosity are derived, using a separate measurement of the bulk volume. Gas at a few atmospheres pressure is introduced into a sample chamber of known volume containing the rock sample. The pressure is accurately measured. Following, the gas is expanded into a burette at 1 atm, and the gas volume is read directly. From the initial pressure p, and the final pressure p2 and volume v,, the initial gas volume v1 is calculated using Boyle's law; that is, p1v1 = p2v2. Volume v, minus the volume of the empty sample chamber is the solids volume of the sample. The accuracy of the method is limited, unless corrections are made, by deviations of the gas from the "ideal" gas-law behavior assumed in the simple form of Boyle's law. The purpose of the present paper is to describe a method for measuring the gas porosity of a rock which avoids many of these difficulties. Gas volumes are measured directly with the same accuracy as the bulk volumes. Pressures of at least an order of magnitude larger than those of previous methods are employed to insure rapid penetration of the gas into the sample. While special equipment may be built to apply the method, the porosimeter may be constructed as well from commercially available components. For simplicity, the apparatus described will be referred to as the "Constant-Pressure gas porosimeter". THE CONSTANT-PRESSURE METHOD Fig. 1 shows schematically the arrangement of components comprising the present Constant-Pressure porosimeter. Briefly, the method is one of volume substitution and may be considered a null measurement. Omitting (for the present) some of the operational details, the method of measurement consists of the following three steps. 1. After evacuation, the volume of the measuring system (a ballast chamber, a manifold, two gauges and their connections) up to the sample chamber is filled with gas to a high pressure (- 1,000 psi). A sample of the gas at this pressure is trapped in one side of a sensitive differential pressure gauge to serve as the reference pressure for subsequent steps. 2. The evacuated sample chamber containing the rock sample is opened to the measuring system. As the gas expands into the chamber, the resulting decrease in pressure unbalances the differential pressure gauge. 3. The pressure is restored by means of a mercury volumetric pump. The volume of mercury injected exactly equals the free or void volume of the sample chamber (volume of empty chamber minus the solids volume of the rock within). From the injected volume and the known empty chamber volume, the solids volume is obtained and the porosity calculated. The pressure and the volume occupied by the gas are the same before and after opening the sample chamber. Expansion and compression of the gas are incidental operations and do not enter into the calculation of porosity. By the pressure balancing or nulling, the free volume of the sample chamber is merely substituted by an equal and measured volume of mercury. Since the measurements are at constant pressure, there are no compressibility corrections necessary for the sample chamber.

By C. S. Matthews, M. J. Fischer

Distortions caused by dip on five-spot sweep pattern were determined by visual experiments on a model composed of parallel plates of glass with holes drilled at the well positions. This model war designed to simulate a reservoir depleted of its primary oil and containing a gas phase. Prior to interference, it was found that the injected fluid took a nearly circular shape about the injection well for the cases studied. However, these circles moved down-dip considerably under the influence of gravity. It was concluded that pattern flooding should be practical for many dipping reservoirs. A quantitative measure is given of how much the well position should be altered to allow for dip. Experimental deternzinarions of injectivity showed this to be about the same in a dipping reservoir as in a horizontal one. INTRODUCTION This investigation was undertaken at the suggestion of the Pacific Coast area, Shell Oil Co., where a pilot water flood in a dipping reservoir was being carried out at the time this paper was in preparation. Prior to the undertaking of this flood there had been considerable speculation as to how gravity would distort a flood pattern in a dipping reservoir. The thought was expressed in some quarters that a five-spot flood would not be applicable to a dipping reservoir because the flood pattern would elongate in the direction of dip and would never form a "seal" around each five-spot. If the flood pattern elongated sufficiently, both the oil and the water might flow largely downdip rather than into the pattern producers. In view of the importance of the California water floods, it was deemed important to determine the effects of gravity on a five-spot pattern in a dipping reservoir. Previous work had also been carried out in this laboratory,1 but only for the single-fluid case. The present investigation extends these results to multifluid cases. Several sets of pictures of flood fronts are presented so that the process can be visualized. From these pictures, curves are obtained which show the correct drilling position for achieving equal times of oil-bank breakthrough. Several injectivity curves are also presented which may aid in prediction of injection rate. ASSUMPTIONS The present results were obtained with a model composed of two parallel plates of glass spaced about 1 mm apart. When the inertial forces are small, as in these experiments, the flow between the plates follows Darcy's law and is analogous to flow through porous media. In using these plates, it is not possible to scale effects of capillarity correctly as regards a transition zone between injected and displaced fluid. Calculations show, however, that effects of such capillarity should not be important in determining the gross flooding pattern in the field. In fact, in this model the capillary forces tend to circularize the pattern to a greater degree than in the field. To preserve a large ratio between the viscous forces and the circularizing capillary forces in the laboratory model, viswus glycerine was used as the flooding liquid and the rate of flooding was very high. The circularizing force is believed negligible in these experiments. In addition to the effect of the surface tension y, which tends to circularize the pattern, another capillary force must be considered. This is the adhesion tension y cos, which acts at right angles to the circularizing term y. If the plates in the model are "wetted" the same all around the periphery, then the adhesion tension,

Jan 1, 1957