Search Documents

By G. H. Glos

There was a time when an occasional failure of a tailings embankment was considered an acceptable risk associated with mining activities. This is not the case anymore. A major tailings failure today could result in the indictment of mine officials, huge fines and a possible site closure due to the loss of operating permits. As a result of this changed operational environment, tailings disposal facilities must be operated with these potential consequences in mind. Physical controls for tailings embankments are prohibitively expensive. However, the nuclear industry has developed a system of administrative controls as an alternative to physical controls. Administrative controls may be defined as control of actions, through organizations and management, using written procedures, record keeping, review and audit to assure operation of a facility in a safe manner (10 CFR 50). Administrative controls can reduce the cost of operating tailings impoundments, as well as increase the level of confidence in these facilities.

Jan 1, 1999

By E J. Hofmann, J E. Hofmann, A J. Hofmann-Fritz

"With its Head Quarters in Perth, global heavy engineering company Hofmann Engineering is one of the largest family owned engineering companies in Australia. The company has pioneered a number of new manufacturing processes and design improvements for components of large comminution equipment.Some of these innovations are: Innovative design improvements of high pressure grinding roll (HPGR) forged steel tyres. Hofmann Engineering has developed unique patented tungsten carbide stud and edge protection designs which improve wear life and introduce quick and cost effective replacement options to further extend the life of the tyre. Development of forged fabricated mill heads using rolled and welded high tensile steel plate as an alternative to cast conical mill heads. Establishment of Australian capability to manufacture and export up to 42 ft. (12.8 m) diameter SAG mill shells. Techniques for fast track, product improved replacements of cracked or damaged mill gears, including in situ profile & wear measurements of girth gears using portable coordinate measuring machine (CMM) technology, forged rim steel and induction hardening technology. Gear and pinion mechanical drive designs to 24 MW or higher for SAG mills of 38 ft. to 42 ft. (11.6 m - 12.8 m) diameter as reliable alternatives to gearless drives.The paper presents case studies from the Australian heavy manufacturing facilities of Hofmann Engineering on large, high powered equipment for both the domestic and export markets.CITATION:Hofmann, E J, Hofmann, J E and Hofmann-Fritz, A J, 2013. Australian innovations in high pressure grinding rolls and large grinding mill manufacture, in Proceedings MetPlant 2013, pp 141-151 (The Australasian Institute of Mining and Metallurgy: Melbourne)."

Jul 15, 2013

By M. A. Trevits, T. A. Gray, L. M. Crayne, P. Glogowski

Mine seals can he remotely constructed in underground coal mines through vertical boreholes when direct access to a fire is impossible or considered to he too dangerous. This method has great merit because the boreholes can be drilled to specific mine areas, and the seals can be positioned close to a fire zone. The technology however can be largely ineffective if the constructed seals do not provide of effective barriers to airflow or if they cannot he used to impound water and other inert materials. Unfortunately, no viable alternatives exists to sealing the entire mine to the ground surface. Full-scale remote mine seal construction research is being conducted at the National Institute for Occupational Safety and Health (NIOSH) Lake Lynn Laboratory (LLL). Under this effort. Howard Concrete Pumping Company and GAI Consultants Inc. have joined forces to evaluate a potentially significant improvement to the current state of the art. This paper covers the development o f novel technology or rewrote mine seal construction, the evaluation of rite materials used, the construction practice and the Fellow-up testing.

By J. T. Nel

The establishment of a pilot plant for the continuous production of titania pigment and nano-titania is discussed. The plant was designed, constructed and operated on the Necsa site. The TiO2 production capacity is 7?10 kg.h-1 by the reaction of TiCl4 and O2at>1000°C. Liquid TiCl4 is sprayed under pressure into a plasma reactor where it evaporates and reacts with O2 to form TiO2. Process optimization has allowed for continuous production without any clogging or blockages. The primary advantage of the plasma-assisted process is the possibility to recover and recycle the chlorine values in high yield as Cl2 for use in the manufacture of TiCl4. The HCl formed in most conventional fossil-fuelled processes makes Cl2 recovery and recycling much more difficult. Particle size distribution analysis of the fine, white product revealed a d50= 125 nm. The raw product consisted of 60% rutile that could be increased to 80?100% rutile by post-annealing. The specific surface area (SSA) of the pigment was 8?20 m2?g-1(BET). According to the CIE colour index the ?whiteness? (L*>95) of the plasma-produced pigment was better than that of several commercially available pigments. The particle size could be manipulated by controlling the TiCl4 feed rate, gas flowrates, and the particle quench rate and residence time. Nano-sized TiO2 powders with a particles size (d50) as small as 50 nm were also successfully produced in this way.

