Search Documents

By James H. Boettcher

INTRODUCTION A complex interaction of worldwide economic and political forces is increasingly requiring 3 primary participants for the successful development of large natural resource projects in developing countries. They are the host government of the developing country (the "host government"), one or a consortium of international natural resource companies, (the "Resource Company") and one, or more likely, a syndicate of international lenders. The host government typically controls the rights to its country's natural resources and establishes the economic ground rules for their development. However, because of worldwide inflation and the depletion of easily accessible, high-grade resources, the nominal and real capital development and operating costs of projects have been spiraling. Moreover, the depletion of high-grade and mineralogically uncomplicated resources also raises the technical and operating risks associated with resource recovery. Thus, despite increasing political self-consciousness and economic nationalism in developing countries, international resource companies are still often necessary project sponsors since they can provide, among other things, the required technical and management expertise and equity capital. The presence of these factors, and often the addition of worldwide market access in a project, is often critical in determining whether or not the project will be considered creditworthy by international commercial lenders. The international commercial lender is often the critical third party in project development, due in a large part to the ever escalating-costs of project developments relative to the financial resources of many companies and even a growing number of countries. In addition, lenders have other reasons for taking more than a passive interest in all aspects of a project's development. These reasons include the following: 1. Projects are often financially structured such that the primary security for a loan, once completion and performance standards are met, is the cashflow of the project rather than recourse to its sponsors. 2. Lenders, in addition to the "project risk" just described, are frequently asked to assume various elements of "political" risk which nay put them in a sensitive economic position between the host government and international resource firms (This is discussed later). 3. International banks today have large loan portfolios of developing country debt with final maturities ranging from one to fifteen years. The various current projects undertaken are the building blocks of most developing country economies that will hopefully contribute to the long run economic development of a country as well as generate the foreign exchange necessary to repay foreign loans in the future. In this regard it is useful for lenders to realize whether a given project "makes sense" in terms of a country's natural resource endowments and competitive advantages over other countries. The type of framework discussed in this article would help lenders to better make such assessments. In order to facilitate the continued international flow of capital and technology to developing countries and of readily available supplies of raw materials from them, it is important for each of these three parties to increase its awareness of the goals and objectives of the other parties as well as the bargaining process that results in the legal, financial, credit and fiscal structures by which a joint venture project ultimately proceeds. Each of the parties faces increasingly complicated accept/reject decisions when they are faced with difficult choices among. a wide range of alternative structural combinations each of which has different implications for risk and reward to them. Therefore, by increasing the mutual, 3-way awareness of the decision processes of the other parties, it should be possible to: • provide better information and results in the formation of more rational expectations by all parties going into negotiations, • facilitate the project formulation process by allowing parties to focus more clearly on the areas of common interest and conflict, • result more frequently in final structures and agreements that better reflect the economic and political realities associated with a project and thus • result in more creative, innovative agreements that should reduce the need for changes to them and thus promote the kind of long run political stability conducive to attracting capital and further investment.

Jan 1, 1982

By H. Lyn Bourne

