Search Documents

The Crown Minerals Act provides a new+framework for the management of Crown owned minerals and requires that regard be had for the principles of the Treaty of Waitangi. The Act requires the preparation of Minerals Programmes to achieve the Act's objectives of efficient allocation of rights and the recovery of a fair financial return on its minerals. Efficient allocation requires that a permit is allocated to the party which values it most highly. A fair financial return to the Crown must be considered in the context of the economic rent available from individual mineral projects, with the proviso that a nil return to the Crown as the owner of the resource will not be considered a fair return from the community's perspective.

Jan 1, 1992

By Staia. Mariana H.

Thermochemical equilibria of the heterogeneous system Ti-N-H-CI-Fe-Cr-Mo-Ni-Mn-C were calculated by using a commercial software based on the minimization of the free energy. Equilibrium calculations were done for temperatures varying between 800" and l 100°C and an initial molar relationship between reactants, a=NJTiCI,, of 0.1 to 2. The deposition process was considered to occur at atmospheric pressure. The activities of the alloying elements of the stainless steel were considered equal to their molar fractions. It was shown that for the conditions studied the substrate is no inert with respect to the gaseous phase. The existence of liquid chlorides was predicted for temperatures between 800" and 850°C and a ratio between 0.3 and 0.9. The most important metallic chlorides which appears to be present from the calculations in appreciable quantities in the gaseous phase are iron dichloride and manganese dichloride. In the low-temperatures zone and for low initial nitrogen contents, it result that the appearance of the liquid chlorides seems that partially inhibits the production of gaseous HCI. The highest deposition process efficiency calculated is obtained at a temperature 2 900°C and a N2/TiCI, ratios 2 0.7.

Jan 1, 1995

By M. Ataee-pour, S. E. Jalali

So far, a large number of algorithms have been developed for the optimization of pit limits, most of which follow deterministic rules. In this paper, a new algorithm is introduced, which follows a probabilistic logic using the Markov chain process. The algorithm is implemented on a transition matrix that corresponds to the 2D conventional economic block model of the mining area. The probability of mining a block is proportional to the desirability it may provide. Applying this algorithm, the probability of mining each block is obtained, and finally, the optimum pit limits are defined as the pit, which provides the highest probability of mining. In order to validate the proposed model, a numerical example is given and the results of the algorithm are compared to those of a dynamic programming algorithm. A 2D analysis of the problem is discussed in this paper. Although dealing with 3D problems requires larger size matrices, the optimization method is the same and the algorithm does not need any smoothing.

Jan 1, 2006

By Bureau of Mines Staff

Coal research and related activities in 1964 by the Bureau of Mines resulted in many advances in the production, preparation, and utilization of this valuable natural resource. Mining research during the year was marked by work on hydraulic mining, coalbed degasification, mine drainage control, mine sealing, strip-mine rehabilitation, and quenching of mine-refuse fires. Bureau research continued to produce information leading to increased safety in underground mining. Such research included roof and rib control via bolting and bonding, pressure monitoring of coal and mine rock during mining operations, and development of safety devices for detecting conditions leading to roof falls. Other safety work was continued on methane detectors, coal-dust suppression, and mine venti lation equipment and techniques. Mine explosion and explosives research, mine equipment and materials, safety rule modernization, and safety education and training of mine workmen were continued. Coal storage and preparation work in 1964 ranged from freezeproofing of lignite for rail shipment, through washability studies of coals and washing equipment, to research on new ways to remove ash and sulfur from coal. Electric power production, which ,accounts for a major portion of the market for coal, received further study. Unique power generation systems, inc luding coal-burning gas turbine and megnetohydrodynamic generators, were intensively investigated. Conventional power systems were improved through continuation of work on furnace performance, powerplant stack gas purification, and engineering and boiler water service.

Jun 1, 1965

By Michael Hood

