Search Documents

By C. S. Matthews, M. Prats, R. I. Jewett, J. D. Baker

By mathematical analysis it was found that injectivity history of a uniform five-spot pattern can be calculated by rather simple formulas. These calculated injectivities were found to agree rather well with injec-tivities measured in a potentiometric analog. With this as a basis, a simple method was finally developed which allows prediction not only of rate of injection but also of rate and kind of production at the production well, for a uniform five-spot. Using this method formulas for the eflective injectivity and production behavior of a waterflooded reservoir having a wide range in permeability can be calculated by considering that the reservoir consists of several layers, each having a uniform but different permeability. Rather good agreement is shown by a comparison of the observed production history of a field in the Illinois basin and that calculated by the methods just described. INTRODUCTION The importance of being able to predict the injection rates for a water flood is well known, for upon injection rate depends the life of the flood, the size of pumping and treating facilities, and the rate of oil recovery. One method of determining rates of injection is through a pilot flood. However, it appears that even if a pilot flood is resorted to, a theoretical method for calculating injection rate will be valuable in extending field results to times after inter- ference and to locations of other five-spots. It is also important to be able to predict the production history of a flood. Such prediction must generally be done by theoretical means, since in general a small pilot flood will not furnish much quantitative information in this regard because of the distortion after oil-bank interference. Methods for predicting the behavior of five-spot water floods have been proposed by Yuster and Cal-houn' and by Hurst.? However, these methods do not consider the effect of the mobility of the different fluids in the reservoir. Other investigators3,4,5,6 have determined the effect of the water-to-oil mobility ratio on the production history by means of different experimental techniques. Most of the published work applies to homogeneous sand bodies, does not provide information for determining the injection requirements of a flood, and seldom considers the rate of build-up of the oil bank as it fills up the partially depleted reservoir.* The present work was undertaken because of the lack of predictive method which takes into account mobility ratio and the rate of buildup of oil bank resulting from the void space in a partially depleted reservoir. As water is injected into the reservoir, which contains oil, water, and gas in macroscopically homogeneous saturation, it forms several banks ahead of the injection well. In each bank or region there is assumed to be only one mobile phase—water in the water bank, oil in the oil bank and gas in the gas region. The saturation in each region and a plan view of the banks are shown in Figs. 1 and 4, respectively. FACTORS AFFECTING INJECTION AND PRODUCTION RATES Results in this report are determined as functions of the three parameters, F, M, and M, , and of a fourth parameter, r,/L** which relates the size of the well to the flood spacing. The parameter F is the displacement factor defined by Hurst in Ref. 2. It determines the rate of build-up of oil bank as it fills up the partially depleted reservoir. When gravity effects are neglected, the four parameters just given are sufficient to describe the waterflood behavior for horizontal and homogeneous reservoirs having incompressible liquids. Definition of InjectIvity The dinlensionless injectivity (or injection rate) used in this report is defined by

By R. J. Wygal, J. Naar

An equation for three-phase (water, oil, gas) imbibition oil permeability is developed, assuming the water to be the dominant wetting fluid. Oil isoperms are obtained for consolidated sandstones characterized by 1/pc2 = CS*. The evolution of an oil-gas system imbibing water from Swi to(Sw,imb )max is shown to proceed along a line of constant oil saturation with increasing oil permeability and decreasing gas saturation. When the gas satuation cannot be reduced further, the system evolves along a line of constant Sg with decreasing oil saturation and permeability. The initial gas saturation is shown to reduce markedly the effect of complete wetting by either oil or water on flow performance. INTRODUCTION Imbibition oil isoperms are required for performance prediction when a well is producing water, oil and gas. This situation occurs in multiphase displacements such as underground combustion, steam injection and the water flooding of highly depleted reservoirs. In a recent paper,1 a model was presented for the prediction of two-phase imbibition characteristics. This paper extends the imbibition model to the case of three phases by assuming that the water is the dominant wetting fluid. The following results were obtained from the model: (1) an analytical expression of oil isoperms; (2) oil isoperrns as functions of reduced water, oil and gas saturations, valid for all sandstones having a capillary pressure curve which can be approximated by 1/PC = CS*; and (3) evaluation of the three-phase flow performance as dictated by complete wetting by either oil or water. The agreement between predicted and observed oil recovery in the presence of a gas phase, reported in Ref. 1, is a partial support for the present development. However, experimental data are not available at this time to check fully the model predictions. Perhaps this paper will stimulate the collection of such data. The imbibition model of a porous medium has been described previously, and the reader is referred to the paper of Naar and Henderson1 for details. In brief, the model is formed by the random interconnection of straight capillaries, with a provision for the blocking of the non-wetting phase by the invading wetting fluid. The following expressions were derived from the model. 1. The amount of oil blocked by the invading water is

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 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 W. F. Hower, W. Brown

