Search Documents

By J. Douglas, W. T. Ford, G. E. Henderson, I. F. Roebuck

An implicit numerical method is presented for simulating the differential and algebraic relations governing one-dimensional three-phase flow in porous media. The method is based upon compositional representation of the hydrocarbon system. Variable physical properties and water-oil capillary forces are included in the formulation; however, the effect of gravity is ignored. The effects of changing composition and mass transfer are considered through the use of known phase behavior concepts and correlations. The differential and algebraic equations and the numeric approximations to these equations are presented, and the computing algorithm is discussed in detail. This method is used to simulate a single-phase, two hydrocarbon-component gas displacement laboratory experiment. The method is also applied in simulating a solution gas drive and a gas injection problem using a reservoir oil represented as a nine-component mixture. For comparative purposes, the latter two problems also are examined with a one-dimensional volumetric model. The results of the simulation of the laboratory experiment demonstrate the relevance of this method to problems wherein the effects of mass transfer significantly enter the displacement process. The comparison of the results of the volumetric and compositional methods applied to a reservoir system illustrates not only the advantages and benefits derived from examining the hydrocarbon system composition, but also the necessity of this approach in application to a general reservoir problem. INTRODUCTION The advent of large memory, high-speed digital computers enabled reservoir engineers and mathematicians to pool their knowledge in the development of sophisticated techniques for the prediction of oil and gas reservoir performance. These techniques have progressed from volumetric and trend analysis calculations to the development of "mathematical models", or reservoir simulators. These latter models, in general, utilize finite difference approximations to the rather complex partial differential equations that mathematically describe the physics and thermodynamics of fluid flow in porous media. Many forms of this "model" type are available today, and much of the recent literature has been devoted to either improving known methods or developing new techniques for use in these prediction methods. The more advanced techniques in use at this time consist of the solution of differential and algebraic relations that describe a reservoir geometry in one, two or three dimensions and that incorporate reservoir hydrocarbon and water fluid systems that obey specified physical laws, having properties that might easily be measured in the laboratory. This type of model utilizes formation volume factors and other volumetric data, expressed directly as functions of pressure, to represent the reservoir fluid system. Experience gained with this type of simulator, commonly known as a "beta-factor model", has indicated that the approach has many outstanding defects, most of them due to the assumption that those properties measured in the laboratory completely describe the behavior of the fluid system throughout the life of the project under study. The effects of varying hydrocarbon composition and mass transfer between the phases of this hydrocarbon system must, by definition, be ignored. These defects, which may often be ignored in

Jan 1, 1970

By W. F. Rochow, W. C. Burke

Numerous designs of linings and accessories, including dams, lifters, and heat exchangers, contribute greatly to kiln efficiencies. Greater conductivity is achieved with basic brick than with fireclay and high-alumina refractories, depending largely upon the use of steel plates in radial joints between adjacent brick to increase heat transfer through the lining. Refractory lifters expose all portions of the load to radiant or convective heat and raise some portions to higher temperatures by cascading the load down through the path of the hot gases.

Jan 3, 1955

By H. J. Gerritsen

From presently available components a portable, rugged, reliable apparatus can be built which will be able to detect methane concentrations of 0.1% and lower in air. Sensitivity and design considerations of the apparatus are given. Ways to build even simpler and cheaper apparatus in the future are indicated. Traditionally the problem of detecting and monitoring the presence of dangerous gases present in small concentrations in air has been a difficult one. Whether canaries,' or sophisticated chemical or thermal conductivity methods are used, they seem to be either slow, cumbersome, or otherwise not completely satis factory. A seemingly natural way of gas detection is to measure the strength of absorption undergone by light of a wavelength that the gas to be detected absorbs. This method has been difficult to apply in view of the low power output over a narrow wavelength range in the infrared from traditional glow bar light sources. Besides, the apparatus to select such a narrow band of wavelengths is quite costly and complex. The use of a laser solves these problems simultaneously. Although a variety of lasers are available to detect a variety of gases (such as carbon monoxide or methane) whose early detection would be desirable, this paper deals exclusively with the use of lasers for methane detection. ABSORPTION SPECTRUM OF METHANE IN AIR The absorption spectrum of pure methane is very well known.2-4 There are four fundamental vibrations, all occurring in the infrared of which only two absorb strongly. They occur at 3.3 µ and 7.7 µ. The absorption spectrum is quite complicated due to the fact that the methane molecule can rotate in a variety of rotational quantum states. These rotational states are furthermore split by anharmonic forces such as the Coriolis force.4 When methane is present in air at atmospheric pressure, much of the very fine structure due to the anharmonic force is washed out since the pressure broadening is of the same magnitude as the splitting due to anharmonicity. For most of the spectrum at 3.3 µ, the rotational lines are typically separated by 0.01µ while the individual linewidth is 0.0001 µ. It is thus very important that the laser emit at a wavelength that coincides with one of these rotational lines. Before elaborating on the point that such a laser exists and how it performs, it is necessary for later discussion to point out that in the