Jan 1, 2010

By G. L. Oldwright, F. K. Shelton, W. E. Keck

"INTRODUCTION Zinc must be of high grade to satisfy the present demand. Even when the amount of premium paid is small, the purer metal as the buyers' preference. The following specifications for spelter were taken from the American Society for Testing Materials for 1927, Part I, pane 524. (The specifications for 1930, p. 582, are tile same.)Another specification for zinc of a much higher grade than any of the above has been established lately by the same society. The committee on zinc die-casting (Proc. Am. Soc. Test. Mat., vol. 31, 1931, p. 263) reached the following, conclusions: I. ""That the purest of zinc is necessary to produce stable zinc-base die castings, and as a result requested specifications for zinc spelter .with total impurities not over 0.02 percent."" Their tests showed that the die castings made of the purest zinc had the least elongation and the greatest tensile strength, and exhibited the least changes with temperature.Still more far-reaching claims are made for the properties of zinc containing over 99.99 percent metal. Zinc foil, for instance, may someday be made as a common commercial product. It would seem likely that the demand for zinc of high grades rill continue and that even purer varieties may be desired. This demand fits in with the modern trend in manufacturing plants of reducing the number of variables to a minimum by using pure materials, and of controlling the processing of metals during their various stages of fabrication by mechanical means instead of by expert laborers."

Jun 1, 1933

By Ludovic Charpentier

Production of metallic wires directly from liquid jets would avoid drawbacks of extrusion (space needed, damaging, cost?) and enable recycling. But liquid jets tend to break, due to capillary instabilities, unless a solid layer (coating or start of solidification) surrounds it. We recently developed a process to produce tin wires, cooled by water flow parallel to the jet. Superficial oxidization of tin produces a glass coating around the jet. Studies showed the length of wires increased when increasing speed of tin jets, and that water average speed should be the closest possible to the jet one to have smooth fibers. Producing 1 mm diameter wires is possible. The length of wires is limited by the experimental setup, which does not provide possibilities to coil the wire. Most of other metals do not generate glass-like oxides. The lack of coatings makes it necessary to reduce the capillary instabilities to obtain fibers.

Jan 1, 2006

The Canada Centre for Remote Sensing (CCRS) has developed an airborne sensor, the SWIR Full Spectrum Imager (SFSI) which detects in the short wave infrared range of the solar spectrum from 1220 to 2420 nm. In June of 1995, a demonstration flight was undertaken in the state of Nevada with the support of the companies listed above. The focus of the flight was to collect image data over well known mineral deposits with the objective of evaluating the applicability of this sensor to alteration mineral mapping. Some of the Nevada sites chosen include Cuprite and Goldfield, volcanic hosted acid sulfate systems; Yerington, a structurally complex porphyry copper deposit; and Virginia City, a quartz-sericite-pyrite vein system with advanced argillic pods. Ground spectral data was collected with the PIMA-II spectrometer. Solar reflectance calibration was done by Desert Research Institute personnel at Virginia City and Cuprite. This paper will discuss the configuration of the sensor and present image and ground data for Goldfield and Virginia City.

Jan 1, 1995

By Timothy P. Smirnoff

Sustaining the region?s mobility needs, achieving projected economic growth, and maintain Midtown Manhattan as a center of regional, national and global importance has led to continued investment in commuter rail and transit infrastructure in New York City and resulted in design and construction of large underground caverns. The current plans for stations in Manhattan include New York Penn Station Expansion (NYPSE) caverns for Trans Hudson Express (THE) Tunnel Project for New Jersey Transit (NJ Transit) and the Port Authority of New York and New Jersey (PANYNJ), Metropolitan Transportation Authority Capital Construction Long Island Railroad (MTACC/LIRR) twin caverns for East Side Access (ESA) Project, and station caverns of the No. 7 Line Subway Extension Project and Second Avenue Subway Project for Metropolitan Transportation Authority Capital Construction/New York City Transit (MTACC/NYCT). This paper provides a general overview of these large caverns now in design or construction in Manhattan, the area?s geologic setting, the geotechnical factors considered and approaches commonly applied for cavern design.

By Clarke P. E