Large-scale sand consolidation tests were conducted in an effort to determine the reasons for the successes and failures of this method of sand control. Several different consolidating materials were used in treating both clean and bentonitic sands that were packed in a chamber having a capacity of 3.3 cu ft. The results were essentially the same for all of the different consolidating materials, The data show that low-viscosity consolidating materials pumped at a relatively slow rate gave the best results. Where the formation has produced sand, the treating fluids can compress the formation, thus permitting the channeling of fluids to another horizon. Pressure-packing these zones before attempting to consolidate is recommended. Sands containing more than 4 per cent of water-swelling clays are not good candidates for consolidation. It is indicated that loose sand, particularly when it is bentonitic, can be fractured during the placement of the treating fluids. INTRODUCTION Sand production in oil and gas wells has plagued the industry for many years, and numerous cures for this problem have been suggested. Most methods have been successful to a certain degree, but the great variety of well conditions that exist in the different areas has magni- fied the problem and limited the successful use of the various systems. Four review papers1-4 present a wealth of information concerning the degrees of success that have been obtained by the different sand-control methods. The bridging of sand grains by the use of gravel packs and screens has been quite successful. However, these methods do not leave the casing clear for all types of multiple completions, and the cure does not last for the production life of the well in some instance:;. The control of loose sands by sand consolidation with resins has never been as successful as desired. It has always been hoped that such a treatment would eliminate all sand problems for the life of the well, but. initial applications, starting in the middle 1940's, were only moderately successful. Lott, et a1,3 reported a success ratio of approximately 50 per cent and made the following conclusions. The highest percentage of successes were obtained where: a. Consolidation of a zone was made at the time of initial completion or prior to the production of sand. b. The interval treated was less than 12 ft in length. c. Between 30 and 50 gill plastic/ft of producing interval was displaced through the perforations. REASONS FOR SAND CONSOLIDATION FAILURES Our own experiences in the field of sand consolidation point toward the following conditions as the major reasons for the failure of sand consolidation attempts. 1. Mud-plugged perforations and mud invasion of the formation. 2. Sand in the casing covering all or part of the perforations. This sand could be either formation sand or one of the coarser sands used as propping agents in hydraulic fracturing. 3. Holes in the casing. 4. Channels behind the casing. 5. Attempting to treat too long a perforated section. 6. Too high a percentage of water-swelling clays in the formation. 7. Formations that have produced sand. Recent attempts were made to treat perforated sections ranging from 10 to 30 ft, in wells that have produced sand, by using a straddle packer that was raised and lowered through the perforations as the consolidating material was being pumped. In most instances, the pressure required to pump fluid into the formation varied considerably as the tool was raised and lowered. This suggested the possibility that significant differences in permeability were present or that only part of the formation had produced sand. There were times when a sudden break in pressure indicated that a fracture was being formed. Research conducted several years ago concerning the problem of the control of water in air and gas drilling indicated that shale sections could be fractured quite easily. In addition, it was determined that it was easier to pump fluids into shale bodies by fracturing the shale itself, or the interface between the shale and sand, than to pump into a fluid-saturated formation. Formations that produce sand are usually adjacent to shale bodies and frequently have shale streaks of various thicknesses inter-bedded in the sand. Therefore, where shale is exposed to fluid pressure it

By R. V. Higgins, A. J. Leighton