Daily everyone depends on the great variety of glass products, so much so that glass is often taken for granted. In fact most people do not realize how versatile glass has become. Consider the various uses and then try to imagine a day in which we are not influenced by glass. Common uses include container ware, table ware, window glass, lead crystal, automobile glass, and fiber glass. Several less common, but important, uses include laboratory ware, pharmaceutical, TV bulbs, light bulbs, glass ceramics, optical glass, fiber optics, and laser glass. Corning, Inc., a leader in specialty products, uses nearly 1 000 different compositions to manufacture about 60 000 different products (Edwards and Copley, 1977). Glass is such a complex product that-definitions vary and exceptions can be found for most definitions. Glass is an inorganic amorphous (non-crystalline) solid. Most glasses are produced by melting of a mixture of oxide raw materials, and then cooled to room temperature. Soda-lime-silica composition.s account for about 90% of all glasses melted (Anon, 1973). The properties of the glass product come mainly from its chemical composition. All of the different glasses require melting a combination of raw materials and forming the molten material into the desired shape. Both the melting and the forming processes use sophisticated technology and these technologies require experts to manage these production systems. The manufacturing process is continuous and takes place in tonnage quantities, so adjustments in the batch to achieve the desired finished product requires a great deal of expertise. Raw materials are fed to the batch mixing area in very large quantities (tons in most cases). As a result, impurities in the range of 0.1% result in addition of that impurity within the molten glass in kilogram amounts. More than twenty different industrial minerals are consumed in the manufacture of various kinds of glass (O'Driscoll, 1990). This chapter describes the major and minor ingredients of the various glass batches. It discusses the roles of the various oxides in the glass batch and most importantly considers the mineral raw materials which supply the glass industry. Each of the raw materials is described in detail in other chapters so the geology and mineralogy sections are kept brief here. Container glass, by far, accounts for the most production; followed by flat glass, fiber glass, and specialty glass of which table ware accounts for the greatest tonnage. [Table 1] shows the general production data for 1987 through 1990. Statistics for many of the uses do not appear because production volumes are small compared to the major uses. The glass industry is organized in four categories: containers, flat glass, fiber glass and specialty glass. The US Department of Commerce, Bureau of Census, publishes production data about the glass industry in three different categories: 1) glass containers, 2) consumer, scientific, technical and industrial glassware, and 3) flat glass. The Bureau has very complete statistics about the glass industry in these three categories but they report production data in different units according to industry standards. Therefore, [Table 1] gives the production data in dissimilar units. The production of most glass articles follows similar steps. The raw materials are mixed and the resulting batch is fed into the furnace. In soda-lime-silica glasses melting begins between 600 and 900°C. At these temperatures CO, and other gasses are released which create bubbles in the molten glass. To remove the bubbles and insure complete melting the temperature is raised to between 1 500 and 1 600°C. This is the melting-refining stage during which the refining agents in the glass batch serve to aid in the release of gas bubbles, homogenize the melt, and prevent the formation of scum on the surface of the molten liquid. At the conclusion of the melting-refining stage the glass is too fluid for working and the melt is cooled to about 1 100°C to attain the proper viscosity for working and forming to begin. After the glass article has been made, it must undergo annealing (slowly and uniformly reheated and cooled) to remove thermal stresses that were created during the forming process.

Jan 1, 1994

Introduction The objective of the design and operation of a mine is to make as much money as possible within certain reasonable and responsible constraints. The imposition of constraints on the objective of maximizing profits can result in an optimum mine design. This paper deals with the maximization of profits. Optimization will be discussed in a future paper. In the mining industry, it is generally accepted that the maximum economic recovery is defined by identifying the reserves and the sequence of mining that results in maximum net present value (NPV) or discounted cash flow-rate of return (DCF-ROR), as indicated in Fig. 1 (see Lemieux, 2000). Current parametric analyses techniques, such as varying the mining costs; profit reservation; sales price or metal content, to define pit limits and mining sequence technically referred to single, double or triple parameterization, do not assure that the curve in Fig. 1 will be maximized. In defining the sequence, the first material to be mined should be the first grouping of recoverable mineral that has the highest after-tax unit value and that forms a feasible mining unit. The assumed mining of the first pit should be followed by the assumed mining of the second most profitable and practical pit or pit expansion. If this process of grouping and sequencing is followed until the profitability approaches zero, the general sequence required to extract the maximum economically recoverable reserve is identified. The curve in Fig. I is constructed by calculating the NPV for a series of sequences, each assuming operations are terminated on successively lower profit pits. If the operation ceases on a pit with too high of a profit, opportunity is lost, and the NPV is lower than the maximum. The NPV curve peaks and then starts to decline before the zero-profit pit is assumed mined. This occurs because the money invested in advanced stripping would have produced a greater return if it was invested at the dis¬count rate. The maximum economically recoverable reserve is defined by the pit limit corresponding to the maximum NPV or DCF-ROR. Pit and phase definition The real challenge of reserve definition is identifying the phasing and after-tax profit per ton, so that the curve in Fig. 1 can be plotted. The after-tax profit cannot be properly addressed until after the pits are designed, the production is scheduled and the cash flow is estimated. Throughout the years, many investigators have ad¬dressed this problem. A major thrust of the famous paper in which Lerchs and Grossman (1965) presented their algorithm was to identify a sequence that would produce a result similar to that identified in Fig. 1. Whittle's (1988) four-dimensional analysis is designed to address this issue. The author addressed this same issue in a paper presented to the 1968 Canadian Institute of Mining national convention (Lemieux, 1968). Without a relatively simple methodology to con¬struct the curve in Fig. 1, the number of iterations in a trial-and-error method becomes prohibitively burden-some. The pit designer seeks a simple pit-planning parameter that will serve as a proxy for the after-tax profit per ton used in defining the curve in Fig. 1. The designer's goal is to apply this pit-planning parameter using standard pit-wall location techniques to design a series of pits that will provide a guide to the sequencing. Common techniques used to position pit walls are strip-ratio limits, highwall incremental analysis, floating-cone methods or "maximizer" applications1. Costs, revenues, noncash charges, taxes and profits are usually structured on a unit basis for use in the pit-wall location analysis. Application of the proxy parameter using the pit-design tools should, ideally, result in the definition of pits and sequences that will maximize the curve in Fig. 1. 1 The pit-design tool frequently referred to in the literature as an "optimizer" defines a pit boundary that maximizes the value contained within based on the input parameters and the safe pit-wall angle. The optimization of the pit design involves selection of rate of production, cutoff grade, product quality and other considerations. These considerations constrain the maximization of value. Therefore, pit-design tools that maximize value should be called "maximizers."

Jan 1, 2000

By E. Wainer, K. D. Burnham

A wet beneficiation technique for producing wollas-tonite from its ore in high yield and purity has been evaluated in a pilot plant operation at the rate of 75 tons per month. Finely crushed, unsized wollastonite ore mixed with water is passed through a high in tensity oscillating wet magnetic separator of unique design and in a single pass over 90% of the wolla-stonite exhibiting a crystal purity of 99.8% is obtained irrespective of size of feed fraction. Thereafter, the material is processed by closed cycle standard wet milling and classification techniques to yield a 1µ size product, though any mesh size up to and including 50 mesh may be obtained, if desired. Costs appear to be significantly lower than those available from dry processed techniques. Wollas tonite made by this unusual wet process appears to have potential utility in the following fields: ceramics, paints, plastics, paper, organic finishes, reinforcement of portland cement items, controlled porosity refractory ceramic foams, cinder and concrete block paint, and the like. The 1µ and certain chemically processed varieties of wollastonite may have unusual utility in the paper industry both as a filler and as a coating material and in the organic finish industry. Extensive deposits of wollastonite ore equivalent to an average tenor of 50% to 60% of this latter mineral in easily separable form exists in and around Essex County in northern New York State in reserves of the order of several tens of millions of tons. While important portions of these deposits are susceptible to open pit mining techniques, one operation near Wills-boro, N.Y., involves tunnel mining and dry milling and beneficiation techniques. This mill and tunnel mine is presently being operated by The Cabot Corp. and a variety of particle sizes are now sold into markets in the ceramic, paint, plastic finish, floor tile and in similar fields in substantial tonnages on a repetitive basis. Serious investigation and study over the past several years has indicated that, outside of the obvious economies of open pit techniques for mining purposes, wet beneficiation and milling procedures applied to the ore not only represent a potential for greater economies in the production of a superfine finished product but yields products of improved properties exhibiting increased market potential which may not be available from the dry ground product. The Lewis and Deerhead deposits, controlled by the Adirondack Development Corp., appear to be identical to the Wilisboro deposit. Utilizing ore taken from the Lewis and Deerhead deposits, a pilot plant process for the wet beneficiation, milling and classifying of wollastonite ore has been operated for several months. After scalping a marketable garnet product from a minus 16 mesh dry crushed feed on a high intensity roll magnet the balance of the material is then roll crushed in a closed cycle until it will pass a 50 mesh screen. The product constitutes the feed of the wet magnetic separator. The heart of the new beneficiation process is a cyclically operated wet magnetic separator which exhibits the unique feature that unsized feed is easily handled. Product yields of higher purity are equal to that obtained with dry magnetic separations which use closely sized dry ore and multiple passes, but only produce wollastonite of about 98% mineral purity as determined by sink float techniques. It was anticipated that the wet processing through the grinding, milling and classification stage would yield a low cost 1µ ground product which should make available greatly increased applications for the mineral beyond those presently enjoyed. The improved purity was also expected to provide coarser sizes which might be utilized as a raw material for chemical modification which again would expand the uses of wollastonite. The evidence thus far collected appears to indicate that this premise may be expected to be fulfilled. There are on the surface of wollastonite particle sites at which cheinical reactions may occur1. It is believed that wet ground material provides a better base which will allow wollastonite a deeper entree into the field of chemical raw materials. Wet ground at 20µ and finer permits with certainty a number of chemical reactions, some of which are mentioned later in this article. While well crystallized wollastonite makes up the majority of the ore, the balance consists mainly of very weakly magnetic diopside of the hedenbergite-

Jan 1, 1965

By V. F. Hollister

The Queen Charlotte gold deposit (also known as the Specogna, Babe, or Cinola) was discovered in late 1970 by Efrem Specogna and Johnny Trico. They were prospecting along the trace of the Sandspit fault, and they sampled jarositic material that included veins from the fault zone. Their assays included some good gold values. They located the Babe Claims in 1971 and optioned them to Kennco Explorations Ltd. Kennco conducted stream sediment and soil geochem studies, geological mapping and drilling of 55.2 m (181 ft) in two drill holes. Kennco withdrew from the property, and by 1977, Cominco, Placer Development, Silver Standard, and Quintana Minerals had each successively worked on the deposit. Consolidated Cinola Mines bought the claims in 1977 and completed exploration of the deposit ultimately in a joint venture arrangement with Energy Reserves. Energy Reserves entered in late 1979, and by 1983, a total of 200 surface and 12 underground drill holes had been completed, totaling 28 600 m (93,832 ft). In addition, a 461.9 m (1515 ft) adit was driven in the ore body. Reserves published by Cinola aggregate 34 Mt (37 million st) of 1.7 g (0.060 oz) Au per tonne. All industry geologists working in the area agree that the Queen Charlotte gold is an epithermal gold deposit in porous volcaniclastic and clastic rocks that is genetically related to a Miocene-Pliocene rhyolite plug. The geologic description here follows the oral presentation made by G. G. Richards, J. S. Christie, and M. R. Wolfhard at a 1976 CIM meeting. Additional comments by the late C. S. Ney of Kennco Explorations are included with the Richards, Christie, and Wolfhard data to complete the geologic set- ting. The attached reprint by Champigny and Sinclair (1 980) provides additional data on the deposit. The Queen Charlotte (Babe, Specogna, or Cinola) gold deposit is located on Graham Island in the northern Charlotte Islands, at the fault boundary of the Skidegate plateau and the Charlotte lowlands. The Sandspit fault marks the physiographic boundary, and it is intruded by rhyolite porphyry at the deposit. Because displaced geochem anomalies, drain- age patterns, and topography suggest dextral and east side down movement for the fault, exposed east block rocks are younger than rocks west of the fault. The east block, which included the ore deposit, is a lowland largely masked by unconsolidated Pleistocene and Recent clastic rocks. Rocks west of the rhyolite plug, on the Skidegate plateau, are Cretaceous carbonaceous and calcareous shale unconformably overlain by a thin veneer of rhyolite tuff. The deposit is about 100 m (328 ft) above sea level. GEOLOGIC SETTING At the ore deposit, the Sandspit fault is a complex structure. However, lithologies on either side of the fault-controlled rhyolite porphyry plug are distinct. West of the rhyolite, the upper member of the Cretaceous Haida formation shale is the most common rock type. The member is composed of dark grey to black, poorly consolidated, and thinly bedded carbonaceous and calcareous shale. The sequence is silicified to an argillite or a hornfels near the rhyolite. West of the rhyolite porphyry plug, but not importantly involved in the mineralization, is a thin masking cover of rhyolite tuff. The tuff unconformably succeeds the cretaceous shale, but it is largely eroded in the mine area. East of the fault-controlled rhyolite porphyry plug is a thick section of conglomerate and sandstone of the Miocene- Pliocene Skonun formation. The ore body occurs within the coarse clastics, the rhyolite intrusive as the western boundary. Ore is entirely contained within the clastics. Original permeability of the clastic rocks was a major control for ore deposition and alteration, and for that reason, the Skonun formation is described in greater detail. The Skonun formation clastics unconformably overlie Haida shale. The Skonun is at least 300 m (984 ft) thick in the mine area, where it trends northerly and dips easterly 15º to 25º. The sequence is about 62 percent conglomerate, 3 1 percent coarse arkosic sandstone, and 7 percent interbedded siltstone and sandstone. Contacts between adjacent beds are generally sharp. Stratigraphic correlations between drill

Jan 1, 1985

By J. E. Warren, J. H. Hartsock

Because of the extensive utilization of hydraulic fracturing for the stimulation of low-productivity wells, the two related problems of fracture design and evaluation have become economically significant and, as a consequence, helve motivated this investigation. The producing characteristics of horizontally fractured wells were studied to determine the fracture configuration that should be employed as the basis for the design of the treatment and to develop a method that can be used to establish the degree to which the design objectives have been achieved. The equations which describe the steady-state flow of a single-phase fluid into, and through. a finite-capacity fracture were solved numerically for an idealized reservoir-fracture model. The numerical results were used to obtain an apparent skin effect for each combination of the parameters considered. Based on the computed results, subject to the limitcitions implied by the assumptions that were made, the following general conclusions were drawn. 1. For a radius of drainage at least four times us large as the radius of the fracture, an apparent skin effect that is independent of the radius of drainage can be calculated. 2. The productivity of the hydraulically fractured system, relative to that of the unfructured well, can be determined from the apparent skin effect and can be used to establish design objectives. 3. In the evaluation of a fracture job, it is not Possible to determine both the radius of the fracture and its flow capacity uniquely from the apparent skin effect; an independent determination of one of the quantities is necessary. lNTRODUCTION Although hydraulic fracturing has been employed as a method for stimulating the productivity of literally hundreds of thousands of wells during the past 10 years, it is only in the last few years that improvements in fracture design1-6 and fracturing technique' have combined to increase the probability of obtaining a successful treatment to such an extent that the mechanics of the method may be considered to be standardized. From an economic point of view, however, two related questions must be satisfactorily answered before hydraulic fracturing can be used in the most profitable manner. The two questions are the following. I. For a particular well in a given formation, what are the optimum design specifications for the fracture treatment? 2. Have the design objectives been achieved by the fracture treatment? The significance of these questions has been recoguized, and some attempts to obtain answers have been made. Howard, et al,8 endeavored to determine the optimum treatment, based on maximizing profits, for any given formation; unfortunately, this work was based on a crude method for approximating the productivity of a well. Carter and Tracy9 utilized the same approximation to study the effect of fracturing on the behavior of a well producing by virtue of a solution-gas drive. Electrolytic models were used by van Poollen10 to investigate the variation in productivity due to fracturing; however, only a limited number of results were presented. Later, from the same model results, van Poollen, et al,11 attempted to justify an approximate expression for determining the productivity of a fractured well. It is quite apparent that there is a definite lack of the practical information necessary for specifying the optimum fracture configuration to be considered for design purposes. The only detailed attempt to develop a procedure for evaluating the result of a given fracture treatment appears to be that of van Dam and Horner.12 These authors described a technique for analyzing pressure build-down data, obtained immediately after fracturing, to determine the final fracture volume, the final fracture porosity, the fracture area, the fracture thickness and the in situ fluid loss of the fracturing fluid. While this approach should be useful whenever acceptable pressure measurements are available, it does not yield a value for the flow capacity of the fracture. Since the problems of fracture design and evaluation are inversely related, it should be sufficient to study the effect of the fracture configuration on the performance of a well. The primary objective of this investigation is to evolve a technique for computing the desired solutions. The secondary objective is to analyze these computed results in order to prescribe a method for evaluating fracture treatments. Because this study is exploratory in nature, its scope

By P. G. Winchell, M. Cohen, J. Woodilla

The kinetics of stabilization have been studied with respect to the isothermal component of the martensitic reaction in ivon-nickel alloys. Although the carbun (or nit-vogen) content may be very low in these alloys, it plays an important vole in the stabilization phenomenon. In fact, if these interstitial elements are removed by moist-hydrogen treatment, no evidence of stabilization is found. In the pvesence of 0.007 wt pct C, the activation energy for the stabilization process is comparable to that -for the diffusion of carbon (or nitrogen) in ferrite) rather than in austenitc. This suggests that the interstitial diffiision controlling the stabilization occurs within the martensitic embryos rather than in the matrix of the parent phase. The kinetics furthev indicate that the interstitial atoms tend to diffuse -from the embryos toward the suvroundings, thereby immobilizing the interface. This appears to be the origin of the stabilizing effect -found in these alloys. PHILIBERT' has shown recently that interstitial elements play an important role in the stabilization of martensitic transformations in iron-base alloys. However, determination of the temperature dependence for this stabilization process has produced values (- 11,000 cal per mol) that cannot be readily identified with the activation energy for diffusion of carbon (or nitrogen) in either ferrite or austenite. In the present work, the isothermal mode of the martensitic reaction was employed to derive the activation energy for stabilization, and this was found to agree fairly well with that for diffusion of carbon (or nitrogen) in ferrite. An iron-nickel alloy with 30.8 wt pct Ni and 0.007 wt pct C was used in these experiments, the specimens being in the form of rods 0.074 in. in diam and 2 1/2 in. long. The specimens were sealed in evacuated vycor tubes and austenitized at 1100°C for 1/2 hr prior to quenching in water. A second series of specimens was decarburized (and denitrogenized) prior to the above austenitizing treatment by annealing in moist hydrogen for 182 hr at 1215°C. In order to achieve a "standard" condition in the austenitic state, all the specimens were transformed to martensite by cooling to - 195°C before the final austenitizing. Both the 0.007 pct C and the decarburized series then had an Mi temperature of about - 30°C. The isothermal martensitic transformation was carried out at a reference temperature of -88"C, and was traced by means of electrical resistance measurements. Each resistance value was normalized by dividing by the resistance of the austenitic specimen at O°C, thus giving a resist- ance ratio R. The rate of decrease of the resistance ratio (-dR/dt) provided a measure of the transformation rate. It was found that a plot of R vs log t at the reference temperature resulted in a slightly curved line which could be closely approximated by two straight lines intersecting at a time of 10 min, as shown in Fig. 1. If G = dR/dlogt is the linear slope of either straight segment, the rate of transformation at any given time is: The ratio of G for 1 less than 10 min to that for t greater than 10 min was 1.19 for the 0.007 pct C series and 1.24 for the decarburized series. These ratios made it possible to compute the transformation rate at any time greater than 10 min, knowing the transformation kinetics prior to 10 min. The stabilization treatment was introduced by up-quenching from the reference temperature, usually after 10 min of isothermal transformation: to one of four stabilizing temperatures between 0" and 22"C. After holding for various times, the specimens were down-quenched to the reference temperature for additional isothermal transformation. In this way, it was possible to compare the transformation rate (- dR/dt), immediately after stabilization with the transformation rate (- dR/dt)U in the unstabilized condition. The degree of stabilization was defined as: Ordinarily, Gs and Gu would be measured after the same time of prim isothermal transformation.