Jan 1, 1967

By L. W. Pommier, F. B. Brien, A. K. Bhasin

In the field of mineral processing, although many operations are applied to suspensions of ore particles in fluids, there is a lack of fundamental knowledge relative to the flow resistance encountered by immersed bodies moving through solid-liquid suspensions. Previously proposed equations, based on the Ostwald-de Waele (power-law) model for non-Newtonian fluids, are exam ined in the turbulent flow region relative to a number of pertinent variables and their effects on the movement of immersed bodies. As more knowledge is acquired in regard to the effects of suspension variables on the rheological properties, useful relationships may be proposed to allow predictions of conditions to achieve optimum results for engineering purposes. In the field of minerals processing, many operations are applied to suspensions of ore particles in fluids. The differential movement of solids of various sizes and densities in suspensions is of utmost importance to gravity concentration, wet classification, flotation, hydrometallurgy, and pipeline transport of solids. The subject is also of importance to other branches of engineering; for example, those involved with soil mechanics, ceramic slips, and oil-well drilling fluids. There is a lack of fundamental knowledge relative to the flow resistance encountered by immersed bodies moving through solid-liquid suspensions. The problems encountered in the direct application of the well-known Stokes and Newton's Laws to non-Newtonian systems have been examined in a previous paper.' Equations were proposed based on the Ostwald-de Waele (power-law) model for non-Newtonian fluids. This present paper examines these equations in the turbulent flow region relative to a number of pertinent variables and their effects on the movement of immersed bodies. It is thought that, as more knowledge is acquired as to the effects of suspension variables on the rheological properties, useful relationships may be proposed to allow prediction of conditions to achieve optimum results for engineering purposes. The purpose of this paper is to make an initial assessment of the nature and magnitude of pertinent suspension variables. In order to avoid too voluminous a document, the examination of kinetics and mechanisms will be omitted at this time. THEORY A Newtonian fluid is defined as a fluid in which the stress is proportional to the rate of strain, 2 i.e.: H = viscosity; see Table I. for list of symbols The following equations are conventionally used to assess the motion of spherical particles moving under laminar and turbulent flow conditions through Newtonian fluids: When these relationships are applied in practice to suspensions, the liquid density is substituted by a composite suspension density, and viscosity is usually taken as that of the suspending liquid or an approximated "apparent viscosity." However, suspensions are non-Newtonian in character3 and thus, more than one viscosity parameter is required to describe their behavior. A number of models have been proposed for non-Newtonian fluids;2,4 among these the simplest is the Ostwaldde Waele or power-law model for incompressible fluids which for one dimensional flow may be reduced to:

Jan 1, 1970

NEW MEMBERS The following list comprises the names of those persons who became members during the period Oct. 10, 1916 to Nov. 10, 1916. ANCHOR, HANS C., Supt., Dome Extension Mines Co., Ltd., South Porcupine, Ont., Canada. BABCOCK, EARLE J., Dean, College Of Min. Engrg., State Univ. Of N. Dak., Grand Forks, N. Dak. BIDDIS0N, P. McDONALD, Construction Engr., The Ohio Fuel Supply Co., 52 W. Gay St., Columbus, 0. BOISE, PAUL T., Asst. Mgr., Stimpson Equipment Co., Felt Bldg., Salt Lake City, Utah. BOURNE, RALPH H., Vice-Pres Whiting Foundry Equipment Co., Harvey, Ill. BROWN, EDWARD OTTO FORSTER, Min. Engr., 706-707 Salisbury House, London, E. C., England.

Jan 12, 1916

By M. N. Brodie, John Chapman