Interest by petroleum engineers in the flow of three phases—oil, gas and water—in irregularly bounded porous media lies mostly in the performance calculation of water floods of reservoirs that have been partially depleted as the result of expansion of much of the originally dissolved gas. The authors present a method to forecast three-phase flow in complex geometry and explain the details by the use of a specific example of a five-spot water flood of a partially depleted stratified reservoir. In this example, the fluid and rock properties of a field given by Prats, et al,' were used. The computed results fit the field performance more accurately than Prats, et al, and Slider.' The time required for the high-speed digital computer to make the calculations, including the contributions from the different layered zones, is about one minute. INTRODUCTION The declining rate of discoveries of new oil fields in the United States makes the recovery of more oil from known reservoirs more attractive than previously. Water flooding has an excellent proved background for recovering additional oil economically. Accordingly, the main interest of this paper is in this type of recovery. In petroleum engineering studies, commercial interest in three-phase flow is mostly in the water flooding of reservoirs in which the oil has been partially produced by the expansion of dissolved gas. When water is pumped into these reservoirs, three-phase flow takes place. Although this paper is concerned with these conditions, the principles involved could be used for other conditions when and where they come to the fore. Many investigators have made contributions to non-empirical forecast methods using basic scientific engineering principles. Several of these used the oil in place at the start of the water flood and the oil remaining after a large quantity of water has passed through a core for key values in their calculations. Recently, Prats, et al,' and Slider have reduced assumptions by adding the third phase—gas—in their forecasting methods. Slider uses the mobilities in the immediate vicinity of inlet and outlet wells as an aid to simulate the resistance effect in the five-spat pattern of the flow of oil, gas and water. Prats, et al. minimize assumptions by using previously determined laboratory data for determining sweep efficiency in the five-spot pattern. In neither of these papers is the saturation profile continuously affected by permeability-saturation curves. Sheldon and Dougherty- recently described a method that employs continuously changing saturation profiles using permeability curves, has a minimum of assumptions and needs no prior sweep efficiency. The Higgins-Leighton method, described in this paper, has all of these features; however, many of the techniques are different from those of Sheldon and Dougherty. The Higgins-Leighton method, tested in May, 1961, is direct and easy to apply and requires very little computer time to calculate a forecast. The short computer time is especially helpful in the study of a reservoir containing many layers of different relative permeabilities. In the Higgins-Leighton method, the individual pressures do not have to be calculated, as the resistance to flow in each cell in the flow pattern is readily determined without the use of any iterative techniques. The saturation and permeability distributions are readily determined. These data and the shape factor, which is measured only once from the potentiometric model when mobility ratio is one, determine the resistance to flow in each cell. THEORY The authors showed in a previous paper4 that, as an aid to calculating performance, the reservoir can be divided into channels using the streamlines of a potentiometric model as a guide. See Fig. 1. This procedure also was used in this paper. The authors also showed that, by treating conduits as approximately one-dimensional and neglecting pressure gradients transverse to the main flow, the Buckley-Leverett equation may be expressed as The principles expressed by this equation are employed extensively in the three-phase flow, as they were in two-phase flow. In the three-phase flow, the channels (the size and shape of which are taken from a potentiometric model)

By G. C. Wallick, W. R. Foster, J. G. Savins

A comprehensive review of the salient features of the differentiation method of rheological analysis in Poiseuille flow from its inception circa 1928 is presented. Here no initial assumptions regarding the nature of the function relating rheological parameters to observed kinematical and dynamical parameters are required in the data-analysis process. In contrast, the integration method involves interpreting flow properties in terms of a particular ideal model. It is shown that, although both methods represent modes of solution of the same integral equation, being relatively bias-free, the differentiation method offers a more discriminating procedure for rheological analysis. The application to problems involving plane Poiseuille flow is also described. INTRODUCTION In most instances, the approach to the problem of interpreting the rheological properties of various compositions as they are affected by changes in chemical or physical environment, assaying the characteristics of a particular constituent of a suspension, analyzing flow behavior in terms of interactions between components in a system, to cite but a few examples, has been in terms of what Hersey1 terms the integration method. Briefly, it consists of interpreting flow properties in terms of a particular ideal model. The usual practice of the integration method is to choose a model with a minimum number of parameters because, other things being equal, it is desirable to use the simplest model which will describe the behavior of a real material and yet be mathematically tractable for the requirements of data analysis. This expression is then substituted into an equation which relates observed kinematical and dynamical quantities, such as volume flux Q and pressure gradient ], and angular velocity and torque T, in a capillary and concentric cylinder apparatus, respectively. The rheological parameters appear on integrating, in an expression relating the pairs of observable quantities such as those just given. In many instances a particular model provides a good representation of rheological behavior over a reasonable range of compositional and environmental changes. Just as often, however, it is obvious that the interpretation of rheological changes by the integration method is not providing realistic information about changes in flow behavior. A more general method of interpreting rheological data for a given material is to make no initial assumptions regarding the nature of the function relating rheological parameters to observed kinematical and dynamical quantities, e.g., flow rate and pressure drop in capillary flow or angular velocity and torque in a rotational viscometer. This general method Hersey1 terms the differentia-ation method. Instead of integrating, one differentiates the integral equation with respect to one of the limits, i.e., one of the boundary conditions; the resulting expression contains the same observable quantities just given, their derivatives, and the rheological function evaluated at that boundary. By obtaining these derivatives from experimental data, graphically or by a computer routine, they can be substituted into the differential equation and a graphical form of the function derived. THEORY OF THE DIFFERENTIATION METHOD FOR POISEUILLE-TYPE FLOWS In this introductory paper, two flow cases which are important in viscometry are considered (one for the first time) from the differentiation method of analysis, flow in a cylindrical tube and flow between fixed parallel surfaces of infinite extent, the basic integral equations being formulated in a manner analogous to the way they originally appeared in the literature. In addition, the following ideal conditions will be assumed: (1) an absence of anomalous wall effects, (2) isotropic behavior everywhere, and (3) steady laminar flow conditions. Before proceeding it would seem appropriate to preface the analysis of these flow geometries with the comment that the basic philosophy of the differentiation method could be said to stem directly from the efforts of Herzog and weissenberg2 to indicate systematically the various ways of using

