Search Documents

By E. H. Timmerman, A. F. van Everdingen, J. J. McMahon

The prevent paper contains a method which combines the material balance equation' with the water influx equation' to obtain reliable values for the active oil originally in place and a quantitative evaluation of the cumulative water influx. The method is illustrated by an application to a reservoir without original gas cap. In the absence of an original gas cap. results may be obtained using only field production and pressure data, PVT analyses. and a minimum of subsurface information. CHARACTERISTICS OF RESERVOIR AND RESERVOIR FLUIDS Production in the field under review is obtained from the top of the Wilcox formation of Eocene age, at a depth of approximately 8,100 ft subsea. The structural map, Fig. 1, shows that the accumulation, half elliptical in shape with the long axis in east-west direction, is trapped to the north and on the upthrown side of a normal fault. A number of east-west faults with additional minor faulting in random directions are found in the general area, and seismic information and production performance indicate in this particular case the existence of a second fault a short distance downdip, with a strike more or less parallel to the fault controlling the accumulation. It is estimated that 1,830 acres were originally underlain by oil and that the maximum thickness of the original oil column, about 37 ft, was only slightly greater than the maximum thickness of the oil-bearing sand. Maximum net sand thickness was estimated at approximately 26 ft. From logging and subsurface information the gross sand volume has been placed at 37.400 acre-ft. which reduces to about 27,500 acre-ft after probable non-productive intervals have been deducted. The following table presents laboratory analyse; data for the Wilcox sand cores. Average Average Samples Considered Permeability Porosity All Samples 236 md ..19.9percent Samples with k>l md 256 md 20.5 per cent Samples with k>8 md 275 md 20.9 per cent Interstitial water was estimated at 15 per cent. The initial gas-oil ratio was approximately 900 cu ft/bbl. and the initial shrinkage factor was determined by the laboratory to

Jan 1, 1953

By W. J. Joslin

The frontal-advance equation can determine how the fluid withdrawal rate and subsurface operating pressure influence oil recovery from pressure-maintained reservoirs having characteristics favorable for efficient vertical segregation. For effective vertical segregation in formations having low dip angles, the requirements are: thick, massive sands, low viscosity oil and high vertical permeability with few discontinuous shale members. This approach has been successfully applied to the LL-370, which is a massive Eocene sandstone reservoir containing crude with an average gravity of 26" APT, a dip of only 3" and vertical permeability of ½ darcy. The average gas cap saturation was calculated for various rates of vertical gas-oil contact movement corresponding to a wide range of reservoir producing rates. This method also evaluated the effect of different pressure maintenance levels on oil recovery. Applicability of this procedure for calculating the present oil saturation in the LL-370 gas cap has been substantiated from field performance, where the oil saturation was computed by comparing oil recovery with the gas cap volume vacated. INTRODUCTION During the natural depletion history of the LL-370 reservoir, located in Lake Maracaibo, Venezuela, excellent vertical segregation resulted in the formation of a thin gas zone over a substantial portion of the oil-productive area. Subsequently, gas injection forced this thin gas-saturated zone slowly downward, reducing the oil saturation in the gas cap to a low value. The oil recovery has been raised considerably above earlier predictions where vertical segregation effects were under-estimated. In order to determine the maximum efficient producing rate and optimum operating pressure for this reservoir, a vertical gas frontal-advance analysis has been made through application of the fractional flow equation. LL-370 RESERVOIR CHARACTERISTICS The reservoir consists of a truncated monocline of massive. highly-permeable Eocene sandstone, called the B-b formation, covering an area of more than 13,500 acres. The original oil in place was 2.17 billion STB. Internal non-sealing faults subdivide the structure into four blocks called the B-6-x.6, 10, 11 and 13, which are effectively bounded by major faulting on the flanks, an unconformity updip to the west, and an immobile aquifer to the southeast. The producing formation occurs at an average depth of 5,250 ft subsea and has a closure of more than 1,500 ft. A typical producing section has a net oil sand thickness of 170 ft. Scattered throughout the section are a few randomly-positioned shale lenses which increase in thickness and areal extent downdip. The reservoir rock and fluid properties vary with depth. The permeability and produced oil gravity are 1600 md and 28" API at the top, compared to 300 md and 18° API at the water-oil contact. The productivity of updip wells is, therefore, much higher than the reservoir average. The reservoir oil was originally saturated at the crest, grading to a slight degree of undersaturation at the water-oil contact. There was no original gas cap; and material balance and the field performance indicate no water influx. Except for the northern B-6-x.6 block, the fluid and rock properties have already been published.' Average LL-370 reservoir properties are shown in Table 1. RESERVOIR PERFORMANCE UNDER NATURAL DEPLETION Fig. 1 shows production performance from 1939 to date. Early development was gradual, with 30 wells pro-

Jan 1, 1965

By J. C. Trantham, J. W. Marx