Jan 1, 1960

By R. A. Greenkorn

This project originated in a practical problem—we needed five tracers that could be used together to locate flow paths in a pilot flood. While tracers for subsurface liquids have been used since the turn of the century,1-15 none of those reported in the literature appeared to be either consistent or quantitative enough for our purposes. Most were used in field systems without controlled experiments to determine the accuracy and precision of analysis, and many were tracers that could not be used collectively. The ideal tracer, of course, would follow the fiuid of interest exactly, traveling at the same velocity as the fluid front. But the ideal is impractical to attain because adsorption-desorption effects cause the tracer to lag behind the front; these effects, plus diffusion-dispersion effects, cause the tracer front to spread more than the fluid front. Thus, our objective was not to locate a tracer that would be ideal for all circumstances but rather, to find one that would approximately follow the fluid, or one that under controlled conditions could be corrected to calculate the movement of the fluid front. We considered tracers satisfactory—(1) if they were easy to analyze; (2) if their breakthrough-elution curves were not too different from those for the chloride ion, a tracer believed to follow the fluid front closely; and (3) if we could calculate from the curves a material balance, at 1.25-PV (pore volume) injected, within 5 per cent of that calculated from chloride curves. Of a possible 35 materials, we selected 13 tracers that could be quickly and easily identified and whose analysis was claimed to be accurate within 5 per cent. Only one of these, tritiated water, was a radioactive tracer. Radioactive tracers are easy to detect even in small quantities, but they require special handling and special equipment. Also, those that can be used together are limited because special equipment is required to separate emissions from the various tracers. Two of the original 13 tracers were eliminated in static tests to determine how accurately and precisely they could be analyzed, and to check on gross adsorption. The remaining 11 were flowed through a 9-ft linear sandstone model, and breakthrough-elution curves were obtained. Finally, three tracers were field-tested as breakthrough tracers. These tests are described in the following sections. The 13 tracers considered in these experiments were EDTA (ethylene diamine tetra acetic acid), fluo-rescein, picric acid, salicylic acid and ammonium, boron (as borate), bromide, dichromate, iodide, nitrate and thiocyanate ions, plus chloride ion and tritiated water. All but the chloride ion and tritiated water were subjected to static sand tests to eliminate the tracers that could not be analyzed quickly and accurately and to eliminate those that showed excessive adsorption. All but two of the tracers, EDTA and salicylic acid, qualified for flow tests on this basis. The procedures used in the static tests (see Appendix A for analytical details) were as follows. Each tracer (in the amount shown in Table 1) was dissolved in 1,000 ml of water. Then, 800 ml of this solution was added to 500 gm of sand, and the remaining 200 ml was reserved as a control standard. The sand-tracer mixture was agitated once a day over a 10-week period, and during this period we analyzed five duplicate samples from this mixture, along with five duplicate samples from the control solution. From the control-sample analyses, we constructed a control chart on which the sand-tracer analysis results were plotted. The control chart, and its use to determine accuracy and precision of analytical results, as well as the amount of adsorption, is described in Appendix B. The results of the static tests are summarized in Table 1, which lists the concentration c (concentration at time t) divided by c, (initial concentration). If there was no adsorption, c/c, = 1. The acceptable deviation from this value varied from tracer to tracer, depending upon the limits of analytical precision established for the different tracers. The best results were obtained with boron, bromide and dichromate, with all values falling within limits of accuracy and precision. Ammonium, iodide, nitrate and picric acid also were satisfactory, although the ammonium results were erratic. Flourescein and thiocyanate were unsatisfactory ini-