INTRODUCTION - DEFINITION & HISTORICAL PERSPECTIVE While this symposium is concerned with comminution which, by definition, relates to all stages of size reduction of an ore body to a form suitable for mineral separation, the electrical drive aspect of the grinding section will receive most attention in contrast to the crushing, screening and materials handling elements of the process. The reason for this is that today's large mills, thirty feet diameter and upwards, present mechanical and electrical problems which require a degree of sophistication in their solution not previously required. The ability to satisfy these demands, so far as electrical technology is concerned, largely results from the development of semi-conductor devices. These can now handle typical electrical distribution power levels within remarkably s m all volumes and with high efficiency when compared with the motor-generator sets, mercury arc rectifiers, rotating amplifiers, magnetic amplifiers and vacuum tubes which were commonly used to provide speed and torque controls as recently as the early nineteen-sixties. The hardware available to the electrical engineer has changed dramatically whereby even the simplest control system formerly using hard-wired relays is now superseded by the program m able logic controller with its computer-like memory containing the sequence of operation, in other words, the software. Reverting now to the mineral process itself, it is assumed that the ore body has been mined using explosives and is delivered to the concentrating plant in its 'run of mine' form. The problem becomes one of reducing the material to a size which will permit an acceptable recovery of the mineral with the most economical use of energy. The size of 'run of mine' material will range from about minus five feet for high tonnage open pit operations to minus two feet for the smaller underground mine. The final reduction in size may range from minus ten mesh to minus two hundred mesh (exceptionally as fine as minus six hundred mesh) depending upon the concentration method and dissemination of the mineral values within the ore. There is no single circuit that satisfies all applications, but the pattern of metallurgical equipment use shown in Table 1 appears discernable from current practice, in general terms.

Jan 1, 1982

By H. D. Outmans

A method has been developed for determining the relative water wet-tability (fraction of the surface wet by water) of porous media. This method involves the adsorption of methylene blue dye from an aqueous solution onto the solid surfaces of rock contacted by the injected dye solution. Experiments on Berea snndstone cores have shown that water-wet cores adsorb a large amount of dye from solution. In oil-wet cores, however, the injected aqueous dye solution is prevented from contacting the adsorptive surfuces by a film of oil and virtually no dye is adsorbed. This method was tested by experiments on mixtures of water-wet and oil-wet sand. Using the dye adsorption method, it was possible to determine quantitatively the fraction of water-wet and oil-wet sarzd in such mixtures. The dye adsorption test is readily applied to determinitrg the relative water wettability of fresh oilfield cores. The cores are not extracted before testing. Observed values of the wettability of field cores tested to date have varied from 100 to 6 per cent water wet. INTR OD UCTIO N Recent articles1 ' have emphasized the important role that wettability plays in almost all phases of reservoir rock behavior. It has been shown that wettability has a profound effect on the interstitial water saturation, residual oil saturation, capillary pressure, relative permeability, waterflood behavior, and resistivity index of oilfield cores. A quantitative investiga- tion of the effect of wettability on these core properties is severely handicapped by the lack of a simple technique for measuring the wettability of a porous medium. If such a technique were available, it would provide investigators with a means of gaining insight into the distribution of phases within a core. A number of qualitative methods of estimating wettability are known. These include contact angle measurements, flotation tests and studies of displacement pressure: imbibition7 and relative permeability data.' Using these tests, it is usually possible to classify a core as water-wet, oil-wet, or some intermediate degree of wettability. One recently proposed method involving measurement of nuclear magnetic relaxatione offers promise of determining quantitatively the fractional wettability of cores. The equipment is rather elaborate and necessary treatment of the cores before testing probably alters the natural wettability. This latter point also applies to most of the cited qualitative methods. This paper presents a new technique for determining water wettability of cores; it is based on measurements of the dye adsorption capacity of the core. This method has the advantage of requiring relatively simple equipment and involves minimum handling of the core before testing. It is important to distinguish between the various concepts of wettability referred to in the literature. The term "wetting" is defined by Bartell as that phenomenon which occurs when a solid and liquid phase come in contact in any manner so as to form a solid-liquid interface. A preferentially water-wet solid is defined by Jennings in terms of the advancing contact angle in the solid water-oil systems. Brown and Fatt regard fractional wettability as the fraction of the total surface area that is preferentially oil-wet or water-wet. However, none of these definitions adequately describes the type wettability measured by the dye test. This calls for the introduction of a new term, "relative water wettability", which we define as the fraction of the total surface area of a core contacted by injected water. DESCRIPTION OF METHOD The ability of a porous medium (especially the clays) to adsorb large quantities of surfactants, dyes, etc., is well known. It is reasonable to assume that a substance can be adsorbed from solution onto a solid surface only if the solution contacts the surface. In this connection Shapiro" has shown by experiments on silica gel that completely dry gel will adsorb methyl red dye from a benzene solution. On the other hand, if a film of water more than several molecules thick is present on the silica gel, the adsorption capacity of the gel for the dye is drastically reduced. Thus, a study of the dye adsorption capacity of cores should give some clue as to the fluid distribution on the surface of the core. TEST OF METHOD As a qualitative check on this theory, dye adsorption tests were run on similar Berea sandstone cores (Table 1, Cores 1-5). Cores saturated initially with water and flooded with 0.01 per cent aqueous methylene blue adsorbed 0.58 mg dye/gm of core. Similar cores saturated initially with water, driven to irreducible water with oil, and then flooded with the same dye solution also adsorbed 0.58 mg dye/gm of core. Another core, treated with a 5 per cent Dri-Film solution, dried and saturated with oil (Soltrol) adsorbed no dye when the dye solution was injected. These tests confirm