Subsidence profiles above two adjacent panels in Illinois are compared with pro- files predicting subsidence behavior obtained using the (i) National Coal Board (NCB, 1975) method, (ii) the profile function method and (iii) the influence function method. The NCB method predicts the maximum subsidence values at the center of the troughs accurately, but, the overall shapes of these predicted profiles do not match the profiles from the measured data well. Consistent values for the angles of draw were measured but these angles were different in the transverse (average value 43 degrees) and in the longitudinal (average value 17.5 degrees) directions. Time dependent subsidence effects are shown to be small but measurable. These displacements continue at a linear rate for at least a 12 month period. A comparison of measured horizontal distances, interpreted as horizontal strain, and the NCB predictions for strain shows that the peak measured strains are greater by a factor f about four than the predicted strains. The relationship between surface curvature and strain is investigated. Problems associated with calculation of surface curvatures from vertical displacement data are highlighted and a recommendation is made for future studies to consider direct measurement of this parameter. Surface curvatures above a moving face are found to be about three times less than the curvatures at the stationary end of the panel. On the other hand, a hyperbolic tangent profile function is shown to serve as an accurate predictive tool for subsidence behavior in two adjacent longwall panels at Old Ben Number 24 mine in Illinois. This function predicts not only the vertical displacements but also the surface curvatures above both panels. Influence functions are shown to be more problematic, although potentially more flexible, in their application than the profile functions.

Jan 1, 1982

By Larry L. Snyder

A compact 4.87 m (16 ft) diameter tunnel boring machine (TBM) is capable of boring a 12.2 m (40 ft) radius curve while gripping horizontally, or vertically, as may be required for the excavation of turnouts. The machine's 5.0 m (16.5 ft) length and a patented machine structure allows steering while boring without causing harmful eccentric loads. An extremely stiff thrust system, rigid front and rear grippers, and the torsional stiffness of a front end electric cutterhead drive combine to create superior machine stability and boring efficiency. A 2.2 m (7.2 ft) diameter prototype, complete with a soft ground shield attachment, was used successfully on a tunnel project in Scotland.

Jan 1, 1989

By Pierre Marcotte

Jackleg rock drills expose their operators to high hand-arm vibration levels. To address this problem, an antivibration handle has been developed. The antivibration properties of the handle were achieved through a suspension based on a hinge design, achieved with a torsion spring and o-rings. To provide some bench marking for this new handle and to follow the evolution of its performance over time, a test bench was developed to characterize the vibration emission values of jackleg drills under controlled operating conditions. Testing of four antivibration handles on the test bench have shown reductions of the overall weighted acceleration level in the order of 54 % to 60 %, with most of the attenuation occurring in the zh-axis at the rock drill percussion frequency.

Jan 1, 2008

By Steven J. Van Kauwenbergh, Talaat A. B. Lawendy

Batch laboratory-scale beneficiation tests were conducted on four high-iron phosphate samples. These samples included sedimentary ores (Tilemsi Valley, Mali and Tahoua, Niger), a sample from an igneous phosphate deposit (Matongo, Burundi), and a fine-grained phosphatic iron ore processing residue from Chile. The Matongo phosphate deposit occurs as a thick (up to 80 m) residual blanket of P205-enriched soil over solid parent rock (carbonatite and associated rock types). Although the Matongo deposit is an igneous deposit, ore from the deposit can be subjected to beneficiation treatments commonly used to treat sedimentary rocks. The beneficiation flowsheets designed to process these diverse ore types involve attrition scrubbing, size classification, and flotation. At the best flotation conditions tested, phosphate concentrates assaying 32.0%-39.0% P205 and 1.0%-1.5% Fe2O3, were obtained at 64.2%-82.2% overall P2O5 recoveries from run-of-mine (ROM) ores containing 9.2%-28.6% P205 and 5.9%-16.1% F,20,. A cleaner phosphate concentrate assaying 33.1% P205, 1.8% A1203, and 3.3% Fe2O3 at 53.4% overall P205 recovery was obtained from the phosphatic iron ore residue that originally contained 4.5% P205, 8.0% A1203, and 17.6% Fe2O3.

Jan 1, 1993

By William J. Campbell