From 1955 to 1958 the Phillips Petroleum Co. conducted a series of small scale counterflow combustion field tests in a tar sand about 60-ft deep and 6 to 12-ft thick near Bellamy, Mo. A total of seven different well patterns, including conventional five-spots and seven-spots plus a 15-well line drive pattern and a 10-well radial pattern, were employed. Electric heaters, gas burners and wellbore fuel packs were tested as counterflow ignition devices. Direct drive ignition followed by reversal to counterflow operation was not possible in this case because of formation plugging by the semi-solid tar during the direct drive phase. Sustained counterflow burning was achieved with wellbore fuel packs under conirolled conditions which permitted close correlation of fire front velocity, air velocity, air-oil ratios and yields for a wide variation in process parameters. The tar in place was about 10' APl with a reservoir viscosity greater than 500,000 cp. The composite oil prodlcct was 26' API with a viscosity in the range of 5 to 15 cp. Vertical sweep efficiency in the line drive pattern was essentially 100 per cent, with no evidence of gravity segregation. Original tar in place ranged between 800 and 1,000 bbllacre-ft, specific rock permeability averaged about 800 md before burning and effective air pemeability with residual oil and water in place (at the time of ignition) averaged 250 md. The volume of oil recovered during stabilized line drive counterflow burning represented about 67 per cent of the volume of tar originally in place in the burned out zones. A minimum was found when air-oil ratios were plotted against formation air velocities, and a limiting air velocity existed below which the counterflow fire front reversed and burned back over its own trajectory, using the residual coke left behind. INTRODUCTION Early in 1955 the Phillips Petroleum Co. began preparations to field test counterflow underground combustion as a method for producing oil from high viscosity hydrocarbon deposits, such as bituminous (tar) sands. At that time all published information on the use of underground burning io increase oil recovery was based on the direct drive process, in which the fire front moved through the formatim in the same direction as the injected air stream. Such direct burning is not applicable to many heavy oil sands or to most tar sands because the heat-thinned native hydrocarbon congeals in the cold rock ahead of the fire front forming a gas permeability block which smothers the combustion. Laboratory experiments first demonstrated that this permeability blocking would not occur if the fire front and the air stream moved in opposite directions. In this counterflow (reverse) burning process the heated heavy oil or tar is driven back through the combustion zone. Here thermal cracking produces relatively light hydrocarbon products which pass through the heated rock behind the fire front (largely in the vapor phase) and are incapable of causing a gas permeability block.' Exploratory wring revealed that consolidated bituminous sandstones with a combination of tar saturation, permeability, thickness and depth favorable for experimental field testing existed at several points in western Missouri. By the fall of 1955, a test site for determining technical feasibility of the counterflow process under field conditions had been selected near Bellamy, Vernon County, Mo. — about fifty miles north of Joplin. The pay zone was too thin and too shallow to constitute a commercial prospect, but these same factors permitted controlled experimental operation with a small compressor installation and with a large number of monitor wells. These were needed for adequate observation of the underground phenomena at reasonable cost. RESERVOIR CHARACTERISTICS The reservoir chosen for the experimental field tests was a 124 thick section of tar sand extending from 49 to 61 ft subsurface, with shale and siltstone laminations sealing the top and bottom. This section was part of a larger 30-ft thick tar sand deposit which extended both above and below the test zone. Core analyses showed that the 12-ft sand was effectively subdivided into two approximately equal sections, with the lower zone being consistently more permeable than the upper zone. The most significant laboratory core data are summarized in Table 1. The line drive test which is the principal subject of this report utilized only the lower 6-ft thick zone, with the upper zone regarded as part of the 55-ft overburden. WELL PATTERNS The over-all field test program involved a total of seven different well patterns, including conventional five-spots and seven-spots plus a 15-well line drive pattern and a 10-well radial pattern. In some of the tests the full 12-ft zone was employed; in other cases only the more perme-

Jan 1, 1967

By J. Jones-Parra, R. S. Reytor

The porosities of fractured limestone reservoirs can be divided into two broad types in accordance with their effects on fluid distribution and fluid flow. In the coarse porosity, gravity segregation lakes place freely and the resistance to fluid flow is very small. In the fine porosity there is no segregation and a high resistance to flow, and it has relative permeability characteristics similar to tight sandstones. By analyzing the affect of the two porosities it is concluded that in some cases to recover the maximum amount of oil it is necessary to remove large quantities of gas from the reservoir by producing at high gas-oil ratios. In this manner the fine porosity is drained of its oil and the gas-oil contact drops slowly permitting higher production rates from the oil leg. Since this conclusion is contrary to widely accepted principles of conservation, a mathematical model (Fig. I) was constructed to duplicate the conditions desired. The behavior of the model indicates that under certain conditions it is possible to recover more oil by producing at high gas-oil ratios than by production at low gas-oil ratios, and that the rate of production is affected more by gas shut-offs than by the decrease in pressure (Figs. 2, 3, 4 and 5). INTRODUCTION It has been noted that certain limestone reservoirs behave in such a way as to indicate that there are two distinct types of porous spaces available for fluid storage and fluid flow. Regardless of the nature of these porous spaces the distinguishing characteristic is that in one type of porous space there is definite gravity segregation between the oil and the gas; while in the other, the gas evolved tends to remain distributed equally throughout the reservoir. For the sake of simplicity the porosity in which gravity segregation takes place at a fast enough rate to affect the behavior will be called "coarse po- rosity", while that in which this is not so will be termed "fine porosity". The behavior of a reservoir of this type would be different than that of a sandstone reservoir, especially if most of the oil is in the fine porosity. It has been observed that fluid flow to the wells takes place mostly through the coarse porosity. The large fractures which

By J. G. Giddings, R. Kobayashi

Residual viscosity, the ViSCOSitY at a given pressure and temperature minus the dilute gas viscosity at the same temperature, has been found to be independent of temperature for pure components and mixtures for simple molecules and light hydrocarbons. The residual viscosities for methane, ethane, propane and four mixtures were correlated as a function of reduced density and molecular weight. INTRODUCTION Earlier efforts in correlating the viscosity of light hydrocarbons and associated gases under pressure, primarily by some form of the corresponding states principle, have been reported elsewhere.7 The formal theory of corresponding states principle for viscosity from non-equilibrium statistical mechanical arguments for simple molecules has been developed by Helfand and Rice.13 The most serious obstacle to the development of a generalized corresponding states representation of viscosity data, even for simple molecules, has been the limited quantity of precise viscosity data of pure components and mixtures over a wide range of pressure and temperature conditions. Viscosity data on light hydrocarbons have recently been contributed by several laboratories.2,6,8,23,25,26 In the study of the viscosity of nitrogen, Michels and Gibson" observed that the viscosities of nitrogen as a function of density gave parallel isotherms. The same behavior has since been reported for other gases.".= The observed behavior indicates that the residual viscosity or the viscosity at a given pressure and temperature minus the dilute gas viscosity at the same temperature is a unique function of the density. Abas-Zade,1 Michels and Botzen,18 and Michels, Bot-Zen, Friedman and Sengers,19 observed similar behavior for thermal conductivity, the residual thermal conductivity being defined in a manner parallel to that of residual viscosity. Thodos and co-workers 4,15,24 applied the residual viscosity concept to the viscosity of pure monatomic and diatomic gases. Ellington and co-workers10,25 made the residual viscosity vs density plots for ethane and propane, but omitted the region below a residual viscosity of 15 micropoise and showed a systematic separation of the ethane isotherms in the low density region. Eakin correlated the residual viscosity of methane, ethane, propane and n-butane empirically as a function of density and molecular weight: Carmichael and Sage6 recently presented ethane residual viscosity vs density values in the gaseous liquid, and near critical regions as a single univariant function of density. Flynn, Hanks, Lemaire and ROSS" have contributed extremely accurate viscosity data to confirm that the residual viscosity concept is essentially correct, but may not be accurate within 0.1 per cent of the viscosity value for the gaseous state- The recent studies by Giddings12 indicate that the residual viscosity of methane, propane and each of four methane-propane mixtures are unique functions of density Over the temperature and the density ranges studied.