By H. A. Pearse

NORMALLY, the Britannia ore mixture contains chalcopyrite and pyrite as the chief sulfide minerals, together with minor amounts of gold and silver and a low zinc content. Reduction is accomplished by jaw crushers underground, followed at the concentrating plant by Symons cone crushers, rolls, and ball mills. The standard flotation practice, evolved during several years of experimentation, provided for bulk flotation of the two sulfides, after primary grinding to only approximately 35 mesh, followed by regrinding of the bulk concentrate, with subsequent separation of the chalcopyrite and pyrite to make a marketable concentrate of each mineral. This so-called standard milling practice has previously been described in detail in U. S. Bureau of Mines Information Circular 6619. A modification will he discussed here. In the face of a constantly increasing world surplus of copper and a steadily falling price during 1930 and 1931, it was obviously undesirable to continue unprofitable production of large quantities of this metal. On the other hand, company officials were

Jan 1, 1934

By Norman F. Schulz

The technical excellence of separation achieved in a mineral concentration process, or any other process where two constituents of any kind are physically separated from each other, is expressed uniquely and quantitatively by the separation efficiency (E,): E. = (R - R,), where R is the percentage of the valuable constituent and R, the percentage of the waste appearing in the concentrate. The separation efficiency meets all the usual requirements of a true efficiency and is independent of the nature of the process used. It can be visualized as the percentage of the feed material that, in effect, actually undergoes complete separation while the rest of the feed is distributed randomly and unchanged into the two separation products. The problem of quantitatively assessing the results of a mineral concentration process recurs whenever one tries to decide which of two similar test results is more desirable. Although economic factors may greatly influence the eventual choice of process, they usually have very little to do with the technical excellence of a separation and are, therefore, not a consideration at this point. Concentrate grade and recovery, used simultaneously are the most widely accepted measures of performance, and if both the grade and the recovery are greater for one case than for another, the choice is simple. However, when the results of closely competitive operations are compared, it often occurs that one set of test results shows a higher grade and a lower recovery than another and the choice is no longer obvious. Laboratory investigations of a mineral concentration process applied to a given feed material can often be designed to yield successive paired values of concentrate grade and recovery over the practical range of the course of a separation process.' These results may be plotted as some combination of the cumulative values of grade, weight recovery, or unit recovery in the concentrate. The resulting correlations for two slightly different sets of conditions can then be compared qualitatively, or even quantitatively in particular cases.2 Such comparisons are not applicable, however, to individual concentration test results. The reduction of a separation test result to a single efficiency figure, preferably as a percentage of ideality, would afford the simplest and most precise means for comparing the results of competitive tests. The effective proportion of feed material that undergoes complete separation in any real concentration process, hereafter referred to as "separation efficiency," is such a figure. It is essily calculated from the usual feed and product data and can even be evaluated, in many cases, from data in publications describing various concentration processes. Mineral Concentration Processes and Their Evaluation Mineral Concentration Processes: The overall concentration process may be subdivided into a liberation process which severs the bonds between unlike mineral species, and a separation process, which sorts the resulting liberated particles into two products—a concentrate rich in mineral A and a tailing lean in mineral A. A perfect mineral concentration process is represented by the idealized scheme shown in Fig. 1 in which liberation and separation are not differentiated. The feed mixture consists of two mineral types which are completely liberated early in the process and then sorted by type. All the mineral particles of type A which respond uniquely to the selection forces acting in the ideal separator, report to the concentrate, while all the mineral particles of type non-A report to the tail-