At the end of the recent war, reports became available from Europe of the use of tungsten-carbide inserts in percussion rockdrill bits. These bits were used with small size rockdrills, so that the only advantage gained from the use of the harder, abrasion resistant, tungsten-carbide appears to have been the greater life of the bit. Larger machines, such as used in Australia in conventional practice, had to be discarded to protect the relatively brittle tungstencarbide inserts and the vulnerable bonding with the bit body.As it had been found in machine shop practice that some grades of tungsten-carbide had been produced with considerably improved toughness, it was decided by Gold Mines of Australia Ltd. that drilling practice with tungstencarbide bits warranted detailed examination under Australian conditions, and that drilling tests should be undertaken at Morning Star (G.M.A.) Mines N.L., for which it acts as General Manager and Consultant.These investigations progressed through three stages:1. Preliminary Investigations and Tests.When the investigations were initiated no commercial bits were available in Australia, and manufacture was commenced with improvised equipment from virtual scratch. These early tests indicated that the use of 3 in. bore rockdrills in conjunction with tungsten-carbide insert bits was a 'possibility, and that it was unlikely that a small bore...

Jan 1, 1950

A study of literature on the epidiorite dykes which transgress the Broken Hill lode and its wall rocks reveals a considerable diversity of thought regarding the time of emplacement of such dykes. In this paper it is deduced that the epidiorite dykes were formed at a late stage and were probably post-ore. Detailed mapping, sampling, and correlation of data from over 100 exposures substantiates this deduction. Physical disruption of the epidiorite dykes is indicated and some interaction of the epidiorite and ore at the time of intrusion is inferred.Further petrological study may confirm this.

Jan 1, 1968

By Jochen Franke, Derek D. Lichti

"This paper gives a brief introduction and presents recent enhancements in the MillMapper technology developed for thickness gauging and condition monitoring of liners in grinding mills used for mineral comminution. It also introduces the recent expansion of MillMapper functionality aimed at mill operational optimization. MillMapper uses calibrated laser scanner data captured inside a mill to determine thickness at several million points on all liner surfaces at feed, shell, and discharge ends with millimeter accuracy. Data capture takes nominally 5 minutes, hence interruption to mill operation is negligible.The paper is structured as follows. First, a review of the technology and discussion of recentlyimplemented innovations such as liner change out forecasting, automatic ball charge volume and weight measurement, and automatic mill liner volume/weight measurement for mill operational control or load cell calibration are presented. Combined weight information provides control over load related issues such as load fill levels for relining, or start-up load after reline. Feature descriptions are followed by a listing of the numerous advantages over existing techniques and the benefits to operations. Results from the application of MillMapper are presented throughout.INTRODUCTIONThe need for thickness gauging of liners of grinding mills used for mineral comminution is well known but the task is rarely completed to a satisfactory standard. Liner failure resulting from excessive wear can have expensive consequences. Furthermore, mill performance as well as liner wear are known to be correlated to lifter bar geometry (Napier-Munn et al., 1999; Rajamanni et al., 2005). At present, the following techniques are routinely used to gauge the thickness of mill liners:1. Visual inspection;2. Ultrasonic thickness gauging (for steel liners); and3. Physical inspection (e.g. nail inserted in rubber liners)"

Jan 1, 2008

Entry of personnel under a vertical raise bored hole is not standard practice in Australia. Where access is required below a rise, mining operations must be confident that they understand the capacity of the ground control specified and installed. Physical modelling was conducted on Callie mine to evaluate the performance of 5 mm steel wire weldmesh secured to the base of a vertical raise bored hole of 1.4 m diameter. Joint orientation data was used with Unwedge to calculate a series of potential wedges likely to form within the rise. The mass of the most likely wedges were then replicated by collecting rocks in 5 kg increments (¦10 per cent), which were allowed to freefall onto installed mesh at the base of a 40 m vertical rise. The results indicate that two sheets of mesh can retain a freefalling rock of at least 20 kg mass although rupture of strands can result. This paper outlines the test conditions, expected results based on theoretical calculations and actual outcomes of the testing.

Nov 25, 2010