The Federal Bureau of Mines provided five test materials--amosite, short-and long-fiber chrysotile, crocidolite, and nonfibrous tremo1ite--to the National Institute of Environmental Health Sciences (NIEHS) for studies on the biological effects of orally ingested asbestos-related materials. An inter-agency agreement between the Bureau and NIEHS made possible this program to characterize these test materials. The Bureau in turn received support through contracts with the Ontario Research Foundation for the milling of amosite and crocidolite, and with the University of Maryland for microscopic measurements. The five test materials, which were supplied in quantities of 960 to 1,200 pounds each, were characterized by chemical-instrumental analysis for elemental composition, surface area, density, X-ray diffraction, infrared spectroscopy, thermogravimetric analysis, mineralogy, and particle size analyses. Particle size distributions were emphasized because of the current theories relating adverse biological effects to length-width characteristics of durable elongated particles.

Jan 1, 1980

By R. F. Abernethy

Recent investigations in the United States and in other countries have produced considerable information on the nature and occurrence of minor elements in coal. The possibility of recovering certain elements, particularly germanium and uranium, from coal ash and other coal products has received considerable attention. Some flue dusts are potential sources of germanium, and certain lignite deposits are potential sources of uranium, but industrial recovery of these elements has been limited to small quantities for technical and economic reasons. Considering the large amount of coal mined annually, a small percentage of a rare element may have industrial significance if an economical method of recovery can be developed. In 1944 information on the occurrence of rare elements in coal was reviewed in Bureau of Mines Technical Paper 669 (46).3 This report is a revision of that paper with the view of summarizing published information available for United States coals, and pertinent information concerning coals of other countries. Selected references from the original paper are retained. New references do not include all preliminary reports when more permanent literature describing complete projects is available. Consequently, the list of references should not be considered complete.

Jan 1, 1963

Gravity concentration of tin at Renison Limited uses tables and spirals to treat the various coarser fractions produced from sulphide flotation tailing. An acceptable final tailing cannot be made from one pass over the tables and spirals and therefore the tailing from these separators is reground and combined with the sulphide flotation tailing as the feed to the classification section. The gravity section became increasingly overloaded in the years preceding 1976 due to efforts to maximize tin recovery to the gravity concentrate. In April, 1976, two 2DSV Reichert cone concentrators in series were installed to treat the table and spiral tail- ing stream, prior to regrind, to scavenge a small tonnage of hard, coarse low tin grade material for rejection to final tailing. Gravity section loadings were decreased markedly, allowing operator attention to be directed more to table performance than to spillage collection. An increase in mill throughput of approximately 15% was possible following the commissioning of the Reichert cones. The cones are currently operated to make the minimum tonnage of tailing consistent with reasonable- gravity section loadings.

Jan 1, 1978

By Charles B. Kenahan

The chemical and galvanic corrosion properties of commercially pure titanium, a titanium-manganese alloy, a titanium-aluminum-tin alloy, 2 titanium-copper alloys, and 2 titanium-aluminum alloys were studied. Tests were conducted in synthetic ocean water, tapwater, 1-percent sodium hydroxide, 5-percent ferric chloride, and 5-percent ammonium sulfate solutions, various concentrations of sulfuric and hydro-chloric acids, and 1-percent oxalic acid. All tests were carried out at 35° or 60° C., with a flow of air or helium. Dissimilar metals involved in galvanic-couples experiments were aluminum, magnesium, and stainless steel. Electrode-potential measurements on single (uncoupled) specimens were also measured. Commercially pure titanium and the titanium-base alloys were generally quite similar in their resistance to corrosion, whether alone or in. contact with a dissimilar metal. Contact with the titanium metals caused aluminum and magnesium to corrode by galvanic action. Coupling with stainless steel in sulfuric acid solutions decreased the normal corrosion rate of the titanium metals. Contact with aluminum in air-aerated oxalic acid and 0.1N sulfuric acid solutions resulted in cathodic corrosion of titanium and the titanium-base alloys; however, there were a few exceptions. The copper alloys showed better corrosion characteristics than titanium in certain sulfuric acid solutions but were more susceptible to chemical attack in some hydrochloric acid solutions.

Jan 1, 1958

By G. L. Zhu, Y. F. Wu, H. Nagaumi