Jan 1, 1971

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 W. Hurst

The authors of this paper, in substituting values in a material balance equation associated with a gas reserve and a water drive, cannot find sustainment of what the gas in place should be; therefore, their conclusion is that this type of problem dealing with gas reservoirs is insoluble. They have introduced the word, "cybernetics," based on the work of Norbert Wiener, that further extends this line of reasoning, namely, to the effect that the system is a "black box", and that the gas reservoir and aquifer can take on infinite aspects 'based upon their past performance, to quote the authors. The inference is that this black box is an unknown. which permits in computer programming to subscribe as many variables as possible to observe by adjustment of one to another those combinations of conditions that fit field behavior. This may satisfy those in this discipline whose dedication is to programming, but to the others this is a comparison to the rationalization of the combinational results of the throwing of dice—a probability procedure. In reservoir engineering we are more direct. We learn from the examination of individual we'll and field data exactly what is in this black box. I now come to the impasse that has confronted these authors in their machine calculations. In this I have reproduced from their own production data (Tables 1 and 2 of their paper) pressure vs cumulative gas produced. shown as Figs. D-1 through D-5. A perusal of any of these five plots, for the corresponding fields in the Po Valley, shows the straight-line relation of pressure vs gas produced, evidencing pressure depletion. The point I make is that in their fundamental production data, there is no variance from a straight line that can support a water drive.

By W. T. Swensen

The Anaconda Company, through its wholly owned subsidiary The Anaconda Company (Canada) I,td., has investigated, optioned and partially developed a large tonnage of iron ore amenable to magnetic concentration. The reserves are located in north central Ontario and occur within a twenty-two mile belt of layered and intricately folded metamorphic rocks, close to the north contact of an igneous complex of granite and pegmatite. There are two main ore deposits. Briarcliffe is a steep to vertically dipping body 100 to 500 feet wide and over a mile long. The Two Mile Lake deposit is flat dipping and confined to an open, asymmetrical syncline; it is approximately a mile wide and 250 to 4.00 feet thick. On its flat easterly plunge, it is saucer-like in shape and is amenable to open-pitting for an east-west distance of several thousand feet. Magnetometer surveys, confirmed by drilling, indicate at least a three mile eastward extension of this deposit. Location and Access: The property is located in the Kowkash Mining District of northwestern Ontario, roughly forty-four miles north of the community of Nakina on the main line of the Canadian National Railroad. (Fig. 1) Access during the exploration period was entirely by airplane from Nakina. (Fig. 2) More recently, Nakina itself was connected to the highway network in the Longlac-Geraldton area. The property is roughly 190 miles by existing and proposed railroad from the Great Lakes waterway system. The Trans-Canada pipeline moving natural gas from western to eastern Canada passes through Geraldton roughly 70 miles south of the property. Topography: The surface features in the area are typical of a sizeable portion of the Precambrian Canadian Shield area. Many lakes, low swampy areas, and glacial moraine-type features are common in this recently glaciated area. The average elevation is between 800 and 900 feet above sea level, and the relief is probably less than 100 feet with the greater portion of the area less than 50 feet. The topographic "highs" are formed either by outcrops that were resistant to glacial erosion or by moraines, eskers and drumlins which were deposited upon the recession of the glaciers. The better mineralized zones, composed chiefly of quartz and magnetite, commonly form ridges resistant to erosion and weathering. Previous Work: Because of the obvious outcrop of the magnetite-rich beds and knowledge gained from early developments within the general region, it is conceivable that the presence of the iron formations has been known for a good many years. The area is included in studies made by L. F. Kindle* in 1929 and 1930 for the Ontario Department *Kindle, L.F., Kowkash-Ogoki Gold Area, Ontario Department of Mines Annual Report, Vol. XL, Part IV, 1931; map 40f

Jan 1, 1961

By J. H. Swisher