The method of investigation was to study the distribution of ore in great detail, to make a thorough structural study of the lode and the mine area, and then to compare the results of the two studies. Such a comparison usually reveals a very close relationship between geological structure and ore distribution, and can often lead to important conclusions which would not be reached in the ordinary working of a mine.The Triton mine lies 30 miles E.N.E. of Cue, or 380 miles N.W. of Kalgoorlie, in typical pre-Cambrian rocks, consisting of gTeenstone lavas and tuft's, interbedded with slates and some sediments, and intruded by granite and porphyry. Fig. 1 (at end of paper) is a longitudinal section of the mine and shows that the workings consist of the main North Emu and the smaller South Emu, separated by a distance of about 1 of a mile. Both orebodies occur on the same lode, which strikes north and south and dips practically vertically. On the same line, and of 1/2 a mile and a mile respectively to the north of the North Emu, are the small Rand and North Rand workings, where a big width of lode material has proved to be very low grade. In the North Emu the lode width averages 12 to 15 feet; it is seldom less than 5 feet and in places reaches 30 feet. In the South Emu the lode is wide, but only portion of it is payable.The lode is a mineralised bed occurring between massive lava flows on the east and a mixed series of tuft's and thin lava flows on the west. The lode material is hard, fine grained, and dark coloured, and generally has a faint banding which is believed to indicate a tuft'aceous origin of the original lode rock. Microscopic examination shows that the original lode rock has undergone considerable albitisation;' there has also been considerable formation of biotite which gives the lode material a characteristic sheen. The mineralisation consists of abundant fine pyrite, which is auriferous, and a lesser development of pyrrhotite which is, however, barren of gold.THE OREBODIES AND THEIR CHARACTERISTICSThe distribution of ore is best studied from gold contours, as these reveal features that cannot be detected in an ordinary assay section. Gold contours (see Fig. 3) are developed' by plotting the width and grade of all development and stoping assays throughout the mine on a longitudinal projection, representing each sample by a coloured dot according to a suitably chosen range of assay values, and then outlining areas of similarly coloured dots (see H. J. C. Conolly, Econ. Geol., April, 1936). Both width and grade must be taken into account and the simplest method of gold contouring is to use dwt.lft. values, but difficulties arise here with big widths of low grade and small widths of high grade ore. To choose an extreme case, 30 ft. of 2 dwts. and 2 ft. of 30 dwts. both equal 60 dwts.lft., and would therefore be given the same value in dwt.lft. contouring...

Jan 1, 1943

By D. S. Dihalu

A generalization of Gy?s model for the fundamental sampling error introduced the new ?parameter for the dependent selection of particles?, denoted as Cij. This allows for modeling deviations from the ideal situation where the selection of a pair of particles is composed of two independent selections. The generalized model potentially leads to more accurate variance estimates in the case of clustering of particles, differences in densities or sizes of the particles or repulsive inter-particle forces. A straightforward and practically applicable method is needed for the determination of this parameter for miscellaneous mixtures in industrial settings. In this contribution, the feasibility of using digital image analysis to determine this parameter Cij, will be demonstrated. Line transect sampling of a digital image was used to construct a transition probability matrix. A new algorithm to derive quantitative estimates for Cij will be presented and discussed. The applicability was verified by a photograph of zirconium silicate particles of sizes typical for industries dealing with pharmaceutical, food/feed, and environmental applications. Conditions affecting the practical applicability are identified and potential pitfalls will be discussed, including e.g. how a potential unrepresentative surface can affect the quality of the estimate of Cij.

Jan 1, 2010

By Frederick W. Slack

This paper will discuss the inherent value engineering associated with the design and construction of approximately 184,000 SF of anchored tieback wall, installed as part of the widening and drainage improvements for KY9 in Campbell County, Kentucky. The total project is 2.2 miles in length. The anchored walls, excavation grading and drainage total to a complete project cost of $28.4 million. This is one of the largest design-build projects undertaken by the Commonwealth of Kentucky. It also features another first for the Transportation Cabinet - the general contractor is responsible for quality control/quality assurance. The project is located just south of Cincinnati, OH and is situated along the east bank of the Licking River. The slopes along this stretch of the river have had a long history of movement and instability. The Kentucky Transportation Cabinet advertised the anchored wall portion of this project to be designed and constructed by pre-approved specialty contractors. The Transportation Cabinet provided the basic design parameters to each specialty contractor. Aspects to be discussed include the site and project history, previous attempts at slope stabilization along existing KY9, the technical proposal and bidding process, and the actual design and construction.

Jan 1, 2003

By Beardsmore G. R