A novel eutectic-type alloy Al12Si-HTC with high thermal conductivity has been developed by strict chemical composition control and heat treatment. The thermal conductivity of the alloy is about 30 % higher than that of current commercial Al12Si alloy. The Al12Si-HTC alloy has a silicon content of 10~12 wt.%, which provides excellent castability. Trace element strontium or calcium is added to modify eutectic silicon, and contents of transition elements such as copper, manganese are controlled in a low level. By appropriate heat treatment, thermal conductivity of the castings increases by 7.5 % and finally reaches 172 W/(m·K). It is found that eutectic silicon modification which makes average particle size finer than 2µm plays a dominant role in improving thermal conductivity, as well as further spheroidization of eutectic silicon and reduction of lattice defects in a-Al microstructure in subsequent heat treatment process. The alloy is successfully applied to produce a certain communication enclosure with complex finned structure by die casting process in industrial scale. The results of numerical simulation and thermal analysis experiments show that when manufactured by Al12Si-HTC alloy, the surface temperature of chips in the enclosure is reduced by 2.0~4.5 ? and the temperature distribution is more homogeneous compared with traditional Al12Si alloy.

Jan 1, 2015

By Henry P. Cole

A series of studies are described that identify critical health and safety skills needed by underground coal miners for coping with mine emergencies. First aid and self-rescue and escape were selected as two broad performance domains that include many subordinate skills needed for preventing, recognizing, and coping with mine emergencies. The performance domains and their objectives were used to direct the development of simulation exercises that permit miners to experience mine emergencies in the safety of a training room. The exercises are based upon research with similar simulations for the training of flight crews and medical personnel in information gathering, judgment, and decision making skills. Exercise content is drawn from actual mine accident and disaster case studies. The exercises incorporate recent principles of cognitive psychology and instructional design. Initial studies in annual refresher training classes suggest the exercises are stimulating for miners and their instructors, easy to use and read, sufficiently complex and authentic to challenge experienced miners, and useful devices for both teaching and testing critical problem solving skills required for preventing and minimizing emergencies.

Jan 1, 1988

By J. D. de V. Louw

Introduction Zircon is a natural occurring mineral and may exhibit natural low level radioactivity. It contains trace levels (90 to120 pCi/g and 110 to 115 pCi/g respectively or 3.3 to 4.4Bq/g and 4.1 to 4.3 Bq/g) of naturally occurring radioactive uranium and thorium. This mineral is often found in a suiteof minerals (rutile, ilmenite and monazite) Due to the overlap in physical properties of monazite (containing 4–7% thorium) with some of the zircon (slight magnetic susceptibility) it often contaminates the zircon, adding to the levels of thorium in the product. Hazards associated with the products from zircon and its processing is inhalation of respirable dust (generated by milling) containing radioactive uranium and thorium which may cause cancer of the lungs. Legislature regulates the transport, use, workplace environment exposure and storage of radioactive material. In America the general regulation for the unlicensed use of low levels of radioactive materials limits the concentration to 500 ppm (U plus Th). Europe (IAEA) specifies a limit of10 Bg/g (Bequerels per gram) on the parents for transport only which is a higher level of U+Th. Companies selling into these markets need to comply with these regulations and certification is therefore required, not only of the major constituents determining the price, but of the U and Th levels as well. Exceeding the limits for U and Th may result in material being returned to the country of origin. However, for the analyst the requirement is to analyse for these elements throughout the production process as well as the analysis of the ‘tails’ that is returned to the mine area. This requires a wide range of concentration to be analysed. Methods of analysis A wide variety of methods exist for the analysis of U and Th in a range of matrices. (Table I) All these methods have been used with varying degrees of success in terms of accuracy. Overall precision for the different techniques varies from 5 to 25% relative. Time to produce results from these methods vary from a few minutes to several days. The Production Chemist is expected to produce between5 and 10 analyses per hour for process control purposes. A suitable technique that fits this requirement is X-Ray Fluorescence. Both Wavelength Dispersive (WDXRF) and Energy Dispersive X-Ray (EDX) instruments are suitable. X-Ray analysis in general has a few major advantages over competitive techniques. •It is fairly simple to set up and interferences are known and easy to correct •It has multi-element analysis capabilities •Clean spectral lines for U and Th (in zircon specifically) •Uncomplicated and easy to measure background •Low detection limits (3 to 5 ppm) •A very high degree of precision •Sample preparation is uncomplicated and fast (pressed powders and fused beads). In Zircon (ZrSiO4) processed sands the matrix is constant(>99% ZrSiO4) and the background is thus essentially constant. A typical set up for the analysis of U and Th by WDXRFis given below: X-Ray Tube Tungsten (W) Target, 2.5 to 3kW operated at 50mA/50kV Crystal LiF200 (or LIF220 for higher resolution) Analytical Lines (2?)U La1,226.14Th La1,227.47 Detector Scintillation Counter Analysis time, depending on the concentration can vary from as little as 20 seconds up to 5 minutes. Standards As with any other comparative technique XRF analysis relies heavily on the quality of the standards used. Assuming proper set-up and operation of the instrumentation as well as consistent sample preparation, accuracy is determined by the quality of the initial calibration standards. A search for the available zircon certified reference materials (CRMs) shows that there are only a few CRMs