An experimental study has been made of the possible mechanisms for foam stability in the system CaO-SiO2-Cr2O3, where Cr2O3is the foaming agent. The degree of lowering of surface tension by Cr2O3 was determined at 1600 "C for high -silica melts. Measurements were also made of foam stability under standardized conditions. The results of the two sets of measurements then were correlated by applying the Gibbs and Marangoni theories of film elasticity. It was demonstrated that the Marangoni elasticity effect is probably the largest single contributor to foam stability, although the high uiscosity of' silicate melts makes some contribution in controlling the rate of drainage of liquid from bubble lamellae. ThE tendency of metallurgical slags to foam has frequently been observed during steel-refining operations. In the open-hearth furnace,this tendency is considered undesirable because it slows down the rate of heat transfer from the burners to the molten-metal bath. On the other hand, it is sometimes desirable to have a "foamy" slag in the basic oxygen converter.' A more effective blanket over the steel bath is made by a "foamy" slag, and loss in yield of metal due to splashing is minimized. In addition, a flush slag can be removed more easily in this condition. Since slag-foaming problems in the field have been attacked largely on an empirical basis, it appeared desirable to obtain a better basic understanding of the mechanism of foam stability. The experimental work was directed almost entirely toward the behavior of Cr2O3 as the foaming agent in CaO-SiO,-Cr2O3 melts. In the first stage of the investigation, the degree of lowering of surface tension by Cr2O3 was determined at 1600°C. The second stage consisted of measuring foam stabilities under standardized conditions. The two sets of experimental results then were correlated by applying various theories of foam stability. For convenience in interpreting the surface-tension and foam-stability results reported here, a portion of the CaO-SiO2-Cr2O3 phase diagram2 is reproduced in Fig. 1. EXPERIMENTAL 1) Surface-Tension Measurements. Very little information is available in the literature on the lowering of surface tension of slags by foaming agents. Kozakevitch 3 mentioned that the addition of 1.5 pct Cr2O3 lowered the surface tension of FeO slightly. In another system, some data on the lowering of surface tension of foaming slags by P2O5 were reported by Cooper and Kitchener.4 The technique used for the surface-tension measurements was the maximum bubble-pressure method,' using a single tube immersed to various depths in the melt. This method has been used extensively in high-temperature systems, such as in liquid metals, slags, and glasses.6,7 It has the advantage that the contact angle between the liquid and the tube material need not be known. The equation used for calculating surface-tension values from maximum bubble-pressure measurements is the Schroedinger equatioq5 which can be written as follows:

Jan 1, 1964

By D. R. Parrish, F. F. Craig

In some respects, forward combustion (alias: fie-flooding, underground combustion, in-situ combustion, dry combustion) is an efficient oil recovery method. For example, its displacement efficiency (as much as 95 percent in reservoirs with high porosity and high oil saturation) is good. In its use of heat, however, forward combustion is grossly inefficient. Most of the heat generated during forward combustion is wasted. About one-fifth of the heat is picked up by the incoming air and re-used (regenerated); the remainder is left behind the combustion zone and ultimately is lost. Compressed air is expensive, and heat generated from compressed air via underground combustion is too valuable to be thrown away. Improvement is needed. Since the early 1950's, improvements in forward combustion by addition of water injection have been discussed, particularly in the patent literature, with increasing frequency.'-'' These concepts are based on the fact that relatively inexpensive water can be used to pick up much of the heat that is ordinarily wasted during forward combustion, and transport it forward to a more useful place, thereby reducing the amount of air required to move a heat front through the reservoir. Several approaches have been proposed. Merriam and Squires1 taught that a combustible (fuel) gas should be injected along with the air and water. This, however, would give at least a tendency toward reverse combustion instead of forward combustion, and would reduce the velocity of the heat bank, defeating the primary purpose of the water injection. Pelzer2 described primarily a batchwise method of forward combustion and water flooding. He advocated burning part way through a reservoir with forward combustion, then switching from air to water injection alone to push the heat toward the producing wells. Though possible, this method has serious disadvantages, the major one being that, by the time water injection is started, much of the heat has been irretrievably lost to the overlying and underlying formations. We have described what is considered a practical Combination Of Forward Combustion And Water-flooding (COFCAW).3,4 Air and water are injected simultaneously (or alternately) after a small heat bank has been formed by forward combustion. Some of the injected fluid fills up the region behind the heat bank; the rest enters the heat bank. Upon encountering heated rock, the water is converted into steam, which then flows ahead of the combustion zone and displaces much of the oil in place. Only enough air is injected to provide heat for vaporization of the water and to compensate for reservoir heat losses. Since the size of the heat bank is minimized, heat losses also are minimized. Because of more efficient oil displacement ahead of the combustion zone, fuel consumption is reduced. Too, unlike a forward combustion zone that cannot move until all fuel at a point has been consumed, COFCAW leaves some