By David Yob

INTRODUCTION This paper is intended to fulfill three major purposes. The first of these purposes is to narrate the improvements to the ventilation of Northeast Churchrock Mine and the subsequent reduction of the radiation levels. The second major objective is to pass on to persons unfamiliar with the ventilation of radon and radon daughter producing mines some of the most important characteristics of airborne radiation and the control thereof. The third objective is to describe the use of a computerized digital mine ventilation simulation. This description is not only of the usage, but also some of the important methods and techniques involved in the usage. This paper is not intended for the persons with extensive experience in these areas. However, to those persons who are just becoming involved with either ventilation of mines with airborne radiation problems or persons interested in computer simulation, this paper should be of some interest. It is the author's experience that most of the information on either computer simulation or airborne radiation control either assumes an extensive knowledge on the part of the reader, or does not address the direct application of the information contained in these articles. It is for that reason that this paper intends to concentrate on the actual application of the topics covered. INITIAL STATUS In the first quarter of 1980, the author and others became involved with the ventilation effort of the Northeast Churchrock Mine. At this time, this team began an investigation of the mine ventilation with the intent of reperforming a mine pressure survey. During the course of this investigation, it was determined that the ventilation system was inadequate. During this first quarter, the responsibility for ventilation of the mine was transferred to the author and his team. It was determined that the system that was in use was that of the single entry or haulage return type. In this type of system, there is very little direct control of the airflow, and it is not an effective type of ventilation for mines that experience problems due to airborne radiation. In a single entry, or haulage return, system there is no separate return system provided. In the system, as we found it, the only source of intake air to the ore level working areas was from the track level. This intake air was forced into the ore level, or stope level, via bulkhead fans installed on raises from the track. This air, after being used in the working areas, was then allowed to find its own way back to the exhausting vent holes. There was essentially no control of this air from the bulkhead fan discharge to the vent hole. The only control used was by the sizing of the fans on the raises. The main portion of this air used the haulage and access drifting for return. Due to the large horsepower used in these bulkhead fans and the large resistance both between the areas to the vent holes and inadequate intake ducting, the stope level pressure was higher than the track level. This pressure differential caused a severe recirculation between the stope level to the track level. This recirculation also caused a severe pre-contamination problem to the supposedly fresh intake air, making his air nearly unusable for working area ventilation. [ ] It was found at the time our ventilation effort began, that the radiation levels were high enough to make compliance with the four working-level month per year standard impossible. For this reason, we had to spend several months driving costly development drifts to implement a completely different type of ventilation system. Before describing the system that was implemented to solve the problems that were discovered during our initial investigation, some discussion of the characteristics of airborne radiation due to radon and radon daughters is needed. RADON CHARACTERISTICS AND RELATED TECHNIQUES One of the most notable characteristics of radiation contamination caused by radon and radon daughter decay is that once the contamination has entered the airstream, the radiation levels as measured in working levels will continue to rise without any further contamination. This radiation rise will eventually stop