Jan 1, 1965

By R. H. Adams, A. M. Khan

The Huntington Beach oil field is on the Pacific Coast, southeast of the City of Los Angeles (Fig. 1). It lies both onshore and offshore. Most of the onshore section is operated by Standard Oil Co. of California, WOI. The Tar reservoir is of upper Pliocene age, about 200 ft thick, and consists principally of fine-grained, unconsolidated, somewhat shaly sand, with thin inter-bedded shale stringers throughout (Fig. 2). The reservoir is a broad, highly faulted, elongated anticline, with an areal extent of about 370 acres. Most of the faults bisecting the structure are normal, trend north-south, and are the main cause of reservoir discontinuity (Figs. 3 and 4). From 1964 to Dec., 1967, 259 steam stimulations were performed on a total of 125 Tar sand wells. Eighty-six new wells were drilled and steam stimulated. Of the total, 96 have been stimulated twice, 34 three times and 3 four times. During this same period, a two-injector steam drive pilot has been in progress. The response obtained from this project has been compared with the response to cyclic stimulation. This paper presents an analysis of the field observations for the two methods of thermal recovery. History of Development Drilling to the Tar pool commenced in 1924 with an initial well production of 40 to 50 B/D of 12o API gravity oil. Complete primary development of the pool has been spread over several years. Approxi- mately 40 wells were drilled during the first phase of development, which was completed in 1945. The subsequent primary development started in 1960 with initial well productions of 30 to 40 B/D oil. Thirty-seven wells were completed in the Tar pool during this phase. For the most part, slotted, gravel flow packed liners have been used in the wells in this pool. A cyclic steam stimulation pilot was initiated in March, 1964. Based on the very successful stimulation in the pilot area, the cyclic stimulation project was expanded to include all of the Tar pool wells producing at the end of 1965 (Fig. 5). Also, a development drilling program was begun at this time. To date, 86 new wells have been drilled under this program. The number of wells drilled each year in this secondary phase of development is also shown in Fig. 5. New Well Completion Practices Most of the new wells have been completed to produce from both the Upper and Lower Tar zones. Generally, a casing stub is cemented from the top of Lower Tar to the base of clean Upper Tar sands to exclude intermediate water-bearing sands. The water string is then cemented from the top of Upper Tar to about 1,500 ft with neat cement and 30 percent silica flour. The annulus from 1,500 ft to surface is then filled with a petroleum-based gel through a high-temperature port collar. Once again, a history of sand-free production dictated the use of 60-mesh gravel flow packed liners. The Upper and Lower Tar

Jan 1, 1970

By J. B. Campbell, V. V. Valleroy, B. T. Willman, L. W. Powers

A steam drive of heavy oil was field tested in a shallow, low oil-saturation formation near Deerfield, Mo. The pilot was conducted in the Warner formation, a sandstone containing an 18' API oil having 1,000-cp viscosity at the 60F origind reservoir temperature. The formation war at a depth of 160 ft. Steam was injected into nine input wells arranged in an array of inverted five-spot patterns. In the completely confined center pattern, 14 temperature observation wells were installed to obtain thermal data and observe test progress. Late in the test, slugs of ammonia were injected to trace the flow paths of injected fluids. From the test area about 7,000 bbl of oil were produced. Data were obtained on areal and vertical temperature distribution, steam front advance, reservoir fluid movement and terminal saturations. This field test of a steam drive (I) demonstrated the feasibility of the method, (2) confirmed that the low residual oil saturations observed in the laboratory are obtained in the steam-swept region in the field and (3) provided recovery and conformance data for one set of field conditions. INTRODUCTION The Deerfield steam drive pilot test was conducted in a shallow sandstone containing 1,000-cp oil. The venture was undertaken cooperatively by the research and production departments of Carter Oil Co., which organizations have since been consolidated into Esso Production Research Co. and Humble Oil & Refining Co.. respectively. The production department was interested in steam injection at Deerfield because it appeared to be the most promising method of commercially producing this heavy oil deposit. The research department was interested in applying the new recovery method and in evaluating its performance in the field. At the time the test was begun, the initial oil saturation was not well known. Subsequent air coring and early pilot results confirmed that there was too little oil in place for profitable commercial exploitation by steam. Pilot termination at that time, however, would have been premature for evaluating field performance of the process, and the tert was continued to obtain additional data on steam injection as a recovery method. The test was located in Vernon County, Mo., about 10 miles north of the town of Deerfield and only a few miles from the Kansas border. The pilot site was selected as typical of the area. The location represented neither the highest nor the lowest oil saturation region in the acreage under lease in 1954. The steam drive was conducted in the Warner sandstone of Lower Pennsylvanian age. At the test site the top of the Warner occurs at about 160 ft subsurface and the formation is a fine- to medium-grained micaceous sandstone that dips gently to the northwest at the rate of 12 to 15 ft/ mile. A cross-section and permeability profile of the test location are shown in Fig. 1. At the pilot location the average total thickness of the Warner formation is about 43 ft, but the effective thickness for steam drive is 26 ft. Two distinct types of hydrocarbon saturation are apparent. The lower portion of the total sand, averaging about 17 ft thick, contains a very heavy asphaltic material that will not flow under the influence of a steam drive. This bottom interval, referred to as a dead oil residue, was not considered as part of the net sand undergoing steam exploitation. The initial formation and fluid properties of the upper 26 ft in the test area are summarized in Table 1, and variation of oil viscosity with temperature is shown in Fig. 2. Imbibition tests on preserved core samples taken at the end of the pilot test showed that the Warner sandstone was then neutral or slightly water-wet. Initially, the reservoir may have been strongly water-wet as indicated by low relative permeability to water during both water injection testing and early steam injection. PRIOR HISTORY Initial production tests of wells at the pilot site produced water with only a faint show of oil. No gas was produced except at Well 7-W in the pilot area and at another well about 1/3 mile northeast of the pilot. Prior to the start of the steam drive, a two-well water injection test and a two-well air injection test were conducted. No oil was produced by either. Water was pumped into Well I-W in the northeast corner of the pilot area with simultaneous production from Well 1 (Fig. 3). The air-injection tat was run at input Well 9-W and its offset, Well 2, in the southwest corner. Air and water injectivities were about the same when corrected for viscosity and pressure differences.

By H. J. de Haan, J. van Lookeren