Jan 1, 1970

By Frank E. Woolley, Robert D. Pehlke

Jan 1, 1965

By J. W. Givens, P. B. Crawford

A potentiometric model study has been made of the effect of vertical fractures existing in the matrix of the reservoir on the flooding or cycling performance. Fractures can have unusual flow characteristics. Fluid entering one side of a fracture can emerge on either the same or opposite side of the fracture, depending on the particular streamline. Fracture orientation has a great influence on the sweep efficiency. However, sustained fluid injection may still permit large areas to be swept. INTRODUCTION One can find considerable information available in the literature on the effect of fractures originating at the wellbore on sweep efficiencies of water-flooding or gas cycling programs.1-5 However, there arc few, if any, quantitative data reported in the literature on the effect of vertical or horizontal fractures existing out in the reservoir matrix on secondary recovery performance. Several papers have been presented on the effect of vertical and horizontal fractures on the productivity or conductivity of waterflood and gas injection patterns, but in all cases the fractures initiated at the well and simulated commercial fractures.6,7 It is believed that natural fractures may exist out in the reservoir matrix and would effect the displacement performance. The purpose of this report was to study the effect of a few isolated vertical fractures existing out in the reservoir matrix on the performance of waterflood or gas injection patterns. DESCRIPTION AND EQUIPMENT Since the potentiometric model was described by Lee8 it has been widely used to study numerous types of fluid displacement problems. The potentio- metric model can be used for fluid displacement studies when the following assumptions are valid: (a) steady-state conditions exist; (b) the mobilicy ratio is one; and (c) the capillary and gravitational effects can be neglected. Five-spot and direct line drive square patterns were studied using a copper strip of the desired length and orientation to simulate a vertical fracture in the reservoir matrix. The 20 x 20-in. model was considered to represent only one element of an infinite array of similar patterns. DISCUSSION AND RESULTS FIVE-SPOT PATTERN — ONE FRACTURE BETWEEN INJECTION AND PRODUCING WELLS Fig. 1 shows a quadrant of a five-spot pattem with the vertical fracture existing out in the matrix along a line connecting the injection and producing wells. Length of the fracture was equal to 35.4 per cent of the distance between injection and producing wells. In this particular pattern it is noted that fluid breaks through from the injection well to the producing well when the dimensionless time is equivalent to 26.9. Dimensionless time is defined as volume of fluid injected divided by the volume of displaceable fluid in the pattem expressed as a per cent. The cross-hatched areas of this

Jan 1, 1967

By D. H. Gurinsky, D. G. Schweitzer

The empirical equilibrium constantsd the heat of reaction for the reduction have been determined from 300° to 500°C. The mechanisms of the oxidation of uranium and magnesium from dilute bismuth solutions have been investigated. The kinetics of some of the reactions were affected by air adsorbed on the uranium oxides. THE Liquid Metal Fuel Reactor studied at Brook-haven National Laboratory uses a solution of uranium in bismuth as the fuel. Corrosion and stability problems require that the fuel solution contain corrosion inhibitors and deoxidants. Of the additives tested, magnesium has been found to be the best deoxidant and zirconium the best corrosion inhibitor in bismuth. Some experiments1 indicate that magnesium may play a secondary role in increasing the effectiveness of zirconium in inhibiting corrosion. The work reported here was intended to determine the relative concentrations of Mg and U required to prevent oxidation of the uranium fuel in the event of an air leak. In addition, an attempt was made to identify some of the mechanisms of the oxidation and reduction reactions occurring in the liquid bismuth solutions. EXPERIMENTAL Apparatus and Procedure—The equipment shown in Fig. 1 was designed so that the reaction kinetics could be studied as functions of temperature, melt composition, pressure, and flow rate. At the start of a run, the 10-liter bell jar (F) is removed by opening the bottom Dresser fitting (L), and the bismuth and additives are placed in a pyrex crystallizing dish (I). The bell jar is replaced, the system is evacuated to 10-5 to 10-6 mm Hg and then the heater (J) is turned on. When the selected temperature is reached, the oxidizing gas is admitted to the system to the desired pressure as read on the manometer (D). Gas is introduced through the slide tube (Q) which has a pyrex frit (R) sealed to its bottom. Stirring is accomplished by immersing the end of the tube in the melt. To maintain a constant pressure under flow the input and evacuation rates must be equalized. To obtain this condition stopcock (C) is closed and flowmeters (B) are matched. The reaction is "stopped" by closing the inlet flowmeter and quickly evacuating the system. To obtain a liquid metal sample, the system is evacuated with the sampler positioned near the surface of the liquid. The sampler, Fig. 2, is then immersed until its thermocouple reading matches thermocouple (G) Fig. 1. The system is then pressurized with helium which forces the solution through the frit. Twenty to 30 min elapsed between the time the system was evacuated and sampled. During evacuation and pressurization of the system, the temperature variations did not exceed ± 5°C. While bubbling occurred the temperature remained constant to within 1/2 oC. The volume of the equipment was determined by adding a known volume of air at atmospheric pressure to the evacuated system and noting the resulting pressure.