Jan 1, 1982

By Jürgen F. Brune, R. Karl Zipf, Michael J. Sapko

Seals are barriers constructed in underground coal mines throughout the United States to isolate abandoned mining panels or groups of panels from the active workings. Historically, mining regulations required seals to withstand a 140-kPa (20-psig) explosion pressure. However, the Mine Improvement and New Emergency Response Act (“MINER Act”) requires the Mine Safety and Health Administration (MSHA) to increase this design standard by the end of 2007. This report provides a sound scientific and engineering justification to recommend a three-tiered explosion pressure design criterion for new seals in coal mines in response to the MINER Act. Much of the information contained in this report also applies to existing seals. Engineers from the National Institute for Occupational Safety and Health (NIOSH) examined seal design criteria and practices used in the United States, Europe, and Australia and then classified seals into their various applications. Next, the engineers considered various kinds of explosive atmospheres that can accumulate within sealed areas and used thermodynamic calculations and simple gas explosion models to estimate worst-case explosion pressures that could impact seals. Three design pressure-time curves were developed for the dynamic structural analysis of new seals under the conditions in which those seals may be used: unmonitored seals where there is a possibility of methane-air detonation or high-pressure nonreactive shock waves and their reflections behind the seal; unmonitored seals with little likelihood of detonation or high-pressure nonreactive shock waves and their reflections; and monitored seals where the amount of potentially explosive methane-air is strictly limited and controlled. Figure 1 is a simple flowchart that illustrates the key decisions in choosing between the monitored or unmonitored seal design approaches and the three design pressure-time curves. For the first condition, an unmonitored seal with an explosion run-up length of more than 50 m (165 ft), the possibility of detonation or high-pressure nonreactive shock waves and their reflections exists. The recommended design pressure-time curve rises to 4.4 MPa (640 psig) and then falls to the 800-kPa (120-psig) constant volume (CV) explosion overpressure. For unmonitored seals with an explosion run-up length of less than 50 m (165 ft), the possibility of detonation or high-pressure nonreactive shock waves and their reflections is less likely. A less severe design pressure-time curve that simply rises to the 800-kPa (120-psig) CV explosion overpressure may be employed. For monitored seals, engineers can use a 345-kPa (50-psig) design pressure-time curve if monitoring can ensure that (1) the maximum length of explosive mix behind a seal does not exceed 5 m (16 ft) and (2) the volume of explosive mix does not exceed 40% of the total sealed volume. Use of this 345-kPa (50-psig) design pressure-time curve requires monitoring and active management of the sealed area atmosphere. These design pressure-time curves apply to new seal design and construction.

Jan 7, 2007

By H. S. Price, J. C. Cavendish, R. S. Varga

This paper presents a new technique for solving some of the partial differential equations that are commonly used in simulating reservoir performance. The results of applying this technique to a simple problem show that one obtains accurate pressure values near wells, as well as accurate pressure pdients, which can be explicitly calculated. The method is completely rigorous in that convergence of the discrete numerical solution to the continuous solution for both pressure and pressure gradient is established. High-order, piecewise-polynomial approximations are used near the wells where pressure gradients are steep, while low-order, piecewise-polynomial approximations are used elsewhere to reduce greatly the calculation time. This combination is shown to give a uniformly good approximation to the solution. These approximations, obtained by using a Galerkin process with suitable Hermite subspaces, are shown to be theoretically and numerically superior to the usual approximations obtained from standard finite-difference techniques. Not only are much greater accuracies obtained, but computer times are also greatly reduced. The application of this technique to multiphase flow problems (e.g., single well coning problems) would have considerable practical interest, but such extensions of this technique with full mathematical rigor have not been made as yet. However, the numerical methods presented here are general, and in principle extend to multidimensional, multiphase flow. Moreover, the preliminary results given in this paper are sufficiently encouraging that we feel the effort in attempting these extensions is justified. INTRODUCTION The problem of obtaining accurate pressure distributions and pressure gradients around wells is of considerable importance in the numerical simulation of reservoir performance. The most common approach to solving this problem is to use finite difference techniques (see McCarty and Barfield9 or Peaceman and Rachfordl0). This approach, however, has many disadvantages, the major one being that many grid points are generally necessary for accurately describing the pressure distribution and the pressure gradient around wells. This need for a fine grid results in large computer times and often in prohibitively high costs. Besides investigating the method of finite differences, some authors, such as Welge and weber17 and Roper, Merchant and Duvall,13 have considered a combination of analytical and numerical techniques with some success. These approaches, however, are all nonrigorous and quite often cannot even be applied. In this paper, we present a numerical formulation of high-order accuracy, based on the Galerkin method, for solving this problem. We treat here only the partial differential equation that describes steady-state, single-phase flow. However, the methods presented are general and in principle extend to multidimensional, multiphase flow. Specifically we treat special cases of the problem in two dimensions described by: where G is a rectangle (with sides parallel to the coordinate axis) with boundary ?G, ?/?n denotes the outward normal, and y and B are non-negative constants such that y + B > 0.

Jan 1, 1970

By Robert D. Pehlke, Robert G. Blossey

The solubility of nitrogen in liquid cobalt and a number of dilute cobalt alloys Was measured in the temperature range 1550" to 1700°C at nitrogen pressures to 1 atm. At 1600°C and 1 atm nitrogen pressure, the solubility in pure cobalt is 0.0047 wt pct [3.8 cu cm (stp) per 100 g Co]. The standard free- energy change for the solution reaction 1/2 N,(g) = gout pct in Co) is iF - 10,000 + 5.24 T cal per g atom N. The alloying elements Cv, Cu, Mo, Nb, Tu, W, and V increase the nitrogen solubility, while Al, Fe, Ni, and Si decrease the nitrogen solubility. THE solubility of nitrogen in cobalt was studied previously by Sieverts, Busch and Dodd,2 and Schenck, Frohberg, and Graf.3 Sieverts investigated the solution below the melting point of cobalt, and Busch and Dodd worked with the liquid metal; however, neither study found any solubility for nitrogen in cobalt. Schenck et al, reported a solubility of 0.005 wt pct N in cobalt at 1600°C. Seagle et al. 4 reported residual nitrogen in cobalt after electron-beam melting, although no reasons were ventured for the presence of this interstitial impurity. This investigation was undertaken to clarify these inconsistencies, and to provide additional data on atomic interactions in dilute metallic solution. EXPERIMENTAL PROCEDURE This study was well-suited to the Sieverts' method, which has not changed significantly since its introduction by Sieverts. 5 The apparatus was of the same design as that described by Humbert and Elliott, 6 with a reaction bulb similar to that reported by Weinstein and Elliott. 7 The buret was not water-jacketed, but the measured gas volumes were corrected for small changes in room temperature. The melt temperature was measured by a disappear ing-filament optical pyrometer sighted through a prism onto the melt surface. Temperature calibration was made based on a measurement of the melting point of pure cobalt taken as 1495°C. The difference between the reciprocal of the observed and the true temperatures was assumed to be a constant, neglecting the changes in emissivity with small additions of alloying elements. The total error in gas-volume measurement caused by accuracy limitations in reading the mercury level in the buret, variation in ambient temperature, and variation in melt temperature was about 0.5 cu cm. This represents approximately 10 pct of the total solubility, and is a larger relative error than in previous studies6-' because of the small nitrogen solubilities encountered in the present investigation. The source and purity of the materials used are listed in Table IA. Recrystallized alumina crucibles were used throughout the investigation. Examination of the crucibles after solidification of the melts showed no reaction or penetration. Alloying of the cobalt melt was accomplished with the dropping technique described by Humbert. 6 The investigation was started using cobalt from African Metals Corp. in shot form, but this material was not of sufficient purity and did not yield consistent solubility results. Sherritt Gordon Mines Ltd. cobalt was then used throughout the investigation. In experiments involving alloying elements which could react with oxygen dissolved in the melt, the cobalt was further refined by annealing at 1400°C in a dry-hydrogen atmosphere. This treatment, described by Woolley and Pehlke, 8 markedly reduced the oxygen content. The oxygen and nitrogen analyses of the cobalt used in this study are presented in Table IB. RESULTS Hot Volume. The hot volume was measured prior to each run and was 50 to 55 cu cm (stp) at 1600°C. It increased linearly with increasing pressure and decreasing melt temperature. Pure Cobalt. Measurements of the solubility of nitrogen in pure cobalt and the temperature coefficient of solubility are summarized in Table 11. At 1600°C and 1 atm nitrogen pressure, the solubility was found

