Search Documents

By Richard Burt

"In seeking to supply the ever increasing needs for raw materials, the mining industry has been required to develop lower grade resources in increasingly remote regions of the world. This has resulted in the extension of many frontiers, both technologically and geographically.The accent of mineral processing research over the last few decades has been to “do existing things better”. There have been many breakthroughs in equipment and processes; however in general the accent has been to increase capacity of existing unit processes, often by orders of magnitude, or to refine the process itself by basic research and development. This has lead to added sophistication, enabling the process engineer to maintain optimum performance of his equipment while reducing the cost, and the amount of labour needed to operate the plant.While much of the increased requirement for resources is handled by the development of ever larger, lower grade orebodies, a significant portion still comes from the small miner - even the artisanal miner - operating in remote and often politically sensitive parts of the world. For them, innovation and simplicity are the keys, not sophistication.The tension between sophistication and simplicity - or research and reality - is specially evident in the area of gravity concentration of rare metals, such as tin and tantalum. This paper looks at this tension, primarily from the eyes of an operator and seeks to put into perspective the two paradigms."

Jan 1, 1999

By Enrique Hauser

Discussion of the paper of ENRIQUE HAUSER, presented at the New York meeting, February, 1916, and printed in Bulletin No. 110, February, 1916, pp. 521 to 534. GEORGE A. BURRELL, Pittsburgh, Pa. Mr. Hauser determined the influence of inert gases on the properties of fire-damp mixtures, and one might take some issue with his work on the effect of hydrogen, on those mixtures. It is valuable information to have, but as far as the practical application of it goes, it is not of much significance, because one does not find it in mines, at least not in any appreciable amount. The same Fig. 1.-APPARATUS FUR DETECTING COMBUSTIBLE GASES IN AIR. thing can be said of ethylene. Ethylene is formed in considerable quantities as a product of heat reaction, but normally, in the coal mines, as a result of a great many analyses we have made at the Bureau of Mines, we cannot say that ethylene is a combustible gas found in the fire-damp in coal mines.

Jan 5, 1916

By Richard Akerman

(Translated by FREDERICK PRIME, JR., Professor of Metallurgy in Lafayette College, Easton, Pa.) [THE attention now being paid both in this country and Europe to the greatest economy in the working of the blast furnace, and the eagerness with which all thoughtful men in the iron business look for any new information which will tend to throw light on this sub¬ject, is a quite sufficient excuse for presenting this paper to the Amer¬ican Institute of Mining Engineers.] I. It must be borne in mind that when two bodies unite to form a chemical compound, heat is almost always evolved. In order to

Jan 1, 1873

By Jerry D. Pfarr, Gary Loving, David G. McDowall