By H. L. Hartman

The drop tester has proved an invaluable tool for the investigation of percussion drilling in the laboratory in "slow motion". It has allowed the process of rock penetration by impact to be studied a single blow at a time. In the present work, the relation of cratering to indexing has been determined under conditions simulating those in the bottom of the hole at atmospheric pressure. Indexed blows form craters which are influenced by others adjacent to them. Two unique effects are involved: (1) the provision of additional free faces in proximity to the point of impact and (2) the creation of subsurface damage by the previous blow(s). Both have a pronounced influence on the volume of rock removed per blow, which governs the rate of penetration in actual drilling. The emphasis in this study was placed on determining effects of index distance and impact energy on crater volume when blows were struck on a previorusly "drilled" or damaged surface, a situation most representative of down-hole conditions. The results indicated that (I) the optimum index distance on a damaged surface is greater than on a fresh, undamaged rock surface; (2) the volume of rock broken at optimum index distance is about the same for both surface conditions with die-shaped chisels but tends to be greater on a cratered surface with wedges; (3) the optimum index distance and maximum crater volume are proportional to the energy level; (4) at a given energy level, the volume of rock broken by dies generally exceeds that by wedges; and (5) maximum crater volume varies, approximately, inversely with width of die but is nearly independent of included angle of wedge. INTRODUCTION This is the third of a series l,2 of papers reporting on one phase of a continuing program of drilling research being conducted in the Rock Mechanics Laboratory at The Pennsylvania State U. The phase under consideration is crater geometry — blow energy relations in percussion drilling. While the previous papers have discussed mainly craters produced by single blows, indexed-blow studies are the subject of this paper. The basic mechanism of drilling under study here — that of penetration of a chisel-shaped bit subjected to impact — is the one employed in probably 90 per cent of all rock drilling in mining, percussion drilling. It also forms the fundamental action in nearly all oil well drilling with the so-called rotary (roller-bit) method. Here, the action is more complex, however, and may in actuality comprise rotary-percussion drilling. Certainly, the down-hole air and mud hammer drills, which superimpose percussion on rotation, belong in this latter category. In any event, impact blows constitute one of the principal mechanisms in the vast majority of drilling, mining or petroleum, and it is essential in drilling research to study percussion in all its basic aspects. Percussion drilling can be simulated in the laboratory by means of an impact drop tester. Striking one blow at a time with chisel-shaped bits, the drop tester permits study of most of the parameters in the drilling process except blow frequency and the interrelationship of certain variables. Single-blow craters, the simplest to produce, allow the determination of basic relationships, such as volume vs blow energy; while indexed craters most closely resemble the action in real drilling and provide data for the establishment of optimum operating criteria. Indexed-blow craters are those which are influenced by others in proximity. The simplest case to visualize, and the one studied in detail here, is that in which a given crater is formed adjacent and parallel to a previous crater, such that the second crater breaks into the first. A more complex situation usually exists downhole, where indexing is achieved by rotation of a radial-winged or rolling-cutter bit. But, admittedly, the phenomenon of indexing is more likely to be understood eventually if it is first studied in its simplest form, that related to parallel craters. Limited research has been conducted on indexed blows. The most complete study was that by Simon, 3,4 who proposed universal indexing curves to relate the important variables. His work, however, was conducted entirely with rows of craters formed on a

