Search Documents

By Ian Wark

IN an earlier paper1 measurements of contact angles due to the effect of xanthates on mineral surfaces were reported. The solutions in which these measurements were made differed widely from those of modern flotation circuits. It is intended, however, by introducing, one at a time, the other principal reagents of flotation, to approach gradually the solu-tion composition of flotation liquors. The work now described is an attempt to determine the effect of one of the most important of these reagents-namely, sodium cyanide-on the establishment of air-mineral contact. The work has been carried out for the following companies: Broken Hill South, North Broken Hill, Zinc Corporation, Electro-lytic Zinc of Australasia, Mount Lyell Mining and Railway, and Burma Corporation. Since the effect of cyanide is inseparably connected with the alkalinity, the results cited below deal with the mutual effects of alkalis and cyanide. Alkalis alone exert a significant effect on the air-mineral contact but comparatively large amounts are required. Since the addition of a small amount of cyanide greatly decreases the degree of alkalinity necessary for the suppression of contact at surfaces of any of the minerals experimented with (except galena) and since alkalis are rarely used alone in practice, little attention has been paid to cyanide-free solutions. It was shown in the previous paper that, excluding the dixanthogen effect, whenever contact was practicable at all in an ethyl xanthate solution the angle of contact was 60°. There is, therefore, less necessity to measure equilibrium contact angles; in fact, the experiments herein described have generally been conducted with a view to determining whether contact is possible and not to measuring the actual value of the contact angle. The separations which have generally to be effected in dealing with complex sulfide ores are: (1) galena from sphalerite, pyrite, etc.; (2) sphalerite from pyrite and/or gangue minerals; (3) copper minerals from pyrite and gangue minerals; (4) pyrite from gangue minerals; (5) chal-

Jan 1, 1933

By Madgwick J. C