Jan 1, 1967

By J. A. Stachura

Continuous mining machinery provides the coal industry with one way to compete for a larger share of the total energy market. Various types of machines are discussed and some of the problems with continuous miners, encountered by operators, are reviewed. Equipment manufacturers are working with mine personnel to provide solutions for problems that arise. While coal production over the past 25 or 30 years has been on a horizontal plane, coal's share of the total energy market has declined. To participate more effectively in this total energy market, it is necessary to produce coal more efficiently. It is the obligation of all management, employes, and mining departments to gear the deep mining industry to the rapid progress and changing of today's modern industry. This can be accomplished in the near future with the selection of the proper type of continuous miner best suited to each operator's individual situation. In most mining operations there is tremendous incentive to undertake the continuous mining program. It can reduce the size of the mine greatly by permitting a minimum of working places; it makes pillar recovery work more efficient from the standpoint of overall cost, amount of coal recovered, and safety. The work force can be reduced materially permitting closer and more efficient supervision. It simplifies maintenance because equipment can be more readily standardized. The trend of the coal market favors the use of continuous mining machines. Although there appears to be a general feeling that continuous mining is still a relatively new program and will be slow in replacing conventional mechanical equipment, the fact is that tremendous strides have been made since the first machines were installed in 1948. This program is advancing at approximately the same rate that mobile loading machines replaced hand loading. From 1948 to 1955 there were approximately 450 continuous mining machines in service. In October 1959, a survey revealed that there were more than 700 continuous mining machines in service. Many operators have expressed a desire to undertake this program, but they feel that they could not do so at this time because of one or more of the following reasons: 1) the thickness of their coal seams, 2) seam characteristics, 3) soft bottoms, 4) bad roof conditions, 5) size consist, 6) insufficient flexibility in machines, 7) difficult ventilation problems, and 8) high maintenance costs. With the realization on coal about the same today as it was in 1948, or slightly less and since coal is still failing to participate to a greater degree in the total energy market, it is not surprising that the coal industry is desperately exploring more economical methods for deep mining. The manufacturers are aware that the coal industry is willing to invest in continuous miners if the equipment is built for maximum flexibility, will produce higher tons per man, and assure long life between overhaul programs. CONTINUOUS MINING MACHINES Before discussing details regarding the selection of a continuous miner, let us have a preview of some of the continuous mining machines which are available to the coal industry today. Jeffrey Manufacturing Co.: The machine shown is the Jeffrey 76 A.M. Colmol. This is their most widely used miner, and has been particularly successful in central Pennsylvania and in high-wall mining in western Kentucky. One of the outstanding features of this auger-type miner is its portability. The entire mining range can be changed from its lowest point to the maximum height without stopping the mining operation. Jeffrey 76 B.M. Colmol: This machine is similar to the 76 A.M. model; however, it is built bigger and stronger for a mining range of 50 1/2 to 72 in. This is the model that is now available (Fig. 2). Jeffrey has added two arms to the top row, omitted the odd arm in the bottom row, thus permitting a 50 pct larger throat opening. This eliminates one

Jan 1, 1961

By Robert E. Green, Edmund G. Henneke