By American Institute of Mining Engineers

By R. F. Nielsen, D. E. Menzie

The object of this study was to determine if crude oil could be produced successfully by a process of crude oil vaporization using carbon dioxide repressuring. This process appears to have application to highly fractured formations where the major oil content of the reservoir is contained in the non-fractured porosity with little associated permeability. Crude oil was introduced into the windowed cell and carbon dioxide was charged to the cell at the desired pressure. A vapor space was formed above the oil, and the crude oil-carbon dioxide mixture was allowed to come to equilibrium. The vapor phase was removed and the vaporized oil collected as condensate. Samples of all produced and unproduced fluids were analyzed. Tests were also performed to evaluate the amount of vaporized oil that can he produced by rocking from a high to a lower pressure. The carbon dioxide repressuring process was applied to a sand-filled cell to investigate the performance in a porous medium. A test was performed to evaluate how the condensate recovery changes as the size of the gas cap in contact with the oil changes. INTRODUCTION This study has been directed toward a relatively new process of vaporization of crude oil designed to increase ultimate production of hydrocarbons through the application of carbon dioxide to an oil reservoir. Suggested advantages of carbon dioxide repressuring of a petroleum reservoir are: (1) reduction in viscosity of liquid hydrocarbons due to the solubility of carbon dioxide in crude oil, (2) swelling of the reservoir oil into a larger liquid-oil volume with a resulting increase in production and decrease in residual oil saturation due to an increase in the relative permeability to oil, (3) displacement of more stock-tank oil from the reservoir since the residual liquid is a swelled crude oil, and (4) gasification of some of the hydrocarbons into a carbon dioxide-hydrocarbon vapor mixture. Balanced against these advantages are several detrimental factors which must be evaluated; i.e., high compression costs and corrosion of well equipment and flow lines. Some of the more outstanding contributions to the study of carbon dioxide injection have been reviewed in order to furnish a basis for a continuation of research pertaining to this method. The literature reviewed1-8 has been limited to that dealing with carbon dioxide repressuring processes or with carbon dioxide-crude oil-natural gas phase behavior. Articles relating to carbonated water injection and literature published on the use of low pressure carbon dioxide gas injection in water flooding have not been included in this study. In 1941 Pirson5 suggested the high pressure injection of carbon dioxide into a partially depleted reservoir for the purpose of causing the reservoir oil to vaporize and thus produce the oil as a vapor along with the carbon dioxide gas. By reducing the pressure on this produced mixture of hydrocarbons and carbon dioxide at the surface, it was proposed to separate the hydrocarbons from the carrier gas. He theorized that essentially all the oil in a reservoir could be produced by simply injecting enough carbon dioxide to vaporize the residual oil. This present investigation deals with the vaporization of a crude oil by carbon dioxide, the molecular weight and gravity of the vaporized oil product and the characteristics of the residual oil after several repressuring cycles with carbon dioxide. An attempt is made to evaluate the merits of a vaporization process for the crude oil rather than a flow process where the oil recovery is determined by relative permeability considerations. Such a vaporization of crude oil by carbon dioxide repressuring appears to have possible use in a highly fractured formation where the major oil content of the reservoir is contained in the non-fractured porosity with little permeability. The carbon dioxide flows into the fractures, contacts the crude oil in the matrix and vaporizes part of the crude oil; this vaporized oil is produced and recovered and the carbon dioxide is reinjected again. The specific problem of this study is to attempt to answer this question; Can crude oil be produced successfully (technically, but without economic considerations) from a petroleum reservoir by a process of vaporization of the crude oil by carbon dioxide repressuring? DEFINITION OF TERMS AS APPLIED IN THIS STUDY Carbon Dioxide Contact: One cycle in which carbon dioxide was injected and bled off. Condensate: The hydrocarbon liquid which was condensed out of the mixture of hydrocarbon-carbon dioxide vapor upon reduction of the pressure of the vapor. Hydrocarbons Produced (HCP): All the hydrocarbon!, which were vaporized by the carbon dioxide repressuring process and were removed from the cell during any specific cycle or carbon dioxide contact. Hydrocarbons Unproduced (HCU): All the hydrocarbons which were not vaporized by the carbon dioxide

By J. E. Miller, W. E. DeVaney, L. Stroud