Geothermal gradients in the Gippsland Basin commonly exceed 35¦ C/km. Consequently models of thermal maturation may indicate exploration targets with an anomalous age/depth distribution. However regional geothermal models are constrained primarily by observations of surface heat flow. High geothermal gradients may simply reflect the presence of low conductivity coal beds within the basin. Heat flow data are not readily available for the Gippsland Basin. Normally special purpose drill holes are required to provide direct determinations of surface heat flow. These provide core samples for thermal conductivity and allow measurements of a stable geothermal gradient at least one year after drilling. Different indirect procedures are now required to provide estimates for the Gippsland Basin. Bottom hole temperatures can be used with suitable corrections to provide an indication of the regional geothermal gradient. However the corresponding thermal conductivities are more difficult to determine. Geophysical logs can be used to detect the major lithologies but detailed statistical analysis is required for numerical estimates. Complex mixing models are now suggested for geothermal surveys linking physical properties to variations in the mineralogy of each well. Allowing for low conductivity coal measures the data in the Gippsland Basin suggest values for heat flow consistent with regional trends in Australia. Interpretation is complicated by overprinting of transient thermal events with different time scales. Apart from the tectonics of basin evolution geothermal signatures may be complicated by regional ground water migration, the 'hot spot' trace associated with southward younging of Cainozoic volcanism, diagenesis, and basin refraction.

Jan 1, 1992

By T. Wettainen

Malmberget town is located in northern Sweden, approximately 70 km north of the Arctic Circle. Parts of the town overlie more than 20 iron orebodies, consisting mainly of magnetite with smaller quantities of haematite. The mine is operated by the mining company LKAB. Mining started in the 17th century, but not until the railway to the coastal city of Luleå was completed in 1888 did large-scale production commence. Around 1920, mining proceeded underground and today sublevel caving is the only mining method used. Sublevel caving causes subsidence of the ground surface, and buildings and residential areas have been relocated due to the mining activities for more than 50 years. The number of seismic events accompanied by strong ground vibrations is now increasing. In 2008 the mine received a permit from the Environmental Court of Sweden to increase production to 20 Mt of crude ore per year. A prerequisite for the permit was that the mine conducts a number of investigations regarding the environmental impact on the residents of Malmberget. One of these investigations concerned how seismicity will change as production increases and what measures could be taken to reduce inconvenience to the town residents. Today the mine possesses an extensive seismic monitoring system with more than 180 underground and surface geophones. For this study, eleven seismically active volumes in Malmberget mine were identified, and for each of them, a yearly future maximum magnitude interval was estimated based on the current production plan. Relationships between historical seismic events and measured ground vibrations in the town of Malmberget were established, and future ground vibrations caused by expected seismic events were estimated using a probabilistic approach. The outcome was the number of intervals of expected ground vibration per year and per monitoring point. Possible measures to reduce inconvenience for the town residents include blast restrictions, sequencing, and possibly preconditioning. The ultimate long-term solution is an almost complete relocation of Malmberget town. This process has recently been formalized and LKAB is taking an active part in realizing this goal.

Jan 1, 2014

By A. B. Cook, C. J. Walker, G. B. Spencer

Although water will displace oil from a petroleum reservoir to a greater extent than gas will, there are some reservoirs in which gas rather than water should be used for pressure maintenance. This is indicated by the high-percentage oil recovery in the case history reports on the Pickton field1 in East Texas and the Raleigh field2 in Mississippi. In gas cycling, total oil recovery includes vaporized products from the immobile oil in addition to oil produced by displacement. If a large part of the immobile oil is vaporized, total oil recovery may be higher than that obtainable from pressure maintenance by water injection. However, the amount of immobile oil vaporized may range from almost 0 to 100 percent, and no simple and reliable method has been presented in the literature for calculating oil vaporization. Development of such a method was the purpose of this study. Several papers concerning the calculation of oil vaporization have been published,3 6 and each is based on the concept dealing with K values — equilibrium constants for the various components in reservoir oils. (K for a component in a system of vapor and liquid in equilibrium is the mole fraction in the vapor phase divided by the mole fraction in the liquid phase.) This appears to be the most logical approach. Yet, the problem is too great for a perfect solution, even with modem laboratories to analyze oils and high-speed digital computers to perform calculations; the number of components in reservoir oils is too large. Because of their complexity no two reservoir oils are exactly alike. Also, K values change not only with variations in pressure and temperature but also with composition of the reservoir oil. In addition, during gas cycling, the lighter hydrocarbons tend to vaporize first. Thus, the reservoir oil becomes more dense and less volatile as gas cycling continues. Furthermore, the greatest amount of oil vaporization occurs near the injection well. Therefore, a simplified method is required for calculating oil vaporization because a rigorous method is not practical. The present methods of calculating vaporization generally require knowledge of the reservoir oil composition in which the mole fractions of the lighter components through hexanes are given; the remaining oil is described as C, + (heptanes and others of greater molecular weights). Assigning one K value for the C,+ system provides the simplest solution, but this incorrectly assumes that K for C, + does not change as gas cycling continues. This assumption can cause large errors because the K value may become less than a thousandth of the original after a large amount of oil vaporization. Another way to solve the problem is to determine changing K values for C, + according to how much of the C, + has been vaporized. To determine these values a sample of the reservoir oil is injected into a pressure-volume cell, and dry gas is batched in and out of the cell. Appropriate measurements and tests are made to determine K