Introduction The Open Cut is one of the best known landmarks in South Dakota's Black Hills. It is the site of the original discovery of the Homestake ledge, or lead, in 1876. Between that time and 1945, the Open Cut produced about 36 Mt (40 million st) of gold ore, the majority of which was drawn from underground operations. In 1981, a drilling program was undertaken to determine how much ore the old timers left behind. The results of this study indicated an ore body made up of broken and caved ground, and solid blocks in pillars and crowns. Reserves are estimated at about 6.8 Mt (7.5 million st) of gold ore, grading at 4 g/t (0.116 oz per st). Total tonnage of material, including waste, is estimated at more than 60 Mt (66 million st). The mining plan is a combination open pit and underground operations. The bottom 78 m (255 ft) will be mined by the underground panel-caving method. The open pit was designed as a contractor operation up to the expected milling stage or the point of sustained ore production. The pit was designed to avoid disturbance of US Highway 85 near Lead, SD. In 1983, the South Dakota Board of Minerals and Environment approved the mining and reclamation plan. Test mining began in the north of the Open Cut using open pit methods. Upon completion of the test pit, Homestake Mining's board of directors approved an expanded mine plan. They committed an additional $27.1 million to prestrip the ore body that is expected to yield 23 t (750,000 oz) of gold. Project life is estimated at 15 years. Exploration program In early 1980, 10 diamond drill holes were drilled from the surface in the vicinity of the east side of the Open Cut. The reserves delineated by that drilling were too small to justify their development. A continuation of this study attempted to quantify the tonnage remaining in the bottom of the Open Cut in the crushed remnants of old stope-and-crown pillars. This was originally thought to be minable by improved drawhole methods using highly mechanized procedures. By early 1981, however, thinking turned toward surface mining methods. The effort to quantify minable ore tonnage ultimately involved drilling 385 reverse-circulation rotary drill holes. These were drilled from August 1981 to March 1983 in several groups of holes. Each group, or drilling program, had a slightly different immediate goal or purpose. The holes were drilled to various depths using truck-mounted Drill Tech 40 drills with doublewalled 76 mm (3 in.) pipe with 134 mm (5.25 in.) three-cone rotary bits. Initially, wet drilling was employed using a cyclone collection system. Since cutting recovery was not as high as it should have been, all subsequent drilling was dry. The cuttings were collected and deposited by way of a cyclone into the sample bag. The pickup for the cuttings using the three-cone bits was about 50 mm (2 in.) above the bottom of the hole. Overall cutting recovery for the drilling program was estimated at 75%. Samples were taken on each 1.5 m (5 ft) of hole drilled. Each sample was split into three 48-g (1-2/3 assay ton) charges. The accepted value was the straight average of the three charges. Because of the known "nugget effect," assays were nearly impossible to duplicate. Therefore, to validate the regular assays, a series of screen fire assays were run. Based on this method of assays and checks, it was concluded that 10% more ounces would be produced than indicated by mine grade. Test pit The northernmost portion of the mine as planned is known as the Terraville test pit. This portion was mined ahead of the rest and separately. It provided answers to questions on slope stability, reserves correlation, mill recoveries, mining costs, and surface mining conditions in various weather and ground conditions since many boulders and much timber were expected. The mining plan for this portion of the pit was successful. More ore was mined than was expected. Mixing of Open Cut ore with underground ore in the mill feed was successful. Impact on the surrounding communities

Jan 10, 1987

Efficient project planning during the detailed feasibility to final engineering stages of a project requires not only the accurate prediction of the grade, tonnes and recovery, but their degree of variability as well. An understanding of the major parameters involved in the grade and tonnage estimates is essential, and the communication to those receiving the results is most important. The personnel involved in the project require more than just a global reserves classification. Quality parameters other than grade should be evaluated as carefully as the estimate itself. This paper does not present any novel ideas to overcome or enhance any of this. Rather, it takes a step back and tries to look at the global picture of requirements for proper project planning.

Jan 1, 1987

By A. G. Pasamehmetoglu, N. Celebi, K. Erarslan

A deposit evaluation system has been developed for a proposed new mine in the basin around Collolar in Turkey. It includes compositing, gridding and reserve estimation facilities and takes account of the economic value of the pit and the pit bottom elevation

Sep 1, 1991

By Charles L. Allie