During a recent phase study of a natural gas, two stable equilibrium liquid phases were observed at temperatures below —200°F and pressures above 200 psi. This paper reviews the published literature on the occurrence of multiple equilibrium liquid phases and presents analytical data for the vapor and two equilibrium liquid phases of the liquefied natural gas at five experimental conditions. In addition, data for 30 conditions of two-phase equilibria are included. INTRODUCTION The low-temperature phase behavior of gases, most of which contained helium, has been investigated in the laboratories of the Helium Activity for many years. Since 1952, experimental studies of these systems have been continuous as part of the research program at Amarillo, Tex. Because of their value to private industries interested in participating in the Helium Conservation Program, several "Open File" reports containing phase equilibria data for heliurn-bearing natural gases have already been released by the Helium Activity. A paper containing information on the general phase behavior, operating criteria and extensive vapor-liquid data for two helium-containing systems was recently published.2 Additional publications presenting experimental data on the phase relationships of various gas systems are now in process and will be available in the near future. PREVIOUS EXPERIMENTAL WORK Although the formation of multiple liquids has been reported for various systems, to our knowledge this paper is the only substantiated evidence of a vapor-liquid-liquid equilibria in a naturally occurring gas. In 1940, Vink, Ames, and others13 reported the presence of two liquid phases in a hydrocarbon system consisting of mixtures of crude oils, solvents and natural gas. Eilerts and co-workers; published data on the recombined fluids from a gas-condensate well. This condensed gas, containing approximately 76 per cent methane and 24 per cent ethane-plus, exhibited two distinct liquid phases. Weinaug and Bradley 14 observed "unusual" phase behavior in a reservoir mixture. These workers postulated that the anomalous phase behavior was due to the "imminent formation of a second liquid phase". Botkin, Reamer, Sage and Lacey 1 studied two California crude oils that exhibited multiple phases. Kohn and Kuratal0 recently reported two equilibrium liquid phases in the methane-hydrogen sulfide system. Roof and Crawfordl1 and Eakin, et a16 also have reported experiments with binary systems that formed two stable equilibrium liquid phases. APPARATUS AND PROCEDURE A U. S. Bureau of Mines Phase Equilibrium Apparatus was used in conducting this study. The apparatus and procedures employed in its operation have been previously described3 and will not be repeated in detail ill this report. Briefly, the apparatus consists of a windowed cell which can be maintained within ±0.5°F for temperatures between room temperature arid -320°F. Pressure within the cell can be maintained within 0.1 per cent of gauge reading up to 800 psig. Equilibrium vapor and liquid samples are obtained in special containers4 for analysis by a mass spectrometer. Although the accuracy of the analyzer is about ±0.1 mol per cent, the reproducibility of the phase-equilibrium apparatus is considered relatively poor. Values reported from methane and nitrogen are considered accurate to ±l.0 and ±0.6 mol per cent. Data for ethane-plus in the vapor are accurate within 0.2 mol per cent; liquid-phase data for this aggregate compolnent are accurate within 1.5 mol per cent. For helium in the vapor phase, the analytical data are accurate within 0.2 mol per cent; liquid-phase analyses for this component were obtained by the charcoal adsorption method described by Frost8 and are accurate within 0.006 mol per cent. All of these references to the accuracy of reported values are conservative estimates based upon a statistical treatment of reproducibility data obtained with the apparatus.

By William P. Blake

The occurrence of siderite in beds of considerable thickness in the clay formations of Martha's Vineyard, Mass., may have some economical importance, and is at least interesting in a scientific point of view.

By H. D. Yoo, M. R. Tek, D. L. Katz