Jan 1, 1970

By W. G. Payne, A. C. Richardson, B. L. Gandrud

"INTRODUCTION The Connellsville mine of the Yolande Coal & Coke Co. is in the southwestern part of Jefferson County, Ala., about 30 miles southwest of Birmingham. The mine is on the Blue Creek bed, which is about 6 feet thick and dips at a angle of 18º. The bed is worked by the room-and-pillar method and is mined by both machine and hand methods. The Blue Creek bed at this point yields a high ash, low-volatile coal with excellent coking properties. Because of its exceptionally low volatile content and high coking strength, this coal is excellent for blending with Alabama coals of higher volatile content and weaker coking properties to form blends suitable for production of metallurgical coke.A sample of washed coal taken during this investigation gave the following analysis on an as received"" basis: Moisture, 8.3 per cent; volatile matter-, 20.2 per cent; ash, 8.0 per cent; fixed carbon, 63.5 per cent. The high ash content of this coal limits the quantity which can be utilized when low-ash coke is required. In the Birmingham district it is usually difficult to sell as coking coal a product that contains more than 8.0 per cent ash. To produce a. coal with this ash content at Connellsville presents a somewhat difficult washing problem. This paper gives the results of an investigation that was conducted by the Southern Experiment Station of the U. S. Bureau of Mines in cooperation with the University of Alabama and the Yolane Coal & Coke Co. for the purpose of determining the washing characteristics of the coal and the possibility of reducing its ash content on coal-washing tables.The practice at the Connellsville washery at the time of beginning this investigation was to crush the run-of-mine coal to pass a 1 1/2-inch round-hole screen and to treat the satire output on two primary, 3-compartment jigs. A. middling product was recovered from the primary jigs which was re-treated on a single-compartment rougher jig. This system yields a product of acceptable ash content, but the yield is low and the loss of good coal in the refuse is high."

Feb 1, 1933

By E. H. Herron, J. H. Perry

Mathematical simulation of reservoir behavior may be used to help understand reservoir processes and to predict reservoir behavior, thereby leading to the most economically desirable form of exploitation. In addition, simulation can be used as a tool for reservoir description to learn more about the physical nature of the reservoir and the mode of primary recovery. This use is essential in most reservoir studies and represents one of the more significant applications of simulation. Prediction of reservoir behavior provides information concerning displacement efficiency, optimum well locations, and the comparison of alternative producing processes. Since oil, water and gas all commonly occur in many reservoirs, a simulation that takes into account the simultaneous flow of three immiscible phases is a prerequisite for obtaining such information. Furthermore, a simulation that describes this flow in two dimensions provides the capability of evaluating combination gas- and water-drive reservoirs with either an areal or a cross-sectional description. Many reservoir problems require such a simulation. The objective of the research described here was to develop a mathematical model and associated computer program for accurate, efficient, and economical prediction of the reservoir flow of three phases in two-dimensional geometry. Several authors have discussed multidimensional reservoir simulation of two phases.l-3 Gottfried et al.' presented an analy- sis of three-phase simulation in one dimension. Fagin and Stewart" and Garrett6 have discussed the simulation of three-phase flow in two dimensions, but their solution techniques were different from the method discussed in this paper. In particular, previous models and computer programs have usually involved either fewer capabilities or inefficient and less rigorous descriptions resulting in limited applicability. In the following sections, we present a description of the model, an evaluation of the method, and the application to reservoir problems. Additional and supporting material is presented in the Appendix. The Three-Phase Model General Description The three-phase model consists of the mathematical description of the flow of oil, water and gas in a reservoir. Conceptually this description considers flow in all three space dimensions, but in the subject simulator it is assumed that no flow occurs in the third dimension. The model includes the effects on reservoir behavior of fluid and rock compressibility, viscosity, gravity, capillary pressure, relative permeability and gas solubility. Within the two-dimensional restriction, the reservoir considered may be completely heterogeneous and anisotropic. The differential equations governing multiphase, multidimensional flow in a reservoir are too complex to be solved analytically, so they are approximated by a system of finite difference equations and solved

Jan 1, 1970