In view of their size and favorable characteristics, Shell's heavy-oil fields on the eastern coast of Lake Maracaibo, Venezuela — known as the "Bolivar Coast" — are obvious candidates for the application of thermal recovery methods. During the last decade various methods have been tested, such as steam drive,', steam soak and in-situ combustion.3 Steam soak was discovered as a promising production method rather accidentally in 1969, during early steam drive testing in the Mene Grande Tar Sands1 When steam erupted at the surface due to breakdown of the overburden, the injection wells were back-flowed to relieve the reservoir pressure. This resulted in high oil production rates, all the more impressive because the reservoir is unproducible by primary means. It was concluded that injection of limited amounts of steam might be a very effective method for stimulation of heavy-oil wells. To obtain information that would be more widely applicable, testing of the process was continued in the eastern part of the Tia Juana field (Fig. 1). Here even more favorable results were obtained than in Mene Grande. This created sufficient confidence to prompt a large-scale project, which was considered essential for an adequate evaluation of ultimate recovery, steam requirements, mechanical problems and economic prospects. To obtain an early answer on the ultimate recovery a fairly depleted area was selected in the western part of the same field ("A" in Fig. 1). The first steam was injected in March, 1964. This paper describes the performance until the end of 1966. The production mechanism and well inflow performance are discussed in some detail. In view of the favorable results obtained so far, steam soak operations have expanded considerably since the start of this project. Presently, five other projects are under way, of which two are about the same size as the one under discussion, with further extensions being programmed. Description of the Project Reservoir Data Shell's part of the Main Tia Juana reservoir has a proven area of 8,800 acres and an average net oil sand thickness of 130 ft. It contained an initial quantity of 2,750 million bbl of stock tank oil (STOIIP), with an API gravity ranging from 10 to 15 API and a viscosity of 100 to 10,000 cp (at 122F). As a geological unit it is known locally as the Lower Lagun-illas member of the Lagunillas formation, which is early Miocene in age. The sand bodies are believed to be distributary channel sands deposited in a lower coastal plain. Most of the channel sands have lateral and some vertical communication. The field was developed mainly in the period from 1936 to 1950. The ultimate primary recovery is estimated at 18.1 percent of STOIIP, of which 12.5 percent STOIIP had been produced by the end of 1963 from some 540 wells. By then the average pressure had decreased from an initial value of 960 to 550

Jan 1, 1970

By T. K. Perkins, L. R. Kern, R. E. Wyant