Jan 1, 2003

This report summarizes the results of a research project entitled "Fatal Accident Data Base Development." This project is the last project under Contract No. S0110601. Volumes I and 1I, "Industrial Engineering Study of Hazards Associated with Underground Coat Mine Production," summarize the original industrial engineering project and the Fatal Accident Report Analysis project. These projects were completed and accepted by the U.S. Bureau of Mines on December 10, 1971. Both volumes of the report have been made available by the Bureau to the public through the U. S. Department of Commerce technical information service. Another project under the same con tract, entitled "Equipment Operations Safety Study, ' was completed in January 1972 and the report was accepted by the U. S. B. M. at that time. The analysis in this Natal Accident. Data Lose Development report is closely related to analysis in both of the earlier reports. The original project involved an industrial engineering and safety of production operations for a representative sample of underground bituminous coal mines. The contactor was required to make recommendations 'in work practices, mining methods, equipment, and safety practices which can effectively reduce the fatality rate from rib and roof falls.' (Longwall mining methods were excluded from the study.) Time study and work element analyses were performed in over 50 mines as part of the field data gathering effort. Another requirement of the study was to develop sets of multiflow. Process charts, which are diagrammatic illustrations of the time-related interaction of mining, cycle job elements.

Jan 1, 1972

By Z. Y. Ma, M. D. Liu, Z. Z. Liu, J. K. Zhou

"An optimization method of microwave assisted leaching of copper anode slime was proposed based on Taguchi method. The results showed that the solid-liquid ratio was the most important factor affecting SNR(signal noise ratio) of the copper, tellurium and selenium leaching efficiencies, and its contributions to the copper, tellurium and selenium leaching efficiencies were 58.34%, 55.54% and 60.74%, respectively. Solid-liquid ratio was a significant factor to SNR(copper, tellurium and selenium). Reaction period was a less significant factor to SNR(copper and tellurium) and concentration of sulphuric acid was a less significant factor to SNR(selenium) while microwave power was an non-significant factor, its contributions to the copper, tellurium and selenium efficiencies were 3.22%, 11.73% and 12.73%, respectively. The optimized parameters of the microwave assisted leaching of copper were as follow: microwave power of 500 W, reaction period of 10 min, solid-liquid ratio of 0.10 g·mL-1 and H2SO4 concentration of 1.0 mol·L-1. The optimized parameters of the microwave assisted leaching of tellurium and selenium were as follow: microwave power of 800 W, reaction period of 15 min, solid-liquid ratio of 0.10 g·mL-1 and H2SO4 concentration of 1.0 mol·L-1. The validation experiments under the optimized conditions has been carried out, the copper, tellurium and selenium leaching efficiencies were 98.66%, 84.32% and 64.54%, respectively. INTRODUCTION As important valuable by-products, copper anode slimes are collected from the bottom of the electrolytic bath in the electrolytic copper refining process(Kilic, 2013). The value of the anode slime is determined by the content of rare metals and precious metals(Liu, 2014). Generally, pretreatment of copper anode slime is the first stage of the whole process. The aim of the pretreatment process is to remove Cu, Te and Se etc and enrich the valuable metals. The methods of pretreatment process can be differentiated into: pyrometallurgical processing and hydrometallurgical processing(Amer, 2003). The hydrometallurgical processes are more promising compared to pyrometallurgical processes in terms of cost effectiveness, low energy consumption and environmental concerns(Ranjbar, 2014). However, the disadvantages of the hydrometallurgy process are evident, such as lengthy processing times and low efficiency. Therefore, it is necessary to develop a highly efficient and environmentally friendly pretreatment method."