Since Reserve Mining Co. began commercial operations at Silver Bay, MN, some 25 years ago, major improvements in taconite processing technology have become available. In addition, significant new government standards have been imposed on industrial operations. By implementing a $370-million modernization program, Reserve was able to accomplish two major company objectives: the quality of concentrates produced from the Peter Mitchell mine has improved, and the operation meets both present and foreseeable environmental standards. Process modifications at Reserve were designed to: •Reduce the pellet silica content from 9% to 5%; •Increase the pellet iron content from 62.5% to 65.8%; •Implement a tailings disposal system utilizing "best available technology" water quality controls; and •Reduce air emissions by installing "best available technology" air emission controls. By adding dry cobbing, fine screening, and flotation to the concentrating process, pellet silica content has been reduced; the process reclaims and reuses all process waters; all tailings are stored in a zero discharge tailings basin; and all air emission standards are being met. Initial Beneficiation at Silver Bay Reserve's Peter Mitchell mine at Babbitt, MN, is on the eastern end of the Biwabik iron formation. The open-pit mine is approximately 16 km long, 0.8 km wide, and presently has a maximum depth of about 50 m. Principal constituents of the ore body are magnetite, quartz, and iron-magnesium silicates. The company's large commercial plant began operating at Silver Bay, 75 rail km from the mine, in 1955. Facilities have historically had the capacity to process 31.5 Mt/a of taconite, containing an average 24% magnetic iron, into 11 Mt/a of iron ore pellets. Reserve's initial beneficiation process-which existed up to early 1979-consisted of two crushing stages at the Babbitt mine site, two additional crushing and screening stages at Silver Bay, grinding in a conventional rod-ball¬regrind mill circuit, size classification with cyclones, and concentration with magnetic and hydraulic separators. Concentrate produced by that circuit assayed 65% iron and 8.5% silica, with a screen analysis of 94% -44 µm, and a surface area of 9300 cm2/cm3. Pellets produced from this concentrate assayed 62.5% iron and 9% silica (dry basis). Construction began in 1977 to modify Reserve's operation. Process modifications were designed from pilot scale operations that had been operated and researched by Reserve for more than 10 years. There were two design restraints imposed on the process: economic constraints that Reserve continue to use much of the existing processes and structures in the new operation; and a substantial number of government-imposed permit constraints. In effect, the tailings process was designed first and then a concentrating process was developed that would "fit" the tailings disposal system. Operation of the modified process began in June 1980. Crushing Remains Unchanged No changes in process flow were made at either the Babbitt mine and crushing plant or at the Silver Bay fine crusher. Gyratory crushers reduce run-of-mine taconite to -90 mm-diam at the mine. Taconite is then transported to Silver Bay by intra-plant railroad. There are 10 fine crushing lines. Each consists of a double deck screen, followed by a heavy duty short head crusher, a second double deck screen, and second short head crusher. All screen undersize material bypasses the crushing lines, resulting in an open circuit crushing process. Fine crusher product is approximately 98% - 20 mm, 80% -12 mm, and 50% -10 mm in diameter. Each fine crusher line processes about 400 t/h of crude taconite. Cobbing Provides Aggregate for Tailings Dams The first major process change was the addition of dry cobbing, a dry magnetic separation step between the fine crusher and concentrator. Dry cobbing serves two major purposes in the new process. Its principal function is to provide large volumes of dry coarse aggregate for use in construction of tailings dams. This process rejects about 20% of the fine crusher product as tailings with a size analysis of 98% -20 mm and 50% - 10 mm. These tailings are free draining and have high frictional shear strengths-two key requirements for dam building aggregates. Fine crusher product is conveyed to a new dry cobbing building. Here crushed ore is discharged into one of two 30-minute capacity surge bins. Crushed ore is withdrawn by 12 constant speed conveyor belts; each discharges 500 t of crushed taconite into a dry cobbing machine. Each machine consists of a 3-m-long riffle splitter and two vari-

Jan 12, 1980

By H. T. Kennedy, J. E. Little