Many factors influence the results of hydraulic fracturing to stimulate well productivity. Most of these factors have been studied and their effects discussed in the literature. However, the movement of sand through the crack during fracturing has received little consideration. Exactly where the sand goes in the fracture, how it gets there, and how this might influence results have not been studied, judging from published reports. To get an answer to some of these questions we have studied sand movement through a simulated crack segment in the laboratory. Although the work is preliminary, we have answered some questions about sand movement in a semiquantitative sense. We are reporting our results at this time in the hope of stimulating interest and further study in this important aspect of fracturing. LABORATORY EXPERIMENTS Fig. 1 is a schematic diagram of equipment used in this study. A segment of a vertical crack 0.25-in. wide was formed between two Plexiglas plates, and movement of sand in fluid through the crack was observed under various conditions. Sand and fluid were injected into the crack at constant rates. The top and bottom of the crack were closed so that the fluid flowed horizontally. Because of gravity most of the sand settled to the bottom before moving very far through the crack. If the fluid velocity exceeded a certain critical value, then all the sand injected was washed on through the crack even though it had settled to the bottom. When the fluid velocity was less than this value when sand injection was started, most of the sand which settled in the crack remained in place until enough of the area open to fluid flow was blocked by the settled sand to increase the fluid velocity to the critical value.. At this point all injected sand moved on through the crack even though it settled to the lower portion of the moving fluid section. All sand previously settled remained in place. When fluid velocity was increased at this paint by increasing fluid injection rate (without also increasing sand injection rate), some of the previously settled sand was washed out until enough area was open to flow that the velocity dropped to the critical value. Conversely, when the fluid velocity was decreased by decreasing the fluid injection rate (without changing the rate of sand injection), additional sand would settle until the critical velocity was again attained. In short, the critical velocity is an equilibrium value which the system will automatically seek, if possible. The equilibrium (or critical) velocity was measured under several conditions. The results of these measurements are shown in Pig. 2 where the equilibrium velocity is shown as a function of the amount of sand injected per minute. The points falling in the shaded band shown in Fig. 2 were obtained with water. The higher viscosities were obtained by adding CMC (carboxymethylcellu-lose). Data indicate that the equilibrium velocity is not very sensitive to sand size or transport fluid viscosity, or to sand injection rate (except at very low rates). Points connected by the solid CUrve in Fig. 2 were obtained with water and 20- to 40-mesh steel shot. The equilibrium velocity is apparently fairly sensitive to difference in density between particle and fluid. Only one run was made with a gelled fluid and that was gelled water at zero sand injection rate as shown in Fig. 2. (Measurements of equilibrium velocity at zero sand injection rate were made by filling the crack

By H. T. Kennedy, S. G. Dardaganian, J. L. Clanton

This paper presents equilibrium ratios measured on oil and gas samples taken from the Kelly-Snyder field at the reservoir temperature of 128°F. From these data smoothed curves were constructed for the reservoir fluid up to the bubble-point, 1,874 psia. A material balance calculation, based on the use of equilibrium ratios instead of the usual solubility and shrinkage meayr~rements is made on the field and compared to the con\.entional method. DETERMINATION OF EQUILIBRIUM RATIOS The equilibrium ratios, or K values were obtained by recombining separator samples of oil and gas taken when the reservoir was above its bubble point, and thus represented the reservoir fluid before any change had occurred due to the evolution of gas. The well sampled had a gas-oil ratio of 617 cu ft of high pressure sepalator gas per barrel of separator oil. and a bottom-hole temperature of 128°F. Samples were brought to equilibrium by agitating for 24 hours in a thermostat, maintained at the proper temperature within 0.2ºF, mercury being added when necessary to maintain the pressure constant. Samples of oil and gas were then taken by constant pressure displacement, and analyzed by low temperature frac-tionation. Equilibrium ratios were determined at 128°F and 500, 800, 1,000 and 1,400 psig. At the latter pressure, there were indications that the gas samples were contaminated with small quantities of oil, and only data taken at the first three pressures were included in Fig. I. This figure has been extended to higher and lower pressures by the use of the Fluor charts for 3,000 psi convergence pressure, on which the experimental points fit well. MATERIAL BALANCE EQUATIONS The possible errors involved in conventional material balance equations have been analyzed by Muskat1. Where the errors occur in field measurements, or in obtaining effective averages from variables such as gas-oil ratio and pressure, they will evidently be involved in any attempt to estimate the content of the reservoir. On the other hand, substantial errors may be introduced in conventional material balance calculations by: (1) the neglect of the amount of oil dissolved in gas at high pressures and temperatures in the reservoir, but which condenses under separator conditions; (2) the assumption that the separation of oil and gas in the reservoir may be duplicated by a differential liberation in the

By M. E. Hassan, J. C. Calhoun, R. F. Nielson

An apparatus was constructed for the measurement of inter-facial tensions over a range of temperatures and pressures. This apparatus utilized the pendent drop method, and resembles in construction similar apparatus recently described in the literature and in use in some petroleum research laboratories. The interfacial tensions of benzene, propane, n-pentane, n-hexane. n-octane. and iso-octane against water were measured at temperatures ranging from 26" to 82" C and at pressures ranging from 1 to 204 atm. Values of interfacial tensions for the henzene-water system and their variations with temperature and pressure are generally in good agreement with value? of previous investigations. The data in all cases showed a slight decrease of interfacial tension with pressure at constant temperature in the range studied. The effect of pressure became less as the pressure was increased, with an indication of a reversal of the effect at higher pressures. There was a decrease of interfacial tension with temperature at constant pressure in all cases, as would normally be expected. This rate of decrease became greater the higher the temperature. A general equation is presented for the interfacial tensions as a function of presPurp and temperature over the range studied, and the constants calculated for each system. A definite trend was found in the effect of molecular weight on the interfacial tension at a given temperature and pressure, for the homologous series from propane to n-octane. Data for n-decane from the literature fitted well into this trend. INTRODUCTION It has been recognized for many years that surface forces play an active part in the production history of an oil reservoir as well as in determining the amount of unrecoverable oil. The magnitudes of these forces are governed by the values of the interfacial tensions. These vary greatly with composition, pressure, and temperature. The effect of composition is greater than that of pressure and temperature over the ranges normally encountered. In fact, because of the variety of compounds present in a crude, no attempt has yet been made to correlate crude composition with its surface tension or with its interfacial tension against water. The amount of dissolved methane greatly affects the surface properties of an oil, and there are "surface-active" compounds in crudes which have a much greater effect on boundary tensions than do equal amounts of hydrocarbons. The purpose of the present work was partly to contribute data that would eventually be used in correlations applicable

Jan 1, 1953

By A. S. Odeh

During the last few years, several authors1-6 have advanced mathematical solutions, both exact and approximate, to the pressure behavior of a well producing from a region bounded by a circular discontinuity. The most recent solutions were reported by Carter, and by Bixel and van Poollen. Varter solved the appropriate equations analytically. Bixel and van Poollen reported a numerical solution and provided correlation charts for the interpretation of flow tests; their method uses an overlay technique and is valid for wide ranges of mobility and storage capacity ratios. The mobility ratio is defined as the ratio between k/p of the two zones, and the capacity ratio is defined as the ratio between $c of the two zones. We have developed independently an analytic solution, programmed for the CDC 1604, that agreed with Carter's.'1 The same differential equation, assumptions, boundary, and initial conditions were used. We ran 30 drawdowns and 30 buildups of a unit capacity ratio and of varied mobility ratios. We also varied the radial distances to the discontinuity, and the radii of drainage. We examined the results and obtained correlations that were simple to apply and useful in flow test interpretation. The results and the correlations reported here in graphical form allow the determination of the mobility of the two zones and the radial distance to the discontinuity. They eliminate the need to fit field data to a curve as is required in the overlay technique.6 Thus, the advantages of the present correlations are that they require less effort and they provide a unique answer. They also provide guidelines for selecting the parts of the flow test plot to be used in the interpretations. Moreover, the results obtained from the programmed analytic solution define the time intervals where the approximate solution reported by Hurstl and discussed by Larkin2 are applicable. Normal buildup and drawdown plots resulted in an early straight line, a long transition zone and a late (second) straight line. However, the second straight line will not occur if the drainage radius is not at least 10 times the radial distance to the discontinuity. The mobility of the imer and outer zones can be calculated from the slopes of the first and second straight lines, respectively, by the normal methods of the flow test analysis. We found that the radial distance to the discontinuity may be determined from the abscissa value for the point of intersection of the extensions of the two straight lines. Results similar to some of these were reported by Adams et aL7 The approximate solution compared favorably with the exact one during the first and second straight-line periods. However, during the transition zone, deviation occurred between the two solutions. Discussion of Figures Figs. 1 through 7 summarize the results of this investigation. On Figs. 3 through 7 there appears the equation of each respective plot. In this paper the mobility ratio M is defined as klp2/kzpl, where subscripts 1 and 2 refer respectively to the inner and outer zones, and dimensionless time = t, = 6.32 X In the investigation, M ranged from 0.25 to 50;

Jan 1, 1970

By R. F. Jones, N. E. Truitt, M. Prats

In most published discussions and theories of in situ combustion, the combustion fronts are assumed to be vertical. However, evidence from field tests leaves no doubt that combustion fronts often advance more rapidly along the top than near the bottom of a formation as a result of difference in density between injected air and formation liquids. The approximation proposed in this paper to determine the movement of the resultant tilted combustion surfaces states that the vertical rate of movement of combustion surfaces is proportional to the horizontal oxygen flux. Where self-ignition is possible, the proposed method demands that a secondary combustion surface exist around production wells which produce some oxygen. These secondary combustion surfaces may be formed long before the primary combustion surface can advance from injection to production wells. Heat liberated near production wells at these secondary combustion surfaces can contribute to an early increase in production rate. Results indicate that significant oil recoveries cannot be obtained from the usual flood patterns (five-spots, seven-spots, etc.) without producing large volumes of unused oxygen. Ideally, to increase oxygen- consumption efficiency, well patterns should allow oil production from a first line of production wells and gas production from more distant lines of producers. However, it may be desirable to produce some gas at all wells to support (and benefit from) active secondary combustion surfaces. Results indicate that the well spacing through which combustion can be advanced is larger than that predicted by other methods. A large number of production wells may still be desirable to take quick advantage of gravity drainage. From a comparison with results at South Belridge field, California, it appears that this method adequately describes oxygen concentration and temperature histories and combustion-front shapes. However, this method does not accurately locate the most advanced point of the combustion surface. There is some field evidence to substantiate the actual presence of secondary combustion surfaces at South Belridge. Use of the proposed method appears warranted at this time when lay-over of the combustion surface can be anticipated. INTRODUCTION The assumption of vertical combustion fronts has been embodied in all previous publicationsl-6 which use the movement of combustion fronts away from injection wells to determine the temperature and fluid distributions in the reservoir. The only paper concerned with a mathematical model of the combustion process in which a nonvertical combustion front is used was written by Gottfried.11 Actually, nonvertical combustion fronts have been observed in most in situ combustion field tests for which adequate data are available.7, 8, 6 In practice, the typical vertical extent of the burned zone decreases with distance from the injection well, and this burned zone is at or near the top of the sand body. In some cases, such as at South Belridge field near Taft, Calif., the combustion surface is almost horizontal over a very large area.8 Thus, for some years an obvious and serious gap has existed between theory (vertical fronts) and practice (tilted fronts). This is indicated in Fig. 1. Tilted combustion fronts such as observed at South Belridge8 sometimes result from the natural tendency of injected gases to rise to the top of an oil sand. An initially present gas cap, even if thin, would make apparent the effect of gravity segregation from the first, particularly since formations considered for in situ combustion field operations are generally selected for their

Jan 1, 1969

By L. F. Elkins

In 1946 Littlefield, Gray and Godbold published a thorough geologic description of the West Edmond Hunton Lime reservoir, located in Central Oklahoma, and discussion of its early performance.l They stressed implications of the nature of the reservoir rock on performance to be expected with production both by natural means and by pressure maintenance with gas injection. The essence of their analysis was that presence of an extensive interconnected system of fractures and solution channels constituting about 10 percent of the total reservoir void space in otherwise very low permeability reservoir rock would result in severe channeling of naturally encroaching water or injected gas through the fractures with little or no benefit to ultimate recovery of oil. Considerations were in progress to unitize the reservoir for pressure maintenance by gas injection. Some other engineers involved in the project had considerably different opinions at the time regarding continuity of fractures and regarding effects of such fractures on performance. The field was unitized Oct. 1, 1947. Subsequently, pilot gas injection tests were conducted in four areas of the field; large scale waterflooding was conducted in two parts of the field in the Bois d'Arc or Upper Hunton; and pilot waterflooding was conducted in the Chimney Hill or Lower Hunton. Full scale water-flooding is now in progress in the Chimney Hill section where it is oil productive. The primary purpose of this paper is to divine some of the secrets of this reservoir through analysis of all aspects of its performance, with the hope that it will be of help to engineers and operators in developing and operating other low permeability fractured reservoirs. In part, this considers variance of actual performance from some widely used conventional reservoir engineering assumptions. For brevity, other aspects of unit operation and its benefits to the owners are not discussed. Oil and Gas in Place For unit participation, net pay in the Bois #Arc defined by the gross SP anomaly on electric logs totalled 2,100,000 acre-ft in the 30,000-acre reservoir. Development of West Edmond predated modem logging techniques that now would be used to determine "net" or "productive" pay. In Fig. 1 this log pick is compared with detailed core analysis and core descripfion by Littlefield et a1.l for one well identified as A on the index map in Fig. 2. Not all of the gross Bois #Arc section had oil staining, but all cores with measureable intergranular permeability were stained. (The committee selection of the Bois d'Arc SP interval also includes the Frisco section as identified from cores by Littlefield.) The Bois d'Arc is oil productive throughout the entire 600-ft monoclinal accumulation from the eastern upstructure erosional pinchout to the west-side oil-water contact. (See Littlefield et al.l and McGee and Jenkins for geologic detail.)