Jan 1, 2018

By J. Gordon Parr, L. P. Srivastava

DUWEZ1 has shown that pure titanium and pure zirconium transform martensitically during rapid cooling at temperatures about 30° and 15°C re spectively below their To temperatures. Holden et al.2 determined M, temperatures of Ti-Cu alloys by a metallographic technique applied to isothermally treated samples. Their results, superposed on the McQuillan3 version of the partial Ti-Cu diagram are shown in Fig. 1. Sources of metals were: Titanium: Iodide bar; Foote Mineral Co. Ltd. Zirconium: Iodide bar; Foote Mineral Co. LM. Copper: Spectrographic Standard; Johnson Matthey and Co. Ltd. Alloys were prepared by levitation melting under argon to the following nominal compositions:* 0.5, Quenching, at rates to about 8000°C per sec, was performed by cooling in a stream of helium (99.999 pct pure) a 3/64 in. cube of material previously heated by a tungsten coil in a system operating at. 10-5 mm Hg. Cooling rates below 200°C per sec w were obtained by controlling the current to the heater coil. The specimen was supported on a 32 gage Pt, Pt-Rh thermocouple attached through a DC amplifier to an Edin magnetic oscillograph (Model 8082) with a maximum chart speed of 250 mm per sec. Transformation temperatures were read off the oscillograph charts at the location of the beginning of the thermal arrest. The same specimen was quenched several times, and all runs were repeated on new samples. Regular calibrations of the oscillograph against a potentiometer were made between quenching experiments. The criterion of martensite formation was the effect of surface rumpling. Previously polished samples that showed a rumpled surface after transformation are assumed to have undergone martens itic transformation: transformations of a diffusion type produce no surface disturbance, Figs. 2 and 3. Electrochemical measurements were made in a cell described elsewhere4 whose electrodes were the martensitic and the diffusion-transformed materials in the following electrolytes:

Jan 1, 1962

By Jancar T, de Guzman M M, Matthews B, Tu J. Y