The behavior of gas-storage reservoirs subject to water drive is investigated through analog simulation on an electronic differential analyzer. The simulation technique developed on an LM-10 computer permits the prediction of reservoir volume or pressure resulting from the movement of water in the surrounding aquifer. The method developed on the analog computer consists of setting up an appropriate transfer-function circuit and feeding the arbitrary time-varying boundary conditions as an input signal. The input may be specified as gas reservoir pressure, pore volume or the water flux. Several cases studied include an isolated gas reservoir on a limited aquifer, interference among three reservoirs adjacent to a common aquifer and the growth of gas-storage volume on an aquifer. It is concluded that the method developed on an electronic differential analyzer provides an excellent technique to simulate and investigate the behavior of gas reservoirs subject to water drive. The agreement between the reservoir performance as predicted from the simulation technique and as measured from actual field data is found to be better than the range usually encountered in predicting water-drive behavior. INTRODUCTION It is generally known that some gas-storage reservoirs are located on top of blanket sands of large extent, saturated with brine called aquifers. Because the volume of the body of water associated with aquifers is usually very large and water is compressible, the cyclic pressure variations encountered in normal storage service inevitably cause unsteady, compressible flow conditions in the adjacent aquifers. The solution to radial diffusivity equation for a limited aquifer for constant terminal conditions has been know111 since the early 1930's. The solution for the constant terminal conditions for an infinite aquifer was puhlished in 1949.2 The proldem of handling the unsteady flow of water through a porous medium with time-varying boundary conditions has been studied by many investigators through the use of graphical,3 digital4 computing techniques or on an R-C type of electrical analog called a reservoir analyzer.5 The digital computer performs a series of arithmetic calculations to superimpose the solutions of constant terminal conditions, while the R-C type of analyzer permits the simulation of the continuous medium by means of "lumped-cell" electrical resistance and capacitance elements based on finite-difference approximations. The graphical method is based on the convolution integral solution of the partial differential equalion. The purpose of the present study was to simulate the behavior of an aquifer on a general-purpose, electronic differential analyzer to obtain the influx of water or the change in pressure at the reservoir continuous1r and quickly. The use of the electronic differential analyzer was based on well known solutions for the constant terminal conditions already available in the literature.2,6-8 THEORETICAL DEVELOPMENT MODEL. The gas reservoir (R) is described in relation to an aquifer in Fig. 1. The permeability and porosity of an aquifer sand, the viscosity of the brine, and the compressibility of the brine and the formation are assumed to be constant. It further is assumed that the blanket sand is not thick enough to give appreciable vertical pressure distribution. The model in consideration is circular with radial geometry. The aquifer radius can be of limited or infinite extent. The gas-water interface is assumed to be fully segregated, and the displacement of the gas by

By R. Neilson Clark

THE coal-beds of Canyon City are situated six miles below the town, upon the Arkansas River. At this point the Rocky Mountains have thrown out from their main ridge two spurs. The one to the north, containing Pike's Peak and Cheyenne Mountain, is known as part of the "Rim Range;" the one to the west and south is known as the "Greenhorn Range." Against these syenitic mountains, the later rocks are piled and uplifted ; at right angles to the direction of this main uplift is the axis of a lesser plication. On the north, west, and south, the pitch is steep, gradually flattening as it approaches the centre of the basin and rising again to the east at a pitch of about 5", thus forming a perfect basin, the centre of which is only three miles east of the base of the Greenhorn Range.

By H. R. Brannon, R. J. Robinson, G. W. Kattawar

Accurate prediction of storm wave characteristics are needed for design of offshore structures. Statistical methods of random noise analysis provide techniques for predicting required wave properties. These techniques have been used to analyze and characterize storm wave profiles from Gulf of Mexico recording installations in 34, 65 and .98 ft of water. Correlations based on the results can be used to predict wave crest probabilities and wave shapes for a range of Gulf water depths and storm conditions. Example predictions of wave crest probability as a function of water depth for a particular set of storm conditions are given. INTRODUCTION Accurate predictions of wave crests and wave shapes are needed for design of offshore structures. With predictions of wave crests, platform deck elevations giving adequate protection from exposure to storm waves can be selected. Use of hydro-dynamic wave theories with the wave crest and wave shape predictions permits an estimate of the wave forces an offshore structure must withstand. Observation of ocean waves suggests use of statistical analysis for studying wave characteristics for water depths beyond the breaking wave zone; recent wave generation theories lend support to this approach. l-3 One of the most powerful methods of statistically analyzing wave profiles was introduced by Cartwright and Longuet-Higgins; most of the basic relations used were developed by Rice5.6 in studies of random noise theory. This method predicts wave crest elevation probabilities from two parameters characterizing a wave profile. The statistical approach also provides a means for predicting wave shapes. In the random noise model of a wave profile, surface elevation is represented as an infinite sum of sine waves with closely spaced frequencies and random phase angles. Power density is proportional to the sum of the squares of the amplitudes of sine waves having frequencies within a narrow band, and the distribution of power density as a function of frequency is called the power density spectrum, or power spectrum. The power spectrum characterizes an irregular sea, and hence finds use in motion studies of ships, barges and semisubmersible drilling platforms. From the power spectrum, a wave profile of any duration of time can, in principle, be calculated; this long wave profile depicts the many wave shapes to which a structure may be exposed. Thus, the statistical methods of wave analysis provide an approach to selecting wave shapes as well as wave crest elevations needed for design of offshore structures. For practical application of these techniques, power spectra and related quantities must be predicted from storm properties. Barber and Tucker 7 have reviewed correlations of wave properties and storm conditions. Their review summarizes work of Darbyshire,8 among others, in which a correlation of power spectra with wind intensity of the North Atlantic was developed. These results, plus more recent work by Pierson and Moskowitz and Kitaigorodskii,10 establish feasibility of correlation, but there is no theoretical basis for modifications extending the relations to other areas as remote as the Gulf of Mexico. Hence, to predict characteristics of Gulf of Mexico waves, a study of local observations is needed. This paper presents (1) a practical approach to computing wave profiles hat depict shapes of waves for use in force calculations, (2) a summary of relations for predicting probabilities of wave crest elevations, (3) correlations of parameters needed to apply these methods in the Gulf of Mexico and (4) examples of application of the techniques. THEORETICAL BACKGROUND The following sections summarize relations needed to calculate wave profiles and to estimate wave crest probabilities. These relations also identify quantities that must be related to storm