An empirical equation for the prediction of the viscosity of several pure paraffin hydrocarbons and nitrogen is presented. It involves temperature, pressure and six constants of the material, and it applies reliably to both liquids and gases. The equation is similar in form to van der Waal's equation of state. For the paraffin hydrocarbons methane through r-hexane and nitrogen, an average absolute deviation of 1.9 percent was obtained on 1,006 data points described in the literature by 14 authors. When this equation is extended to complex, liquid hydrocarbon mixtures, a correlation was obtained with an average absolute deviation of 9.9 percent. INTRODUCTION Equations describing the flow of gas and liquid through porous media contain the viscosity coefficient of the fluid. If other pertinent variables remain constant, the volume rate of flow is inversely proportional to this coefficient. In dealing with condensate fluids and volatile oils, however, the compositional effects resulting from changing pressure materially affect the viscosity. The effect of compositional changes also may be significant in certain secondary recovery or pressure maintenance processes, notably miscible displacement or gas injection. Early attempts to describe the performance of reservoirs utilized a volumetric material balance method wherein gas and liquid in the reservoir were identified as produced gas and liquid at the surface. This method of analysis proved adequate for reservoirs at moderate temperature and pressure that contained gas with moderately low amounts of condensable materials. The volumetric material balance procedures for "black oil" reservoirs leave much to be desired when applied to condensate and volatile oil reservoirs because phase behavior and compositional changes are relatively more important in these cases. The alternative is a compositional material balance, which in turn, requires a correlation of properties of the reservoir fluid with composition. This paper supplies this correlation in regard to viscosity, for reservoir crude oils. REVIEW OF LITERATURE The literature contains many empirical equations describing the effects of composition, temperature and pressure on the viscosities of pure liquids and binary liquid mixtures. However, the applicability of a majority of these equations is limited to very low pressures and to a small number of systems. Most of them, when applied to complex hydrocarbon systems, are of little value. The lack of utility of the majority of equations results from the fact that they were developed to show the separate effect of temperature, pressure or composition on viscosity, but not to predict the viscosity as a function of all three variables. And with the few exceptions noted below, they were developed to apply to much simpler systems than hydrocarbon mixtures. Eakin1 developed an equation for the viscosity of pure light hydrocarbons applicable in the gas, liquid and dense fluid regions, and in its most general form is written as µ - µo =A(ebP - e-cpk) The pure component viscosity at atmospheric pressure and temperature T, µo, was calculated using the Sutherland equation: µo = BT3/2/T + S The constants in Eakin's equation permitted the term e-CPk to become quite small at high densities (above 0.215 gm/cc), and the general equation would simplify to

Jan 1, 1969

By G. E. Perry, R. D. Seba

The steam soak process is the most widely applied and most successful thermal supplemental recovery process in use today. This process, which consists of injection of steam in various quantities into a well that is subsequently used as a producer, can be repeatedly applied to the same well. The treatment, if successful, will increase the production rate of the well for a given period. For this reason, the steam soak process should be considered primarily as an acceleration method rather than as a method for producing oil that could not be recovered by primary means. However, application of this process to reservoirs containing viscous crude may produce oil that otherwise would have zero present value or would be produced at such low rates by primary methods that it would not be profitable. A major necessity for the successful application of the steam soak process is a source of natural energy. This source of energy may exist in one or more of the following forms: expansion of fluids, compaction, gravity drainage. This study was undertaken to determine the effect of gravity drainage on steam soak production and recovery and to develop a method for predicting production rates, ultimate recovery, and heat requirements when applying this process to a reservoir in which the primary driving mechanism is gravity drainage. (For additional discussion of some of the factors that affect recovery of viscous oils by the steam soak process, see Ref. 5.) Model The model developed for analyzing the performance of the steam soak process in gravity drainage reservoirs is a modified form of that developed by Matthews and Lefkovits.1,2 The model used in this study consists of two concentric cylindrical volumes of productive formation with cap and base rock of impermeable material. The inner cylinder, extending from the radius of the wellbore to some predetermined radius, is referred to as the hot zone; and the outer cylinder extending from the exterior boundary of the hot zone to the outer drainage boundary will be referred to as the partially heated zone. Fig. 1 is a schematic cross-section of the model, indicating its thermal properties. Horizontally the reservoir is divided into two regions, the upper region having a high gas saturation and residual liquid saturation, and the lower region containing a high liquid saturation. The boundary between these two regions is referred to as the free surface, which drops and changes shape as fluid is produced from the reservoir. The model was used to analyze the steam soak process by applying the following assumptions: 1. The horizontal formation is composed of one or more homogeneous noninterconnected layers containing incompressible fluids. Each layer is treated separately. 2. The hot zone is heated by the initial soak and is maintained at the steam temperature corresponding to atmospheric pressure by resoaking. 3. The height of the free surface at the producing well is at all times zero (producing well is pumped off). 4. The position of the free surface at the start of