The use of computational simulation to predict the performance of existing equipment or for the design of improved mineral processing equipment will be referred to as Computational Mineral Processing (CMP) in this paper. The emphasis will be on solving the governing equations for liquid solid flows in complex geometric domains to obtain, typically, the three-dimensional distributions of velocity and pressure for the liquid phase and the velocity and concentration for the solid phase. The underlying technology, computational fluid dynamics (CFD), is very well-developed in other industries (aircraft, car and chemical) and forms a vital part of product design and analysis. The subsequent design of improved mineral processing equipment may also involve finite element structural analysis (FESA). Computational fluid dynamics (CFD) and finite element structural analysis (FESA) are part of the generic discipline of Computational Engineering and Science (CES). Ciment and Scherlis (1993) define CES as `The systematic application of computing systems and computational solution techniques to mathematical models formulated to describe and simulate phenomena of scientific and engineering interest'. Increasing computer power (Kaufmann and Smarr, 1993), is a dominant factor determining the rapid growth of industrial utilisation of CES. This increase is happening at both the supercomputer and the workstation level. Indeed, today's typical workstation, a Hewlett-Packard 735, is as powerful as a CDC 7600, the leading supercomputer in 1970. Future computer architectures are expected to become increasingly parallel (Gartner Group, 1990) allowing a sustained performance of a teraflop (a million megaflops) to be achieved well before the Sydney Olympics in the year 2000. It is recognised that the modem development of most disciplines in the physical sciences rests on the three complementary strategies of experimentation, analysis and computational simulation. As the 21st century is approached the volume of computational simulation, both fundamental and applied, is growing dramatically and the role of experimentation is diminishing. In the area of equipment design, CES provides the following advantages over experimental testing and measurement: 1. lead time in design and development is significantly reduced; 2. CES can simulate conditions not reproducible in experimental tests; 3. CES provides detailed and comprehensive information; and 4. CES is more cost-effective and time-efficient. The ultimate benefit for industry is greater productivity and hence greater profitability. The broader issues in relation to CFD are discussed by Fletcher (1993b). In the five- to ten-year time frame, growth in computer power will make it practical to combine CFD and FESA with optimisation procedures to produce a design process that is almost fully computerised. It is expected that fully automatic

Jan 1, 1995

By K. E. Starling

Phase equilibria data were used to develop an equation-of-state correlation for complex hydrocarbon mixtures, thereby circumventing difficulties associated with use of pressure-volume-temperature data. Equilibrium phase data for two condensate reservoir fluids were used to determine equation-of-state parameters for hydrocarbons as heavy as C22H46 in a modified form of the Benedict-Webb-Rubin equation of state. Comparative tests of K-values from the resultant correlation were made with data for condensate reservoir, separator and gas plant absorber mixtures. Generally, for temperatures above OF and computed liquid densities below 0.55 Ib-mole/cu ft, the modified BIVR equation predicted K-values in close agreement with the experimental data. INTRODUCTION Accurate predictions of thermodynamic behavior for complex hydrocarbon mixtures are necessary for many calculations in the petroleum industry. Because of the relatively high cost of extensive experimental data, many correlations for prediction of phase behavior have been developed. Some of these, such as the correlation of K-values presented in the NGAA Equilibrium Ratio Data Book, are totally empirical.l Others, such as the Benedict-Webb-Rubin (BWR) equation-of-state method, are semitheoretical.2-4 Unfortunately, published correlation methods often do not accurately predict the phase behavior of complex hydrocarbon mixtures, principally because of inadequate representation of the effects of components heavier than decane. This paper presents a new approach to this problem in which basic equation-of-state relations including heavy hydrocarbon effects are applied. Equilibrium phase data for two condensate reservoir fluids containing hydrocarbons as heavy as C22H46 are used to correlate component K-values. The BWR equation was chosen as the prototype equation of state for this study because of its proven capability for accurately predicting phase behavior and thermodynamic properties of light - hydrocarbon mixtures? Research was directed toward development of BWR parameters for the heavy hydrocarbons and modifications of the mathematical form of the BWR equation for application to complex hydrocarbon mixtures. The new equation-of-state approach presented differs considerably from previous methods in that phase equilibria rather than PVT data are used for determination of equation-of-state parameters. It is an explicit approach since the parameters are determined directly from mixture data. As such, it does not encounter problems inherent in the implicit method used by Benedict, Webb and Rubin 2-4 and numerous other investigators—one in which mixture parameters were postulated to be functions of the pure component parameters. The pure component BWR parameters, in turn, were determined from experimental PVT data. This method has been limited to mixtures containing components lighter than decane because of lack of vapor phase PVT data. Studies have been reported in which BWR parameters have been determined explicitly from PVT data for binary mixtures.6-8 However, since small concentrations of heavier components have only a minor effect on PVT behavior, it is doubtful that these explicit methods would yield useful results for condensate systems. Ellington and Eakin9 have shown that the accuracy of K-values predicted by an equation of state developed from mixture PVT data probably would be more than an order of magnitude lower than the accuracy of the PVT data. On the other hand, an equation of state utilizing parameters developed from phase equilibria data should predict K-values with accuracy comparable to the accuracy of the experimental phase compositions. This work applies this explicit approach with the objective of improving hydrocarbon mixture phase behavior predictions. Results are evaluated by comparing predicted and experimental K-values at reservoir, separator and gas plant

Jan 1, 1967