Expressions are derived for the orientation factors associated with dislocation damping due to dislocation motion in the basal slip plane. The results of the calculation of these factors for zinc crystals aye pvesepzted and the relevance of these results to ultvasonic-attenuation measuvements is discussed. In this paper expressions are derived for the orientation factors associated with damping of sound waves due to dislocation motion in the basal slip plane of hexagonal crystals. For zinc crystals, the orientation factors for dislocation motion in one slip direction only, as well as for simultaneous motion in all three basal slip directions, have been calculated at 1-deg intervals. The results of these calculations are presented in the form of contour diagrams in a sector of a standard polar coordinate graph. The calculations have been carried out for both the traveling- and the standing-wave conditions. The relevance of these results to ultra sonic-attenuation measurements is discussed. Associated with the calculation of the orientation factors, the modes of elastic-wave calculation were also evaluated, according to the procedure previously employed by Hinton and Green.' The results of these calculations are not presented here since similar results have been reported earlier by Mus-grave.2 However, it should be mentioned that, although the qualitative features of Musgrave's results are preserved, the quantitative values are somewhat altered since different values for the elastic constants of zinc were used. ORIENTATION FACTORS It is convenient to refer to three different sets of coordinate axes, Fig. 1, in the following calculations. First, there are the Miller-Bravais axes which are used in describing the Miller-Bravais indices of planes and directions in a hexagonal crystal. These axes are denoted in the conventional manner by a,, a,, a,, and c. The second set of axes used is a rectangular Cartesian coordinate system located as shown in Fig. 1. Note that the x axis coincides with the close-packed direction al [2110], the z axis coincides with the c axis [0001], and the y axis is orthogonal to these two lying in the [0110 ] direction. Finally, a form of Eulerian angles is used to locate the position of the specimen axis. The angle between the specimen axis and the c axis is denoted by 0 and the angle between the projection of the specimen axis on the basal plane and the a, direction is denoted by ø. The majority of observations on plastic deformation of metal crystals indicates that such deformation is caused by the motion of dislocations on a limited number of crystalline planes and in a close-packed direction. In crystals with hexagonal symmetry the most prominent slip plane is the basal plane (0001) and the slip directiops- are one of the close-packed directions, [2110], [ 12101, or [ 1120]. To properly set the stage for the following work, several points must be clarified. First, we want to consider the contribution to the attenuation of elastic waves by dislocations only. While the Granato-Lucke, model for dislocation damping is presently the one most quoted, it is not necessary to refer the present work to it. The only assumption that need be made is that the motion of dislocations absorbs energy from the elastic wave and therefore attenuates it. It is immaterial here how this energy absorption occurs. However, because dislocations have been observed to move most freely on the (0001) plane for hexagonal crystals undergoing plastic deformation, it will be assumed that it is the dislocation motion on this plane which attributes to the observed attenuation. The basic orientation factor, Ok, is essentially the same as that used by Green and Hinton for fee crystals,4 defined by the equation

Jan 1, 1968

By J. B. Dorsey, T. W. Brinkley

Performance of a high-pressure, vaporizing gas-drive project in a volatile oil reservoir is discussed. The project has been in operation for over 8 years. Present recovery from the over-all displaced portion of the reservoir is 40 percent, and an ultimate recovery of 60 percent is indicated. Calculations reflect an over-all displacement volume of 50 percent for the entire reservoir; therefore, ultimate oil recovery will be 30 percent of the entire reservoir oil originally in place. Laboratory vaporization tests conducted prior to project initiation also are presented and discussed, along with an analytical method of determining stock-tank oil production and gar injection requirements related to over-all volumetric displacement efficiency from these data. Introduction The Shoats Creek field is located in Beauregard Parish, La., in the extreme western part of the state (Fig. 1). Discovered in 1956, the field is productive in three Cockfield sand zones, ranging in depth from 8,000 to 9.000 ft. The upper two zones, 1st and 2nd Cockfield, are condensate-bearing gas reservoirs. The lowest zone, the 5th Cockfield, is a volatile-oil reservoir. Laboratory work conducted on crude oil from the 5th Cockfield zone shortly after field discovery demonstrated that, when contacted with dry gas, signidcant volumes of the oil could be vaporized at reservoir temperature and at an appropriately high reservoir pressure. Further study led to the conclusion that early initiation of a crude vaporization project in the 5th Cockfield zone would be superior economically to other methods of reservoir depletion. Consequently, the zone was unitized and gas injection operations started late in 1958. This paper discusses laboratory work performed in conjunction with the project, and project performance to the present. History of Development The 5th Cockfield zone in Shoats Creek is composed of a thin, very shaly, low-permeability sandstone. Oil accumulation is controlled by a structural trap bounded by faults along the north and west and by a water contact on the southern boundary (Fig. 2). The zone covers a productive area of 1,550 acres at an average depth of 9,000 ft. The net sand section is thin, having a maximum thickness of 10 ft and an average thickness of 5.5 ft (Fig. 3). Average reservoir rock properties, as determined from core analysis are presented in Table 1. Shaliness of the zone is illustrated by the typical electric log of the 5th Cockfield section in Fig. 4. Because of the thinness and poor quality of the sand, the 5th Cockfield was developed on a wide spacing. Eleven completions were made in the zone between 1956 and 1961 for an average density of 140 productive acres per completion, with the majority of wells being dually completed in the 2nd and 5th Cockfield. Average initial potential for producers in the 5th Cockfield was 130 BOPD. Oil from the 5th Cockfield has a stock-tank gravity of 50.5" API; with a bubble-point pressure of 3,515 psia, the oil was highly undersaturated at the original reservoir pressure of 4,675 psia. Pertinent properties of the oil as determined in the laboratory are presented in Table 2. Because of the volatile nature of the crude, laboratory studies of its susceptibility to vaporization, as well as a vapor-liquid phase study, were conducted in addition to a standard PVT analysis. The vaporization study was conducted in a visual cell; a stepwise illustration of this work and laboratory measurements are presented in Fig. 5. Field

Jan 1, 1969

By T. E. Scott, R. P. Zerwekh

The charactevistics of the deformatiotz of fcc thorium were examined between 78 and 473 K. Activation energy, actiuation volume, and the temperature dependence of the flow stress were determined and conlpared uith several theories for the rate-controlling mechnnistn of the tempemture-dependent component of the flou~ stress. It was established that the Cottrell-Stokes law does not hold, even approximately, for thorium. Through consideration of the above factors and other parameters, it was concluded that the vandom solute barrier theory proposed by Friedel demonstrates the best agreement with experimental observations of the deformation of thorium below ambient temperature. THORIUM, the first metal in the actinide series of rare-earth elements, is interesting and important from two points of view. In the more practical vein, it is a potential reactor fuel. As such, it is important to know its properties. Of specific interest here are the mechanical properties of thorium. Second, no detailed information concerning the deformation characteristics of the actinide elements has been gathered to date. In view of the interesting physical and chemical properties of these metals, their mechanical behavior might reveal characteristics which could improve our understanding of the mechanisms controlling deformation of the more common metals. For example, it has been shown1 that thorium exhibits a yield point which is greatly enhanced by interstitial solute elements such as carbon. This is not a normal occurrence in fcc metals. The goal of this investigation, therefore, was to characterize the deformation of polycrystalline thorium. MATERIAL AND PROCEDURES Thorium, made by the magnesium intermediate alloy process described by Peterson et al.,2 was rolled and swaged to 0.065-in.-diam wire. The metal analysis is given in Table I. After preparation, the oxygen and nitrogen contents had risen to 308 and 60 ppm, respectively. The 0.065-in.-diam wires were cut into 2-in. lengths and annealed in a vacuum of lo-' Torr at 750°C for 1 hr, producing a grain size of 2735 grains per sq mm. The ends of the 2-in. sections were masked off so that a 1-in. gage length could be electropolished. Specimens were polished at 35 v in a solution of 6 pct perchloric acid in methanol which was maintained at 0°C . The gage diameters after polishing were about 0.0625 in. The tests were run on a screw-driven tensile machine at a constant crosshead speed corresponding to a strain rate of 0.005 min-'. The strain rate was changed rapidly after each 0.01-in. extension between 0.005 min-' and one of three higher strain rates which were selected in order to achieve strain-rate ratios of 5, 10, and 20. Tests were run at several temperatures from 78" to 473°K; most were tested at 423", 300°, 215", and 78°K. EXPERIMENTAL RESULTS Flow Stress vs Strain and Temperature. Typical true stress-true strain curves for thorium are shown in Fig. 1 for several temperatures over the range 473" to 78°K. The curves for 78", 215", 300") and 423°K were obtained from strain-rate change tests using the lower (0.005 min-') or base strain-rate values of stress. Additional tests, using the same base strain rate at different strain-rate ratios, gave curves which virtually superposed on those shown, indicating that strain-rate changes did not materially influence the work-hardening characteristics at the base strain rate. The curves for 273") 373") and 473°K were standard tensile tests at 0.005 min-'. The behavior of the flow stress with temperature is

Jan 1, 1968

By C. S. Matthews, G. L. Stegemeier