Jan 1, 1970

By E. L. Anders

The Mile Six pool is located on the La Brea-Parinas Cullcession of International Petroleum Co., Ltd., in northwestern Peru on the west coast of South America. The reservoir pressure in this pool has been maintained within 200 psi of its initial value throughout its history, and gravity drainage has played an important role in the production behavior. It has now produced 95 per cent of its estimated ultimate recovery. It is estimated that this interesting oil pool will ultimately produce 67 per cent of the initial oil in place and that the resulting residual oil saturation may be as low as 19 per cent of the pore volume (29 per cent of the hydrocarbon pore volume). An evaluation of reservoir rock and fluid characteristics and ultimate oil recovery is presented. INTRODUCTION This study of Mile Six pool was made to evaluate its performance according to latest available information. The production performance of this pool has been discussed in various articles in the past. and the reported behavior has been used as an example for application of computation procedures for gravity drainage depletion' and as an illustration of field behavior under gravity drainage or expanding gas cap drive. There have been wide variations in reported values of initial oil in place, reservoir oil volume factor. connate-water sauration, volume of effective sand, and ultimate recovery because of the paucity of reliable basic data. These various factors have been determined as accurately as practicable with the latest available information, and this evaluation is presented herein. The production history of Mile Six is an excellent example of gravity drainage depleion with effective pressure maintenance by gar injection. GENERAL Mile Six pool was discovered by cable-tool drilling in November. 1927. when well 1996 was completed in the Parinas sand. After slow development with cable tools and sporadic production. the pool was opened to continuous pro(iuction in November. 1933. and develpment was completed with rotary rigs. Pressure maintenance was started in December, 1933. by returning gas to upstructure wells. Most of the development was 'completed by 1937, but some additional wells were drilled in the period 1939-1947. and several old wells were deepened. A total of 46 oil and gas wells and 4 dry holes were drilled on approximately 7-acre spacing. Of the producers. 21 are now flowing. 2 are pumping. 44 are gas input wells. 3 are abandoned. I is a gas well shut in. and 15 are shut in because of non-commercial production or high gas-oil ratio. The locations of all wells are shown on the map of Fig. I. Total oil production on Dec. 31. 1952, was 30,867,373 bbl: cumulative gas production was 22,023,777 Mcf; and 26,410,946 Mcf of gas had been returned to the reservoir. These figures do not include oil and gas lost ill a blowout in January. 1940. GEOLOGICAL DESCRIPTION Mile Six pool is located on the northern end of a structural spur projecting from the La Brea-Negritos uplift.' The spur is probably a reflection of a basement structure. It plunges gently to the north, i broken into a complex series of fault blocks. and contain.; the Verdun Alto. Section Sixteen. and Mile Six pools. The Parinas handstone (lower Eocene!. which is the producing formation in Mile Six. occurs at an average depth of 2,200 ft in the pool and dips north and east at from 15' to 20". The pool covers an area of approximately 350 acres. Mile Six is downfaulted about 600 ft from Section Sixteen pool to the soutb. and a major fault forms its western boundarv. The north and east boundaries are formed by the intersection of the sand top with the water-oil contact which occurs at approximately 2.440 ft subsea. An original gas-oil contact probably existed at about 1.875 ft subsea. Fig. 1 presents the latest structural interpretation of the pool. and Fig. 2 is an isopach map showing thickness of the total Parinas formation above the original water-oil contact. The heavy lints of Fig. 1 are contours on the sand top, and the fine lines are contours on the fault planes. This type of straight-line structural map was developed 1)) International's geologists to reflect structural conditions where the bedding planes dip and have no curvature. 'The La Brea-Parinas Concession is highly faulted by normal fault.. The beds are flat wherever exposed. The Parinas formation is approximately 635 ft thick. and it is etimated that 62 percent II the formation is effective sand. The original oil zone was about 565 ft thick. Fig. 211 presents an electric log showing typical Parinas sand development ill Mile Six pool. The Parinas band in Mile Six i. a well-sorted. medium- to-coal-se-grained. cross-bedded sand with minor lenses of shale and small lenses and pockets of pebble conglomerate. The sand grains are subangular to rounded and consist chiefly of quartz with feldspars. biotite hornblende, and augite as accessory minerals. Because of faulting of the Parinas- formation to the east and north of the pool. there is probably l possibility of a significant, natural water drive in Mile six. The faults within the pool. as indicated in Figs. 1 and 2. are of smaller disI,laceInent and seem to act only , partial barriers to fluid movement within the reservoir. RESERVOIR CHARACTERISTICS Core analysis data are available from five wells. The data were obtained from three wells (Nos. 3401. 3586 and 3719) at the time of their completion and from well 1996 when the original liner was sidetracked and the well was deepened ill 1946. Data from well 2779 were obtained in 1943 from old cores taken when the well was deepened in 1934. From these core analyses. the average porosity was estimated to be 22.6 per cent. and the average permeability to dry ail. was estimated to be 780 and Measured productility indices varied from 3.1 to 71.4 B/D per psi differential. Specific productivity indices varied approximately from 0.1 to 0.3 B./D per psi per ft of sand.