By Thomas Egleston

The ores of mercury of North California are composed of metallic mercury and cinnabar. They are found in serpentine, and are very often associated with chalcedony, in masses more or less irregular, often concentrated enough, however, to furnish ores yielding from three to ten per cent., and sometimes richer. This deposit makes its appearance in Vallejo where it has been worked. North of here the mines are more developed. Most of the quicksilver mines, however, are situated in Sonoma and Napa Counties. On its outcrop the serpentine rocks have become decomposed, and have often been washed away to a considerable depth, so that in many cases in their neighborhood, what is apparently nothing but ordinary dirt, will frequently contain from two to three per cent. of metallic mercury, with but a trace of cinnabar, in which case, as at the Sonoma mine, it is made up into adobes and distilled. Such material does not require to be mixed with dirt to be made into adobes. Very often, also, the outcrop of the rock, where it is not decomposed, is filled with metallic mercury, so that by striking a pick into it, as at the Rattlesnake mine, a pound or more of mercury at a time will sometimes spurt out. Such rock as this is found in several localities in every stage of impregnation, and usually makes very rich ore. There is generally, however, very little of it, and it is found only in the first workings; the ore in depth is always cinnabar. At the Rattlesnake mine near Pine Flat, where large quantities of metallic mercury are found, the rock contains so much petroleum that it has been necessary to make special arrangements to burn the carbides of hydrogen, since the distillation of the petroleum causes an extra quantity of poor soot to be formed in the condensation-chambers. At the Geysers the ore is associated with large quantities of sulphur and gypsum, so that in a hand specimen there is often more sulphur than cinnabar, which is a serious impediment to working, especially for the modern style of furnace with iron condensers, and causes so much soot to be formed that it has been known to penetrate as far as the blower, and to so completely clog it as to prevent its revolution. The ore coming from the mine is more or less hand-picked. In one or two cases attempts have been made to treat the poorest, and especially the very fine ores, mechanically as at the California works.

By R. W. Raymond

The specimen of anthracite which I exhibit is from the Ortiz Mine Grant, about fifteen miles southwest of Santa Fe. The beds belong to the lignitic formation of the Galisteo, which Hayden and

By T. Sterry Hunt

The Jones Mine (or Johannes Mine, as it was originally called, from a former proprietor), situated near Springfield, in the township of Caernarvon, Berks County, Pennsylvania, has long been known as a large deposit of magnetic iron ore, associated with more or less copper pyrites, and with malachite and chrysocolla. It occurs in a belt of soft schistose rocks which rest upon the Laurentian gneiss of the Welsh Mountain,' and to the northward are soon concealed by the unconformably overlying Mesozoic sandstone of the region. These metalliferous strata belong to the Primal slates of Rodgers, which, in many other localities in this region, include similar deposits of iron ores associated with more or less copper, and often accompanied by chloritic and argillaceous schists, and, more rarely, by serpentine and small lenticular masses of crystalline limestone. To this same horizon belong the iron ores of Boyertown, Fritz's Island, and Roudenbusch, near Reading, Wheatfield, Cornwall, and Dillsburg (all of which are found along the northern border of the Mesozoic belt of