In one field in South Texas, approximately 72 per cent of the pressure build-up results show a characteris-i.rtic "hump" (i.e., the pressure builds up and then falls off) which makes interpretation by. standard methods impossible Correlation of size and time of hump with formation permeability. well productivitv index, and method of completion led to the tentativc, conclusion that the humps were caused by segregation of gas and oil in the wellbore after closing-in"'. This conclusion was confirmed by performance of simple laboratory bubble-rise experiments, by theoretical bubble-rise time calculations, and by a detailed calculation of PVT behavior in the well of a particular well Oil which accurate Surface and bottom-hole pressure measure-ments were untile. The hump behavior has since been found to occur in many other fields. The cause, however, is not the mine in all cases. In some of these the hump is trace-ahie to leaks in the tnbing which allow influx from the annulus after closing-in. In other cases the hump is traceable to leaks in the device separating pay horizons in dually completed wells. It is concluded that the recording of both surface and bottom-hole pressures is desirable in wells which show an anomalous build-up behavior. A number of field examples is discussed where use of both sets of measurements enables the cause for anomalous behavior to he found, (and a reasonable interpreation of bottom-hole pressure to be made. INTRODUCTION Theoretically the pressure build-up in an infinite reser-voir should be a linear function of In [(t + at)/nere t is the production time and At a the closed-in time. Some of the variations from this behavior arc well known, such as the curved portion immediately after shut-in which results from after-production and skin effect, and the flattened end portion which results from boundary effects in a limited reservoir. The effect of stratified producing zones and irregular geometrical drainage patterns may also contribute unusual characteristics to build-ups. However, the effect of still another phenomenon—that of movement of fluids within the wellbore—has been neglected in most build-up studies to date. Fluid movements within the wellbore after shut-in occur as a result of after-production, packer failures, leaks in the casing or tubing, or buoyancy of the gas phase when both gas and liquid are present in the well. Each of these movements can influence the pressure build-up, sometimes sufficiently to negate use of the data for computing permeability or stabilized build-up pressure. Examples of each of these phenomena are discussed below. FIELD OBSERVATIONS ON HUMP BUILD-UPS In approximately 75 per cent of the wells in a medium-sized field in South Texas the pressure build-up curves rise to 3 maximum and then decline to what appears to be a stabilized reservoir pressure as shown in Fig. I. Similar behavior has also been observed in wells in other fields in Texas and Louisian:,. Because of the possibility of such behavior being a mechanical failure of the pressure bomb. the bombs were run in tanden into one well which had been known to exhibit this behavior. When both bombs recorded identical curves with the characteristic hump, it was inferred that the unusual build-up is a result of well or forma-

By B. A. Eaton, R. H. Jacoby

This paper presents the depletion-performance correlations developed using data from PVT studies and well-completion tests of 27 rich gas-condensate and volatile-oil reservoir fluids. The PVT behavior for each fluid war determined experimentally in the laboratory, and the depletion performance of each fluid was then calculated. Well-completion test data consisting of reservoir pressure and temperature, initial total gas-oil ratio, and stock-tank oil gravity were used with the depletion-performance data to develop general performance correlations. Three correlations are presented. The first equation shows the empirical relation of gas in place per barrel of hydrocarbon pore space with reservoir pressure and temperature, and initial total gas-oil ratio. The second equation is the correlation of stock-tank oil in place per barrel of hydrocarbon pore space. The third equation relates stock-tank oil production (to 500-psia depletion pressure) with pressure, temperature and oil gravity. All three correlations are also presented in convenient graphical form. Gas production to 500-psia depletion pressure is approximately 92.6 per cent of the gas in place for each fluid, so no correlation was developed for gas production. INTRODUCTION The rate of discovery of shallow, low-pressure and low-temperature oil reservoirs in the United States has been declining in recent years. For this reason, the search for oil has gone to greater depths where high pressures and temperatures are encountered. A high-pressure and high-temperature reservoir is the environment in which many rich gas condensates and highly volatile oils are found. Because of this trend, the performance prediction of such reservoirs has become an important task of the reservoir engineer. The need and usefulness of a performance correlation based on field data that are available from the usual well-completion tests have previously been established. These earlier attempts have shown that the phase behavior of rich gas-condensate systems and volatile-oil systems can be correlated with fluid properties, pressure and temperature. However, very little information on gas-condensate depletion performance5 was available at the time of these earlier attempts. The availability of more recent data has made possible the development of an improved correlation. The data, as listed in Tables I and 2, are for 25 gas-condensate and two volatile-oil systems that were studied. The problem was to correlate the data for in-place fluids and the production data with the initial total gas-oil ratio, initial reservoir pressure and temperature, and stock-tank oil gravity. These correlations were obtained by multiple regression analysis, whereby the data were fitted to an equation by the method of least squared error. BASIS OF THE CORRELATIONS 1. The data used in the development of these correlations were obtained from depletion-performance calculations6 and well-test or production measurements. 2. All of the independent variables used in the correlations are available from the usual well-test completion data. 3. The reservoir volume considered is 1 bbl of hydrocarbon pore space (HCPS). 4. The gas-oil ratio used is the total standard cubic feet of gas (at 14.7 psia, 60F) coming off the entire separator train per barrel of stock-tank oil. This minimizes the effect of separator conditions which might vary from system to system. 5. The correlated pressures are the actual reservoir pressures which might or might not be the saturation pressure of the reservoir fluid.

Jan 1, 1966

By F. J. Fayers, R. L. Perrine

For many years the problem of increasing ultimate recovery of oil from a reservoir has been a subject of interest to the oil industry. At present, a standard secondary recovery technique is to flood the reservoir with water once natural forces of the system are depleted. The permeability characteristics of the rock, together with interfacial forces which exist between the two fluids and the porous material. allow only part of the oil in the rock formation to be recovered by the water-flooding process. Laboratory experiments have shown that when detergents are added to the flooding water, the flow properties of oil in a porous rock arc improved. This improvement results in greater oil recovery. Even though higher recoveries are obtained using detergents, the extra cost associated with the method may make it economically unsatisfactory. To establish economic feasibility requires experimental laboratory measurement of fluid and rock properties together with a calculational procedure to predict field results. Pilot tests in the field may also be required. This paper presents a mathematical procedure for predicting detergent flood behavior. A number of reasonable simplifying assumptions are made to obtain a tractable system of equations. The resulting method may be regarded as a generalization of Buckley-Lev-crett' theory to solve simultaneously for water saturation and detergent concentration profiles. Consequently this type of calculation can be used in much the same manner as Buckley-Leverett theory is used to evaluate a water flood. The mathematical method employed is the so-called method of characteristics for solving hyperbolic partial differential equations. This application of the method is instructive because several other problems in reservoir analysis may be analyzed in a similar manner. FLUID FLOW EQUATIONS The flow behavior of oil and water without detergent through porous reservoir rock is usually expressed by Darcy's law and a continuity equation for each phase. If gravitational forces are neglected, vrg = kr (s)µr v p rg.........(1) vrg = kr (s)µr v p rg.........(1) ? pa va = F ?/?t (sd) where v is the volumetric flux, k represents permeability, S is saturation, and P the pressure. Fluid viscosity is denoted by 11, p is the density, 4 is porosity, and I the time. The subscripts w and o refer to water and oil. The capillary forces in the system lead to a difference in pressure between phases. This may be expressed by the relation: Pm -P0 =Pc(S).........(5) The capillary pressure? P,, is determined experimentally as a function of saturation. DETERGENT BALANCE EQUATION Addition of detergent to the water phase affects the flow properties of the system. Capillary pressure and relative permeability become functions of detergent concentration, c, as well as fluid saturation. That is: P--P (S,c) and k,,,.- - -- k,,,. At low oil saturation values, the presence of detergent reduces capillary pressure and reduces the relative permeability ratio. Detergents are surface-active and adsorb on rock surfaces. The quantity adsorbed, a, will depend on the concentration c of detergent solution and may also depend on the saturation; thus, a = a(S,c). In this analysis we have assumed that the detergents used are water soluble but not oil soluble, and their movement is therefore restricted to the water phase. The dispersion of detergent in water is neglected by comparison with transfer by flow within the system. Detergent movement is determined by the water flux and adsorption: ? cv m = -F ?/?t (cs + a) We will assume that the quantity of detergent in solution is too small to significantly change water density or viscosity. When the permeability and capillary pres-

By Gerald R. Coonan

I'm really here to talk about my pet subject of government relations. I've been in Peabody's Environmental Quality Department, working on regulations, for over five years, and I have seen a number of individuals in the coal industry who were completely stunned when a new regulation was brought to their attention. Their first reaction was, "where did this come from"?

Jan 1, 1980