Jan 1, 1953

By H. S. Price, J. C. Cavendish, R. S. Varga

A numerical formulation of bigh-order accuracy, based on variational methods, is proposed for the solution of multidimensional diffusion-convection-type equations. Accurate solutions are obtained without the difficulties that standard finite difference approximations present. In addition, tests show that accurate solutions of a one-dimensional problem can be obtained in the neighborhood of a sharp front without the need for a large number of calculations for the entire region of interest. Results using these variational methods are compared with several standard finite difference approximations and with a technique based on the method of characteristics. The variational methods are shown to yield higher accuracies in less computer time. Finally, it is indicated how one can use these attractive features of the variational methods for solving miscible displacement problems in two dimensions. INTRODUCTION The problem of finding suitable numerical approximations for equations describing the transport of heat or mass by diffusion and convection simultaneously has been of interest for some time. Equations of this type, which will be called diffusionconvection equations, arise in describing many diverse physical processes. Of particular interest here is the equation describing the process by which one miscible liquid displaces another liquid in a one-dimensional porous medium. The behavior of such a system is described by the following parabolic partial differential equation:

Jan 1, 1969

By R. W. Snyder, W. H. Guilinger, B. S. Gottfried

Two numerical methods are presented JOT solving the equations Jor one-dimensional, multiphase /low in porous media. The case oF variable physical properties is included in the Formulation, although gravity and capillarity are ignored. Both methods are analyzed mathematically, resulting in upper and lower bounds for the ratio of time step to mesh spacing. The methods are applied to two- and three-phase waterflooding problems in laboratory-size cores, and resulting saturation and pressure distributions and production histories are presented graphically. Results oJ the tuo-phase flow problem are in agreement with the predictions of the Buckley-Leuerett theory. Several three-phase flow problems are presented which consider variations in the water injection rate and changes in the initial oil-and ujater-saturation distributions. The results are different physically from the two- phase case; however, it is shown that the Buckley-Lez~erett theory can accurately predict fluid interlace velocities and displacing-fluid frontal snturations for three-phase flow, providing the correct assumptions are made. The above solutions are used as a basis JOT et~aluuting the numerical methods with respect to machine time requirements and allowable time step for a fixed mesh spacing. INTRODUCTION Considerable progress has been made in recent years in obtaining numerical solutions of the equations for two-phase flow in porous media. Douglas, Blair and Wagner2and McEwenll present different methods for solving the one-dimensional case for incompressible fluids with capillarity (the former using finite differences, the latter with an approach based upon characteristics). Fayers and Sheldon4 and Hovanesian and Fayers8 have extended these studies to include the effects of gravity. West, Garvin and Sheldon, l4 in a pioneer paper, treat linear and radial systems with both capillarity and gravity and they also include the effects of compressibility. Douglas, Peaceman and Rachford3 consider two-dimensional, two-phase, in compressible flow with gravity and capillarity and Blair and Peaceman1 have extended this method to allow for compressible fluids. No one, however, has examined the case of three-phase flow, even for the relatively simple case of one-dimensional flow of incompressible fluids in the absence of gravity and capillarity. In obtaining a numerical technique for simulating forward in situ combustion laboratory experiments, Gottfried5 has developed a method for solving the one-dimensional, compressible flow equations with any number of flowing phases. Gravity and capillarity are not included in the formulation. The method has been used successfully, however, for two- and three-phase problems in a variable-temperature field with sources and sinks. This paper examines the algorithm of Gottfried more critically. Two numerical methods are presented for solving the one-dimensional, multi-phase flow equations with variable physical properties. Both methods are analyzed mathematically, and are used to simulate two- and three-phase waterflooding problems. The numerical solutions are then taken as a basis for comparing the utility of the methods.

Jan 1, 1967

By J. C. Martin