Leaching column techniques and microbiological methods have been described in detail in Madgwick and Ralph (1980) and Madgwick, Ralph and Sadler (1981).Sixty-nine kilograms of copper porphyry ore (0.4 per cent Cu, 2.7 per cent Fe) (particle size 2.4-11.2 mm; 60 The effect of sodium chloride (O.6M) irrigation on bacterial leaching of Bougainville copper porphyry ore was studied.Viability of suspended Thiobacillus ferrooxidants was reduced to less than 10 cells/ml in NaCI (O.6M) compared with 10' eells/ml in dilute II,SO, (0.01 M, pll 2). Salt irrigation had no significant effect on heterotrophic and sulphur oxidizing bacteria. Increascd short term recovery of copper (145 g/tonne/ 15 days) was followed by a marked decrease in tenors to less than 20 mg/litre/day compared with an acid leach, steady release, of 110 mg/litre/day. Resumption of acid leaching restored average solution copper to 46 mg/litre/day and the viability of suspended T Ierrooxidates to 10 cells/ml.Quantitative removal of sodium ehloride from ore particles was difficult and t races of NaCI remained for more than 200 days.

Jan 1, 1982

By W. Lee Daniels

Significant deposits (> 4,000 ha) of heavy mineral sands were discovered in the USA Coastal Plain of Virginia and North Carolina in 1989. The majority of these lands support highly productive row crop agriculture, and the development of restoration protocols that would return these lands to agricultural use was deemed to be critical to the long term sustainability of mining operations. Virginia Tech worked closely with all stakeholders to develop appropriate restoration protocols and to coordinate their implementation. Full-scale mining operations at the Old Hickory Project in Virginia were initiated in 1997 by Iluka Resources Inc., and restoration protocols have continued to evolve in response to a variety of economic, technical, and social issues. Approximately 750 ha of land are in various phases of backfilling and final reclamation. Return of these lands to agricultural row-crop production has been complicated by lateral variability in mine soil physical conditions, excessive compaction, and limited topsoil return. However, within the past two years, a series of tailings deposition and soil reconstruction practices have been developed and implemented that significantly improve post-mining soil productivity. Deep ripping and appropriate use of organic soil amendments have been particularly effective at restoring post-mining soil productivity.

Jan 1, 2005

By J. W. Bridge

Triaxial compression tests at various confining pressures and temperatures, as well as tangent modulus tests using ultrasonic waves, were conducted for wax-coated granular composite materials. The material used is a surface for Thoroughbred horse racing. The purpose of these tests is to examine how track shear strength and tangent modulus at a range of loads respond to changing temperatures. Previous work has shown that a wide range of operational temperatures are experienced which have been reproduced in the lab. Using differential scanning calorimetry, it is confirmed that these operational temperatures correspond to distinct thermal transition regions for the wax that is used to coat the sand in these surfaces. Preliminary results show that both shear strength and tangent modulus are sensitive to temperature changes with maximum strength and modulus correlating with major DSC wax transition regions.

Jan 1, 2010

By Irven A. Jacobson

This literature survey consists of references dealing with gaseous state thermal decomposition reactions of organic nitrogen and sulfur compounds. The reference source for the survey was the Chemical Abstracts from 1930 to 1956, inclusive. Original articles were not studied. The starting date was set at 1930 because previous to this time much of the work is questionable. This report is in three parts: Part I lists nitrogen compounds; Part II, sulfur compounds; and Part III, compounds containing both nitrogen and sulfur. The organization of the classes of compounds is alphabetical: Part I - amides, amines, amino acids, ammonia, azides, azines, azo compounds, heterocyclic nltrogen compounds, hydrazine, nitrates, nitriles, nitrites, nitro compounds and miscellaneous nitrogen compounds; Part II - disulfides, sulfides (including double and triple sulfides), thiols, thiophenes and xanthates; Part III - thioacid amides, thiazoles, thiocyanates and thiuram disulfides. A class formula is included with each class name. The R's used in these class formulas may refer to hydrogen or to the same or different radicals when more than one R appears.

Jan 1, 1960

By Werner Ruegg

El presente texto describe algunos aspectos y propiedades de índole mineralógico, genético e individual de los diferentes yacimientos fosfatados. Existe una diversidad de fosfatos combinados como calcio, fluor, hierro, etc pero están aquellos también que son usados para fertilizantes limitados a los fosfatos cálcicos. De éstos, la apatita es su forma cristalina primaria, constituyéndose como la fuente de combinaciones. La fosforita, a su vez, es su contraparte amorfa, que establece un precipitado químico o un residuo orgánico donde existen combinaciones entre los dos.

Mar 1, 1984

By M. Copeland

Phase relations of gadolinium alloyed, up to 40 weight-percent (wt pct), with AISI 304-type stainless steel were established. Body-and face-centered cubic iron-rich solid solution phases and several intermetallic gadolinium-containing compounds of variable composition that approximated the formulas Fe9Gd, Ni7Gd2, Ni3Gd, Fe2Gd, NiGd, and NiGd3 were identified. The intermetallic phases Fe9Gd and Ni7Gd2 were noted to melt at about 1,080° C at about 0 to 2 wt pct gadolinium. The melting temperature of Fe9Gd increased with increas-ing gadolinium content to a plateau at about 1,230°C and 30 wt pct gadolinium, and the melting point of Ni7Gd2 decreased to a plateau at 900° C and 30 wt pct gadolinium where Ni3Gd, NiGd, and Fe2Gd coexist. Body-centered cubic iron solid solution was stable from 1 to 2 wt pct gadolinium to higher contents, 30 wt pct gadolinium or possibly more, and face-centered cubic iron up to 12 wt pct gadolinium below 1,200° C. The stability of these phases above 1,200° C was not resolved except for the mel points. A limited number of fabrication variables as well as mechanical and corrosion properties of alloys containing up to 5 wt pct gadolinium were studied. Equilibrating and forming operations were best conducted at about 940° to 1,080° C. Some variability in the yield and tensile strengths of alloys was noted; however, there was a continual decrease in ductility and impact resistance with increasing gadolinium contents. The resistance of gadolinium alloys to corrosion in water at 354° C was about the same as stainless steel.

Jan 1, 1965

By J. J. Fritts, J. L. Patrick, T. L. Wright, C. O. Ensign, W. S. White, J. W. Trammell, J. C. Wright, D. J. Hathaway, R. J. Leone

The copper deposit at White Pine, Michigan, from which a little more than 5 per cent of United States primary copper currently is produced, is a large stratiform orebody, 4 to 25 feet thick and several miles across. The present ore column, containing about 1.2 per cent copper as chalcocite and native copper, is confined to the basal beds of the Nonesuch Shale. This upper Keweenawan (Upper Precambrian) formation is a distinctive gray siltstone unit, about 600 feet thick, overlying and overlain by thick red-bed sequences. All these rocks contain abundant volcanic debris, but the latest known igneous activity within 50 miles antedated the Nonesuch Shale. The rocks are unmetamorphosed and only moderately deformed, mainly by tilting and faults. Pyrite, presumably syngenetic, is the only prominent sulfide throughout most of the Nonesuch Shale. The basal beds of the formation, however, contain disseminated copper minerals instead of pyrite and constitute the cupriferous zone. On a regional scale, the top of the cupriferous zone cuts across bedding, and ranges from an inch or so above the base of the formation in some areas to more than 50 feet in others. However, within the cupriferous zone the distribution of copper shows remarkable stratigraphic control. Certain beds are copper-rich compared to the beds that lie between them and the differences persist throughout the White Pine region. Several darkgray thinly laminated beds in the lower 10 to 25 feet of the Nonesuch Shale are rich enough and sufficiently close enough together to constitute a mineable thickness where they lie within the cupriferous zone. Pyrite is the characteristic sulfide of the Nonesuch Shale as a whole. The principal sulfide mineral of the cupriferous zone at the base of the formation is chalcocite. The succession of minerals outward and upward from the center of the cupriferous zone is native copper-chalcocite-bornite-chalcopyrite- pyrite. The bornite and chalcopyrite are confined to a narrow fringe at the margin of the cupriferous zone. This pattern of mineral zones is believed to represent reaction between syngenetic pyrite and introduced copper. The main copper mineralization is independent of structure and apparently antedates deformation. Details of the geometry of the cupriferous zone indicate the possibility that copper was introduced from the underlying sandstone prior to lithification. The source of the mineralizing solutions and their copper is not known.

Jan 1, 1968

By Joseph H. Bilbrey

COBALT is a silvery gray metal which closely resembles nickel in appearance and properties although it is somewhat harder and more brittle. It is rather widely disseminated and is estimated to make up about 0.001 percent of the earth's crust as compared to 0.020 percent for nickel. Cobalt, one of the four ferromagnetic elements, is used for many magnetic applications. It is oxidation resistant, a property which when coupled with its moderately high melting point of 1,495° C., makes it useful in high-temperature applications. Although cobalt had been used by man as a coloring agent for glass and other ceramic products for 4,000 years, it was not isolated as a metal until 1742 and not established as an element until some 40 years later. Early in the 20th century, cobalt began to be used in alloys, and soon these applications displaced ceramics as the major consumer of cobalt. The use of cobalt in alloys gave considerable stimulus to production, and the world output of cobalt increased rapidly, rising from 400 tons in 1920 to 1,400 tons in 1930 and 5,500 tons in 1940. With the advent of the jet plane in the mid 1940's and the consequent increased need for cobalt-bearing high-temperature alloys, cobalt production again increased, reaching 7,900 tons in 1950. The increased demand for cobalt during the Korean war coupled with a large demand for stockpiling caused world production to more than double within the next 10 years, reaching a peak of 17,300 tons in 1959.

Jan 1, 1962

By C. Wen

In Huoqiu region of Anhui Province, Fe2O3 occurs mainly in the form of specularite. Besides quartz, the gangue minerals include iron-bearing silicates such as chlorite and sericite which account for about 10% of the total gangue. The separation effect of conventional collectors is poor in the flotation of the minerals, that is, the average grade of iron concentrates is only 60-62%. Using the developed amphoteric collector can successfully realize not only the separation of specularite from various gangue minerals with results of iron concentrate grade above 65% and iron recovery of above 80% but also low temperature flotation at 15°C.

Jan 1, 2012

By Phillip J. Mackey, Gerardo R. F. Alvear F.

It is now 38 years since the publication of the first edition of "Extractive Metallurgy of Copper" by Professor W.G. Davenport and the late Professor A.K. Biswas. The text book would be issued in a total of five editions, soon became the standard for students and those in industry. Interestingly it was published at the onset of great change in the industry. This was a time when sulphur dioxide controls in the USA were introduced, and almost simultaneously, world oil shocks caused significant increases in the cost of fuel - all leading to efforts to adopt new technology and move away from the reverberatory furnace, which had no potential for sulfur capture and required large amounts of fuel. At the time of publication of the book, new processes such as the Noranda Process and the Mitsubishi Process were in their infancy, while flash smelting (which was first commercialized in 1949 at Harjavalta, Finland) had made significant in-roads in boosting copper output, with some 22 new flash smelters licensed by 1976. Computer technology and associated process control techniques were rudimentary compared to systems today, available physico-chemical data remained largely in tabulated format and needed to be manually entered into the main frame computers of the era, while modeling techniques common today had not yet been developed. The copper industry has therefore seen immense changes since the time of this textbook: world tonnage of smelted copper has risen from 7.34 million tonnes/year in 1976 to 16.7 million tonnes/year in 2012, over a two times increase. Technologies such as the El Teniente Process, or Top Submerged Lance technologies, typified by the ISASMELTTM, not commercialized at the time, now treat large tonnages of copper concentrate in a single furnace. Advanced heat and mass balance software and thermochemical modeling are now used extensively to examine chemical and physical characteristics for process design and to support process optimization. These changes in the industry have been chronicled in subsequent editions of the book. This paper looks at the issues and problems seen in 1976 and over the intervening years, and compares this with the challenges facing the industry with the available tools and technologies today. Some ideas looking into the future are also discussed in the paper.

Jan 1, 2014

By K. Lombard

"‘Roof support awaiting time’ (RSAT) is a term used at Goedehoop Colliery’s Simunye Shaft to describe the potential production time lost due to the continuous miner (CM) standing idle waiting for roof support to catch up. Investigations revealed that in 2013, Simunye Shaft had approximately1400 hours of RSAT, which suggests that the mine could have potentially produced an additional 280 000 t of coal. This project consisted of two parts. Firstly, the causes of the high RSAT and means to improve the situation were investigated. Secondly, as insisted by mine management, the CM cutting sequence was investigated as a possible cause of high RSAT. Machine-related challenges due to the roofbolter installing support too slowly, geological conditions (mostly hard roof conditions and slips),logistical challenges pertaining to the CM cutting sequence, man-related challenges related to operator fatigue, re-support, operator inexperience, and the absence of support targets were identified as main contributors to RSAT. Furthermore, results showed that the roofbolters in the sections at Simunye Shaft are slower than the CMs. A target of 28% reduction in RSAT was set. Experts from Kennametal and Fletcher were consulted to find solutions for the identified causes. In total, eight solutions for RSAT were identified, but the solution that contributed most significantly to reducing RSAT was to use hard roof drill bits as a standard product at Simunye Shaft. Calculations showed that by using hard roof drill bits, RSAT can be reduced by 43%, which is more than the specified 28% target. The cutting sequences at Kriel, Greenside, and Simunye Shaft, togetherwith three newly developed cutting sequences, were simulated using the UCMS (Underground Coal Mining Simulation) program. A re-aligningprinciple was incorporated into the newly developed cutting sequences toalign the cutting sequences to general support work and to reduce RSAT. A decision matrix revealed that a cutting sequence in which boxing takesplace in R3 (third road to the right of the belt road) and in which the re-aligning principle has been incorporated will be the best option for Simunye Shaft. The recommended cutting sequence will lead to a 5% increase in production.The work presented in this paper was carried out as partial fullfilment for the degree BEng (Mining Engineering)"

Jan 1, 2015

By Alfonso Guzmán da Silva

El análisis del trabajo señala la importancia y aplicación del proceso Sink & Float en la aceleración de la industria minera metálica, la no-metálica y la del carbón del Perú. Dicho procedimiento que implica la separación de una mezcla de partículas de densidades distintas a través de la diferencia de las mismas, tiene como ejemplo más próximo de su ejecución, a la compañía minera Milpo S. A., así como la primera planta instalada para la pre-concentración de minerales en el Perú, la corporación minera Castrovirreyna, en Huancavelica.

Jun 20, 1954

By J. A. Theron

"This paper provides guidelines on performing mass and energy balance modelling involving coal and coal derivatives. Usually, the inputs to a pyrometallurgical process would be specified in terms of elements and compounds. Reliable thermochemical data is more widely available for species involving uniquely defined, relatively smaller molecules. However, in the case of coal, the molecules are extremely large and not uniquely defined. Consequently, modelling processes involving coal and its derivatives involve several potential pitfalls. These are outlined in the present paper.It was found that coal proximate analysis should not be regarded as absolute; it could vary with several parameters, including heating rate. For modelling, the use of ultimate analyses should be considered a preferable option to proximate analyses, where ‘fixed carbon’ and ‘volatiles’ are not defined in terms of chemical composition. Significant errors could be incurred if the larger molecules are neglected during calculation of the calorific value (CV) of coal gas (the gas liberated when coal is heated in the absence of oxygen).For elemental analysis determination, the oxygen content (which is calculated by balance) should be checked to ensure it is within the expected range. For representation of sulphur in coal, one should avoid doublecounting due to SO3 in the ash analysis.Potentially, oxygen in coal could be represented as O2, H2O, CO, or CO2.However, use of some of these species without considering the experimentally determined gross CV leads to significant errors in the energy balance. If coal enthalpy is calculated from elemental analyses without correction, representation of coal oxygen as H2O(l) gives reasonable accuracy. Coal volatiles could be represented by a complex mixture of compounds, even using different oxygen-containing species than these four, provided the enthalpy is corrected.It is recommended that an ‘enthalpy correction value’ be incorporated in energy balances involving combustion, devolatilization, or conversion of coal and coal derivatives, e.g. coke, char, or tar. That would imply that proximate analysis, elemental analysis, as well as the gross CV would be required for all solid or liquid coal-derived substances being modelled. No other correction due to carbon being present in a form other than graphite should be used, as that would imply double-counting some effects."

Jan 1, 2015

By B. J. Barfield

Standard erosion plots (4.6 m by 22.1 m on 9 percent slope) were constructed on the Alston Surface Mine of Peabody Coal Company in Western Kentucky to evaluate the effects of placement technique and density on the erodibility and eroded size distribution from mine spoil and reconstructed topsoil. PLots were constructed with a dozer and scraper pan from shale spoil material, topsoil (A horizon) and subsoil (B and C horizon mixture) using soils stockpiled at the Alston Surface Mine. Plots were evaluated under fallow conditions with a tilled plot surface and two compacted conditions. One compacted condition was generated by tracking the surface with a dozer and rolling with a loaded scraper pan. The second was generated by rolling the surface with a solid wheel compactor. Rainfall was supplied with the Kentucky Rainfall Simulator at a rate of 63.5 mm h-l with the plots in three different soil moisture conditions; dry, field capacity, and very wet. Runoff was evaluated for flow rates, sediment concentrations, and aggregate and primary particle size distributions. Sediment yields and the Universal Soil Loss Equation E130 indexes were calculated from the data for each plot. The erodibility was calculated from the tilled plot data using a weighted soil loss and E130 value from the dry, field capacity and wet runs, and a cover management factor calculated for every other condition. For the topsoil and subsoil plots, sediment yield was greatest for the dozer scraper pan plots and lowest for the tilled plots. For the shale plots, erosion rates from the wheeled compacter plots were the highest and tilled plots the lowest. Average erodibilities were 0.14, 0.38, and 0.41 tons ac h ft-1 tonf-1 in-1 hundred ac-1 for the shale, subsoil, and top-soil plots. Cover management factors for each plot were 1.29, 1.42, and 1.17 for shale, topsoil and subsoil plots compacted with the solid wheeled compacter and 1.29, 3.79, and 2.51 for shale, subsoil, and topsoil plots compacted with the dozer scraper pan conbination. The accuracy of the Wischmeier erodibility nomograph was evaluated on the data collected. Erodibilities were accurately predicted for the topsoil plot and poorly predicted for the sub-soil plot. The accuracy of the CREAMS procedures for predicting the composition (size distribution) of eroded particles was also evaluated. The composition was well predicted for the topsoil plots, moderately well predicted for the shale plots and poorly predicted for the subsoil plots.

Jan 1, 1983

By Chas. F. Jackson

This bulletin is one of a series of Bureau of Mines reports dealing with mining methods, practices, and costs. A study of this subject was undertaken by the Bureau in cooperation with a large number of mining companies in 1928. This study continued for 5 years, during which engineers of the Mining Division of the Bureau visited most of the larger mines and many smaller ones in the United States, as well as several mines in Canada and Mexico. As a result of this study, 75 reports dealing with underground mining methods and costs at different mines were published as information circulars. Circulars also were issued on 6 open-pit mines and 2 gold-dredging operations. These circulars were written by the managers, superintendents, and engineers in charge of the operations, in some instances in collaboration with Bureau engineers. They were prepared according to a standard outline drawn up by the Bureau to obtain uniformity in textual treatment and in the arrangement and basis of setting up cost and operating data.

Jan 1, 1936

By S. M. Spottiswoode

Much as one hopes to predict whether damaging seismicity might occur somewhere in the next shift, this is not possible at present. Recently, an analysis of seismicity at two mines has shown that a small increase in seismicity occasionally occurs before both large and small events. This analysis is extended here to five mines, with the same results. Earthquake seismologists forecast hazard in terms of the average incidence of earthquakes over the last 50 years, adjusted on an hourly basis according to the possibility of aftershocks of events that have happened. The rate of seismicity in any area remains approximately constant when expressed as seismicity per area mined, or better still, per unit of elastic strain energy released. It is suggested that seismic hazard in mines be estimated using historical seismicity and forward modelling, adjusted according to the possibility of aftershocks. In summary, medium-term forecasting of seismicity is viable; short-term prediction of large events is not.

Jan 1, 2010

By R Zhou, B J. Monaghan, S J Chew, R J. Longbottom, D J. Pinson, H Li, P Zulli, G Zhang, L Lu

Steelmakers are increasingly motivated to use iron ore resources with a wider range of grades and mineral types that were previously considered unsuitable or uneconomical for sintering. In this respect, specific issues, which now require further investigation and understanding, include the sintering of ore types with:- overall higher gangue content- elements that cause problems in the steel manufacturing operations- distinctive sintering performance compared with traditional hematite-rich iron ores. Laboratory-based investigations of the sintering performance and the behaviour of gangue impurities during sintering is an important step towards the successful utilisation of these resources in steelmaking.In terms of laboratory-based sintering investigations, there are two generally accepted experimental scales utilised prior to undertaking industrial-scale trials, which include: bench-scale tablet (compact) testing and pilot-scale pot testing. The former has the advantages of flexibility and more precise control of sintering conditions; but it is difficult to simulate the heterogeneity of a sinter bed. The latter, which is carried out in pots with diameters between 150–500 mm, simulate the commercial sintering process and conditions quite well, but is time-consuming and labour-intensive due to operational complexity. It is therefore desirable to establish a smaller scale pot testing facility that avoids the disadvantages of the normal sintering pot but can provide comparable information.To explore the feasibility of small-scale pot testing, a ‘millipot’ facility was established to examine the sintering performance of iron ores and other non-traditional ferrous materials. The millipot has a diameter of 53 mm and height of 400 mm. The facility has the capability to monitor the sintering temperature at a number of locations (axial, radial) within the sinter bed and to adjust the suction (‘windbox’) pressure. This paper describes the development of the millipot facility and operational procedures, together with initial tests undertaken in order to match its performance with industrial sintering performance.CITATION:Li, H, Zhou, R, Pinson, D J, Zulli, P, Lu, L, Longbottom, R J, Chew, S J, Monaghan, B J and Zhang, G, 2017. A millipot set-up for sintering investigation of iron ores, in Proceedings Iron Ore 2017, pp 85–90 (The Australasian Institute of Mining and Metallurgy: Melbourne).

Jul 24, 2017

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 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