Jan 1, 1962

By W. A. Burns

Vertical flow is an important mechanism in many petroleum reservoirs. Yet no adequate method has heretofore been proposed for determining the in-situ vertical permeability in a reasonable length of time. Core measurements of vertical permeability do not necessarily represent effective in-situ values. Cores may be representative of only localized pockets of high or low permeability, and total reliance on cores might cause one to overlook the presence of vertical fractures altogether. Conventional single-well tests, interference tests, and pulse tests commonly are used to determine horizontal transmissibility, khh/u, but generally are not designed to measure vertical trans-missibility or other groups containing vertical permeability, ku. A new well test described here provides a measure of in-situ vertical permeability as well as the previously measurable horizontal permeability. This "vertical test" determines the vertical diffusivity, kv/gcu, and horizontal transmissibility, khd1/u (where dl is the length of injection or production perforations). The test is basically a short-term interference or pulse test normally lasting less than a day. Conventional down-hole equipment is used between two isolated sets of perforations in the same well, as shown in Fig. 1. A "layer-cake" description of a reservoir can be obtained by running two or more vertical tests in the same well, each over different vertical intervals. Each test can be analyzed for the effective vertical and horizontal properties of the interval tested. In addition, a vertical well test can determine the effective- ness of a shale barrier or tight zone by testing across it. Field use to date of vertical testing has confirmed its utility in reservoir description. Mathematical Basis Interpretation of vertical tests for ku and kh is based on comparison of measured pressure response with computer-generated response. Pressure behavior during a vertical test in a homogeneous, anisotropic, infinite-acting reservoir (infinite in radial extent and bounded in both directions vertically) is given by the diffusion equation for slightly compressible fluids with negligible gravity effects.l A schematic diagram of the problem is presented in Fig. 2, and an analytic solution, based on the image

Jan 1, 1970

By K. H. Coats, M. H. Terhune

Analysis and example applications have been performed to compare the accuracy and computing speed of alternating-direction explicit and implicit procedures (ADEP and ADIP) in numerical solution of reservoir fluid flow problems. ADIP yields significantly greater accuracy and requires about 60 per cent more computing time than ADEP, not 300 or 500 per cent more as reported elsewhere. l,2 INTRODUCTION Several recent papers 1-3 discuss an alternating-direction explicit difference approximation (ADEP) to the diffusion equation. Example applications of ADEP and ADIP4 were reported to support conclusions that ADEP is comparable in accuracy to ADIP and requires one-fifth to one-third the computing time of ADIP. Applications of ADEP in calculation of two-phase flow in reservoirs was also proposed. 3 This study was performed to compare further the relative merits of ADEP and ADIP in simulation of two-dimensional flow of one and two fluid phases in reservoirs. Since two-phase flow equations are often essentially elliptic rather than parabolic, the efficiency of ADEP in solving the elliptic equation was also examined.

Jan 1, 1967

By C. G. Dunn

James, Jones, and Leak (JJL) conclude that growth-rate and orientation-distribution data obtained in a steep thermal gradient should be used with caution to account for isothermally produced recrystallization textures. I agree. However, it is not clear whether they conclude a) that the orientation dependence of rate of growth of secondaries in (commercial) silicon iron depends on both temperature and a temperature gradient or b) that, in applying rate data obtained in a steep thermal gradient, the isothermal annealing temperature must be close to the temperature for which the rate data apply. Also, it is not clear what tem-

Jan 1, 1967