Jan 1, 1970

By B. H. Caudle, L. G. Davies, I. H. Silberberg

This paper presents a method of predicting the recovery and performance of a five-spot steam injection project, in which a realistic approach to pattern sweepout efficiencies is made. Published methods for radial systems were modified for the five-spot pattern by approximating the stream lines with straight lines radiating from the injection well and then converging to the producing well. In each radial segment, the position of the steam front and the temperature profile ahead of the steam front were determined by heat balance equations, which included an estimation of heat losses to surrounding formations. The location of the saturations behind the cold water front was determined from a Buckley-Leverett solution to the material balance equation. Results from this program show steamflood recovery in a five-spot pattern to be considerably less than that predicted for true linear or radial flow systems. For a specific reservoir containing 900-cp oil, a steamflood in a purely radial flow system was predicted to recover more than 75 percent of the original oil in place when 2 PV of water had been injected as steam. A five-spot steamflood with otherwise identical properties was predicted to recover 10 percent of the original oil in place when 0.15 PV of water had been injected as steam and to recover almost no oil thereafter. A cold water five-spot flood in this system was predicted to recover approximately W percent of the oil in place with I PV of water injected. For a five-spot pattern in an example reservoir with 10-cp oil, steam injection similarly showed lower ultimate recovery than water injection but no improvement in recovery rule. Introduction The thermal recovery method considered in this study is steam injection in a five-spot pattern. Pattern steam injection has been given little attention in the past, possibly because of some rather obvious disadvantages. Unless steam temperature is maintained throughout the entire swept area, the process will revert to simple waterflood with all heat being lost prior to reaching the oil bank. Also, the oil viscosity reduction takes place near the steam front and not around the producing wellbore, so that low producing rates must still be endured. This is the reason for the success of the steam stimulation, or cyclic injection, in which the fluids are produced from the same well used for injection. On the other hand, steamflooding can offer some advantages. Oil displacement is by four methods: (1) mechanical displacement by the condensed water, (2) viscosity reduction of the crude oil, (3) swelling of the crude oil, and (4) distillation of the crude oil in the steam zone. Although dependent upon the crude, laboratory experiments have shown displacements up to 80 percent by steamflooding. In addition, steam is a good heat transport medium, since it is cheap and has a high heat content. Previous Investigations Previous publications have reported methods for predicting recovery in a steamflood for linear and radial systems. Since no pattern sweep efficiency is taken into consideration, the recovery even from a radial system must be greater than from more realistic geometries such as the very common five-spot patterns. Probably the broadest coverage is given by Willman who reported experimental results and offered a method of predicting recovery for a radial flow system. They concluded that both hot water and steam injection recover more oil than an ordinary waterflood, and that steam injection could yield recoveries "as much as 100 percent greater than by water flood." They also concluded that both the heat requirements for a reservoir and the residual oil remaining in the reservoir after steam injection were independent of the amount of oil originally in place, that short exploitation times were desirable, and that a high percentage of net sand in the reservoir with a high initial oil saturation was desirable. The method assumes that the flood occurs in three concentric cylindrical zones: (1) an inner steam zone, (2) a central hot water zone, and (3) an outer cold waterflood zone. Displacement in the steam zone is based on laboratory-determined residual oil saturations while the hot water and cold water zones use the conventional Buckley-Leverett equations. Although the results shown by Willman et al. are for a radial system flowing out from the well to an assumed external circular boundary, their equations did

Jan 1, 1969

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

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

Jan 1, 1967

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

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

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

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

Jan 1, 1970

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

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

Jan 1, 1970

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

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

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

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

By A. S. Odeh

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

Jan 1, 1970

By 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 H. S. Price, J. C. Cavendish, R. S. Varga

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

Jan 1, 1969

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

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

Jan 1, 1967