Equations for three-phase, three-dimensional, compressible flow (including capillarity) are reduced to two-dimensional relations by a partial integration. This reduction allows three-dimensional flow problems to be treated with mathematics for only two spatial dimensions. The results can be used to formulate flow equations for two-dimensional reservoir simulators in which the effects of capillarity and fluid segregation in the third dimension are represented. Such reservoir simulators would retain many of the aduantages of two-dimensional simulators while simulating three-dimensional effects. The principal restriction of the method is that the thickness of the reservoir should be small, compared to the distance across the reservoir. INTRODUCTION In recent years, computers have been used to calculate the performances of many reservoirs. Most of the detailed calculations, however, are based on finite difference solutions of the flow equations, and present day computers are seldom able to handle a sufficient number of cells to produce entirely satisfactory solutions, even for reservoirs represented by two-dimensional arrays of cells. The situation becomes much worse when one wishes to approximate the reservoir by a three-dimensional array. A great saving in computation or a more detailed solution can be obtained for many reservoirs by using the partial integration of the equations of flow presented in this paper. The integration reduces the three-dimensional equations to two-dimensional relations; and, for studies of two-dimensional flows in vertical cross-sections, the equations are reduced to one-dimensional relations. Most reservoir performance calculations currently are based on one- or two-dimensional flow re1ations.l-3 In some cases flow in the third dimension is approximated by assuming a particular type of vertical saturation distribution, such as gravity segregation.l.2. The relations developed in this paper approach those for segregated flow as the capillary pressures approach zero, and they approach those for uniformly distributed saturations as the capillary pressures are increased. For this analysis, the ratio of the reservoir's thickness to the maximum distance across it must be small. The capillary pressures between the oil and water should also be small compared to the maximum pressure difference in the reservoir. This requirement is met by most reservoirs. It is assumed that the capillary-pressure curves are well defined, whether or not hysteresis effects are included. Also, the reservoir must have sufficient vertical permeability to allow the fluids to segregate. The results presented here provide a firm theoretical foundation to Coats' et aL5 assumption of vertical equilibrium and extend the relations to three-phase flow. Coats' assumption of vertical equilibrium, which he verified by calculations and experiment, is developed here mathematically from basic flow equations. DISCUSSION SATURATION AND PRESSURE DISTRIBUTIONS Appendix A presents a mathematical analysis of fluid flow in reservoirs where the ratio of thickness to maximum distance across the reservoir is small. The results indicate: (1) that the fluids along any line perpendicular to such a reservoir's upper surface are in static capillary equilibrium (vertical equilibrium); (2) that, to a first approximation, the fluid pressures and properties are functions of only areal position in the reservoir and time; and (3) that hydrostatic pressure gradients exist along any line perpendicular to the reservoir's upper surface. The results might be expected after studying several physical considerations. First, no flow is allowed normal to the upper and lower reservoir boundaries, which are relatively close together.

Jan 1, 1969

By H. J. de Haan, L. Schenk

Scope for Thermal Recovery in Shell's Heavy Oil Fields in Venezuela The main heavy oil reservoirs on the East coast of Lake Maracaibo (Fig. I), known as "Bolivar Coast", initially contained some 20 billion bbl of oil in place with gravities in the range of 10 to 15° API. The current total production rate is about 400,000 B/D. Since the recovery to date is on an average only 12.5 percent of the initial oil in place, these fields offer a vast potential for secondary recovery methods. The reservoirs are characterized by moderate depth (generally 1,000 to 3,000 ft; maximum 5,000 ft), good formation properties such as net oil sand thickness (50 to 300 ft), high porosity (30 to 40 percent) permeability (1 to 3 darcies) and oil saturation (initially about 80 percent) and high oil viscosity (100 to 10,000 cp in situ), all of which are favorable for the application of thermal recovery processes. Since 1957, steam drive, steam "soak" and underground combustion have been tested or are being tested in this area.1-3 This paper deals with a major steam drive test in the Tia Juana field. The test was commenced in Sept., 1961, on the basis of encouraging results of laboratory experiments and pilot field tests, which will be discussed briefly. Early Laboratory Investigation of the Steam Drive Process In 1956 a series of model experiments was carried out in the Koninklijke/Shell Laboratorium, Amster- dam, to investigate the displacement of heavy oil by steam.' The prototype studied was a horizontal, un-consolidated heavy oil reservoir, subjected to a linear steam drive. The model, representing a slab from this reservoir along the main flow direction, consisted of a steel tube filled with oil-saturated sand, provided with layers of rock at the top and bottom to simulate conductive heat losses to cap and base rock. Lateral heat losses, perpendicular to the main flow direction, were minimized by insulating the sides of the model. Heat flow and viscosity and gravity forces were scaled by taking the steam injection rate and sand permeability inversely proportional to the geometric reduction factor of the model, and using materials (oil, steam, cap rock, etc.) identical with those of the prototype. The main feature of the process was a frontal displacement mechanism that was observed within wide ranges of injection pressure and rate, oil viscosity, initial oil saturation and sand permeability. The high stability of this front can be attributed to the high volumetric steam flow rates, combined with the reduction in oil viscosity. High recoveries were obtained due to the low residual oil saturation of about 15 percent in the zone swept by steam. The latter value pertains to the specific test conditions (heavy oil, uncon-solidated sand) and may be lower for light oil, due to the effects of steam distillation. The applicability of this experimental approach is limited, as it does not take into account the fact that

Jan 1, 1970

By R. W. Parsons

The over-all apparent single-phase permeability of fracture-rock systems was studied using two different two-dimensional models. In a strict sense the results are applicable only to these models, yet may give some insight into real reseruoir behavior. The regular fracture-matrix model can have any number of regular fracture sets. All fractures within each set have the same orientation, width and spacing. The gross behavior is amlytically equivalent to that of an anisotropic permeable medium. In the heterogeneous fracture systems, individual fracture conductivities vary within a two-dinensioml network pattern. A computer program calculates the overall permeability after the conductivities are randomly placed in the network. The calculated permeability depends upon the distribution of these conductivities and the particular fracture pattern into which they are placed. Although detailed information on actual reseruoir fracture systems is not available, this study shows how some of the possible variables may affect obserued gross reservoir characteristics. INTRODUCTION Characteristics of reservoir rocks are determined by direct (core analysis) or indirect (logging and production tests) examination. Both methods, singly and together, are widely used in reservoir evaluations. For fractured reservoirs a problem exists in relating characteristics of cores to the in situ reservoir properties. Knowing fracture characteristics can be of value in determining production mechanisms. Because of the limitations of observation, the intimate details of a fractured rock system will never be known. One approach in analyzing fractures is to study the behavior of some conceptually simple models. The two models analyzed contain features which may be realistic, yet are mathematically tractable. Flow is considered in two dimensions only, with the plane of interest at right angles with the fractures. This concept represents the horizontal flow in a reservoir system with the fractures oriented vertically. Only single-phase, laminar flow is considered. The regular fracture-rock system permeability is calculable from a simple analytical expression. The heterogeneous system requires a computer program for calculating overall permeability.

Jan 1, 1967