Search Documents

By Piotr Litwa

The paper focuses on actions aiming at improvement of work safety in mines, particularly in those extracting hard coal, taking into account the characteristics of natural hazards in Polish underground mining and the reasons of mining disasters connected with those hazards. It refers to conclusions addressed to scientific and research centres, drawn during last decade by Commissions appointed by the President of the State Mining Authority in order to investigate reasons and circumstances of incidents, which till now, could not have been accomplished because of financial reasons. Furthermore, the paper presents basic assumptions of the Strategic Research Project Improvement of work safety in mines, announced by the National Centre for Research and Development at State Mining Authority?s President request.

Aug 1, 2013

By Raja V. Rarnani, Robert L. Frantz

"The more than 75 research projects funded throughout the first 10 years of the Dust Center's operation have resulted in numerous cumulative findings fulfilling the goals of the Center. These project descriptions present the methodologies and project-to-project relationships leading to their various accomplishments and demonstrate the mesh between project objectives and Center goals.To develop equipment for collecting respirable dust samples.To develop experimental methods for particle-by-particle characterization of morphology, size, and chemistry.To develop experimental methods for overall respirable dust sample characterization of size distribution, chemistry, and crystallography.The project directly addresses the Center's goals to enhance the understanding of the scientific and engineering characteristics of fine and respirable mineral particles found in mining and milling atmospheres and to establish the relationship between natural and cultural factors (mine geology, mineralogy, mining and milling methods, equipment, and control measures) in the generation of fine particles. The project also aims to contribute to the development of control mechanisms for reducing ambient respirable dust concentrations and worker exposures.The work on this project is related to similar work on Coal Mine Dust Characterization (Project No. 2701/MN1) and Characterization of Dust Particles (Project No. 4203/PS3).Airborne dust sampling was performed at the following locations: 1) a sand and gravel quarry in Southern New Hampshire, 2) a welding facility, 3) a ceramic powder facility, 4) a metal powder facility, and 5) the air in Cambridge, Massachusetts. A portable sampling box incorporating two Sierra personal cascade impactors was designed and produced. At a 2-LPM flow rate, the effective cut-off diameters of the Sierra impactors were 10.0, 6.0, 3.5, 2.0, 0.9, and 0.5 mm, respectively, for each stage, and the air passed through a PVC filter after the last stage. In addition, the sampling box contained a total respirable dust filter. The substrates were weighed carefully before and after collection, and the weight increases were recorded. All films, foils, and filters were put in the balance chamber for more than 8 hours to stabilize the temperature and humidity before weighing. A Cash 21 microbalance (±0.001 mg) was used to weigh the Mylar films, SS films, and filters."

Jan 1, 1995

By J. Nelsen

Social License to Operate continues to emerge as a paradigm fundamental to sustainable mineral resources development. It relates to a set of concepts, values, tools, and best practices that represent the view of reality held by each participant (industry and stakeholders). It enables a forum for negotiation through which all parties involved are heard, understood and respected. The exploration phase of the mine life cycle is generally the initial contact with local communities and forms the cornerstone for effective communication and trust building. This paper presents the development of a mineral exploration template, based on field studies, aimed at assisting exploration personnel to establish effective communications and engagement with local communities within a process to build a Social License. The template aims to provide practical assistance in conducting situational analysis - a process of examining and mapping a situation, its elements and their relations - in order to acquire knowledge and understanding of a project, locale and region. Situational analysis is the ?full situation of inquiry?, going beyond just ?knowing the subject? to being fully situated. It forms the initial foundation for a planning process presented in the paper that prepares for a strategy of initial engagement and consultation.

Jan 1, 2007

By W. H. Bassett

ABOUT twenty-five years ago the copper industry had outgrown the Lake Superior production. The electrolytic copper producers had- their process well in hand and the industry was well started in the use of electrolytic copper, both for working into wire and sheet, and for the making of alloys. The trade knew that some brands of copper made good material and that other brands did not. There came together a group of men-a refiner, a rolling-mill man, and a manufacturer of electrical machinery-each of them desiring to know why there was a difference and what .it was. Much investigation was done with wire bars and cakes which involved ' chemical analysis, metallographic examination, physical tests and electrical tests. To find the influence of different kinds of copper in making brass 200 carloads of ingots were made into ,eyelet brass under careful observation; each carload was carefully analyzed and tested before used. The alloy was made and rolled under standard conditions and the results recorded.

Jan 1, 1924

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

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

Jan 10, 1987

By N. B. Prenn

Many recent projects have been developed and mined based upon imaginary reserves. Imaginary or "geo-fantasy" reserves have cost many companies millions of dollars. The significant and varied technological advances in the mining arena have done little to curb reserve calculation problems. These advances may have served to increase the prevalence of reserve problems by putting more tools in the hands of personnel not properly trained or experienced in reserve estimation. This paper will discuss potential problems related to calculating reserves, which include: SAMPLING ERRORS -obtaining a representative sample. -Decay analysis of wet RC drilling -Bias introduced by sample preparation. -Assay bias; Coarse gold. -Specific gravity samples. GEOLOGICAL ERRORS -structural boundary definition. -Separation of high grade and low grade zones. -Incorrect continuity assumptions MODELING ERRORS -Too much influence of high grade samples. -Improper method of calculation. -Assumed range of influence. -Errors in the data. ENGINEERING ERRORS -Inaccurate assumptions. -Inadequate metallurgical testing. -Designed bottlenecks.

Jan 1, 1992

Mining has traditionally been seen as a relatively risky business. Exploration is expensive and the failure rate is high, while start-up costs commonly run to many hundreds of million dollars. To justify the risk and expense, the returns from successful projects must also be high. Accurate estimation of mineral resources and reserves is fundamental to the processes of optimising returns and managing risk in mining developments. Systematic development from initial exploration through a series of development stages enables expenditure and risk to be balanced. As a consequence, large scale modern mine are more akin to manufacturing businesses. This paper examines the life history of a mineral reserve estimate from the perspective of a large mining company against a background of modern estimation and reporting requirements. It will also touch on what constitutes, in the author's opinion, `best' practice in reserve estimation, and how this concept can be spread more widely. [Abstract to be re-written and reduced to 100 words]

Jan 1, 2003

By Kenneth M. Haley

THE development of Reserve Mining Company's pelletizing practices began in 1949 in a pilot plant located at Ashland, Kentucky. Three years later the participating companies -Armco Steel Corporation, Cleveland Cliffs Iron Mining Company and Republic Steel Corporation concluded the work and the small plant was shut down. Oglebay Norton and Company was the operating manager and Armco Steel Corporation furnished labor and services. Arthur G. McKee and Company were the designers and constructors. The equipment flow sheet that was employed in that pilot operation is shown in Figure 1. This was Reserve's first step from the laboratory work per-formed by the Mines Experiment Station of the University of Minnesota in developing a commercial operation. The plant used New York state concentrate, purchased from Jones and Laughlin Steel Corporation's Benson Mine. This concentrate, minus 14 mesh, 64.3 per cent Fe, 5.7 per cent Si02, was processed at Ashland to 80 per cent minus 325 mesh in conventional ball mill, classifier, thickener, filtering

Jan 1, 1961

By Charles L. Allie

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

Jan 12, 1980

By H. E. Stamm, H. H. Spain, R. C. Rumble

This paper presents a reservoir analyzer study of the performance of the Woodbine formation in the East Texas basin. The study was made possible by the compilation of available information on the configuration. thickness and the pressure drawdown of the formation. The investigation was made of the pressure variations in the basin incident to production from the several Woodbine reservoirs. In addition, the apparent compressibility of Woodbine water was evaluated so that the potential water yield of the formation could be determined. The distribution of permeability of the Woodbine sand and the interference between producing areas were also investigated. In conjunction with this basic study of the Woodbine sand, an acceptable match of the production-pressure relationship in the East Texas Field was established on the analyzer. INTRODUCTION When production is taken from a reservoir contiguous to an aquifer, the resultant pressure gradient causes a water influx into the reservoir. For some years it has been the practice to predict the performance of water drive reservoirs with mathematical equations that relate the water influx into reservoirs to their pressure behavior. A rigorous solution of these equations is quite complex: hence, the conventional method requires the simplifying assumptions that the formation have constant thickness, permeability and poiosity, and that it have a known and regular shape. It is necessary to use performance records of the reservoir in question to determine the equation constants that make the results of the mathematical calculations duplicate the reservoir history. The values of the constants determined in this manner sometimes differ so greatly from theoretical or measured values that the constants appear unreasonable. The equation constants comprise such factors as the thickness, porosity and permeability of the formation, and the viscosity and compressibility of water. Usually the thickness, porosity and viscosity are known within a reasonable degree of certainty. whereas the formation permeability and the water compressibility are known with considerably less certainty. Hence. permeability and compressibility values are arbitrarily selected that will make the mathematical equations duplicate the reservoir behavior. In many reservoir studies the value of the compressibility of water thus selected has exceeded the reported value by a factor of from severalfold to manyfold. Whether this excessive value results from the inadequacy of the mathematical approach to take into account the shape irregularities and the heterogeneities of the aquifer, or is due to some type of formation compaction and decrease of porosity has not been established. Because of the significance of the potential water yield of formations with reference not only to oil reservoirs but to water resources as well, it is very important that a representative value for the effective water compressibility be determined. One of the best known examples of a field for which it has been possible in the past to predict the future reservoir pressure is the East Texas Field. Until about 1943 it was possible by means of mathematical equations to duplicate past pressure behavior and to predict with a high degree of accuracy the reservoir pressure of the field. Since 1943, however, the calculated pressures have deviated progressively from the observed pressures to such a degree that in recent years, it has

Jan 1, 1951

By H. H. Spain, R. C. Rumble, H. E. Stamm

This paper presents a reservoir analyzer study of the performance of the Woodbine formation in the East Texas basin. The study was made possible by the compilation of available information on the configuration. thickness and the pressure drawdown of the formation. The investigation was made of the pressure variations in the basin incident to production from the several Woodbine reservoirs. In addition, the apparent compressibility of Woodbine water was evaluated so that the potential water yield of the formation could be determined. The distribution of permeability of the Woodbine sand and the interference between producing areas were also investigated. In conjunction with this basic study of the Woodbine sand, an acceptable match of the production-pressure relationship in the East Texas Field was established on the analyzer. INTRODUCTION When production is taken from a reservoir contiguous to an aquifer, the resultant pressure gradient causes a water influx into the reservoir. For some years it has been the practice to predict the performance of water drive reservoirs with mathematical equations that relate the water influx into reservoirs to their pressure behavior. A rigorous solution of these equations is quite complex: hence, the conventional method requires the simplifying assumptions that the formation have constant thickness, permeability and poiosity, and that it have a known and regular shape. It is necessary to use performance records of the reservoir in question to determine the equation constants that make the results of the mathematical calculations duplicate the reservoir history. The values of the constants determined in this manner sometimes differ so greatly from theoretical or measured values that the constants appear unreasonable. The equation constants comprise such factors as the thickness, porosity and permeability of the formation, and the viscosity and compressibility of water. Usually the thickness, porosity and viscosity are known within a reasonable degree of certainty. whereas the formation permeability and the water compressibility are known with considerably less certainty. Hence. permeability and compressibility values are arbitrarily selected that will make the mathematical equations duplicate the reservoir behavior. In many reservoir studies the value of the compressibility of water thus selected has exceeded the reported value by a factor of from severalfold to manyfold. Whether this excessive value results from the inadequacy of the mathematical approach to take into account the shape irregularities and the heterogeneities of the aquifer, or is due to some type of formation compaction and decrease of porosity has not been established. Because of the significance of the potential water yield of formations with reference not only to oil reservoirs but to water resources as well, it is very important that a representative value for the effective water compressibility be determined. One of the best known examples of a field for which it has been possible in the past to predict the future reservoir pressure is the East Texas Field. Until about 1943 it was possible by means of mathematical equations to duplicate past pressure behavior and to predict with a high degree of accuracy the reservoir pressure of the field. Since 1943, however, the calculated pressures have deviated progressively from the observed pressures to such a degree that in recent years, it has

Jan 1, 1951

By C. Perez-Rosales

A simplified method for evaluating specific surface of porous media is described. The theory takes as its point of departure a previous statistical method. Due to its simplicity, it is considered that the proposed method has certain advantages over methods currently used. By making some simplifications, the procedure for estimating specific surface is reduced to superposing, on an enlarged reproduction of a section of the sample to be analyzed, an evenly spaced grid and counting the number of intersections of the grid lines with the perimeter of pores. An indirect test of the theory is presented and the advantages of the method are discussed. INTRODUCTION The specific surface of a porous medium, defined as the surface area of the pores per unit bulk volume, is an important parameter with regard to flow of fluids through porous materials. In particular, it has long been recognized that specific surface is closely related to permeability. Since the internal surface of any porous material is extremely complex, the measurement of specific surface can only be accomplished by indirect means. It is considered that the best method' in current use for determining specific surface is the statistical method of Chalkley, Cornfield and Parka who originally applied it for estimating volume-surface ratios of morphologic components, such as cells or nuclei. This method will be outlined later. A method which has been used extensively in determining specific surface is based on the adsorption of a vapor by the internal surface of porous media. However, this method does not yield the surface pertinent to fluid flow, since it includes the tiny molecular interstices which do not influence the permeability of the porous material.' This method has been employed by Brooks and Purcell measure internal surface area of sedimentary rocks. The Kozeny equation, which relates the specific surface to the permeability of a porous medium, has found wide use in determination of specific surface. However, the KOzeny equation is not strictly correct. Large discrepancies have been found between theory and experiments.' Consequently. the specific surface estimated by this method is subject to serious uncertainties. Finally, there are other methods for specific surface measurements which are based upon a variety of physical and chemical properties: the heat conduction by a gas, ionic adsorption phenomena and direct chemical reactions of the surface with other substances. This paper describes a simplified statistical method for determining specific surface of porous materials. It is considered that, due to its simplicity, it has certain advantages over the methods presently used. For its validity, it is required that the porous structure of the material under study be of a spatially random nature. This is not a serious limitation since the naturally occurring porous media have this characteristic. THEORY The method proposed in this paper takes as its point of departure the interesting result obtained by Chalkley et al The method shows that if a line of length r is dropped repeatedly at random over a section of a porous medium, and counts are made of the number of times the two end points fall within pore areas (denoted by h for "hits") and the number of times the line intersects the perimeter of pores (denoted by c for "cuts"), then in a very large number of throws where S is the internal surface of the sample (surface of pores) and V, the volume of the pores (Fig. 1). Since the specific surface S, and the porosity 4 are given, respectively, by S. = Sl V, and 4 = V,l V,, where V, is the total volume of the sample, it follows from Eq. 1

By G. W. Thomas

A simplified mathematical model of underground conduction heating in a system of limited permeability is presented. The model applies to underground retorting of oil shale, or to reservoirs containing extremely heavy oils. We assume that heat is introduced at a constant rate into a horizontal fracture which communicates between wells. The radial temperature distribution along the fractured surface is approximated by a step-function. Heat transfer away from the fracture is assumed to be by vertical conduction, and all convection effects are neglected. The model also takes into account the possible temperature dependence of thermal conductivity. A general expression for calculating the growth of the step-function temperature distribution with time is derived. The use of this expression and solutions to the one-dimensional heat equation make it possible to construct isotherms. Expressions for calculating .oil recovery, well spacing and heat efficiency are also given. An example calculation is presented for the conduction heating process in oil shale. Finally, the effect of the heat transfer coefficient between the gas and the fracture boundaries is investigated. INTRODUCTION Thermal recovery processes of oil recovery fall into the four general areas of hot fluid injection, forward combustion, reverse combustion and conduction heating. The first three of* these processes have been rather extensively studied in the past decade from both the experimental and theoretical points of view. As a result, it is possible to make reasonable engineering predictions and analyses of these processes. Little attention has been devoted, however, to the conduction heating process other than to note its possible utility.I To a certain extent, conduction heating cannot be divorced from the other regimes cited above, insofar as these provide the source of heat energy. In the conduction heating process, heat is introduced (either by combustion in the forward or reverse mode or by hot fluid injection) into a small fraction of the total reservoir thickness. This fraction may be either a streak of high permeability or an interwell fracture. Heat penetrates by conduction into the adjacent, less permeable regions of the oil-bearing rock, where the direction of conduction is essentially perpendicular to the streak or fracture. The heated product then drains by gravity or is gas driven to production wells. Conduction heating is probably most applicable to systems containing immobile bitumens such as tar sands and oil shale deposits and perhaps to low-permeability reservoirs containing highly viscous crudes. The mechanism also acts in combination with other thermal processes where fingering or overriding of a bed occurs. It seems probable that in at least one field test, conduction heating of this type was very influential.2 In this presentation we give a first approximation to some of the quantitative aspects of conduction heating. Marx and Langenheim3 treated a similar problem where they focused attention upon the injection interval, which spanned the total reservoir thickness. In their model, conduction heat losses to the bounding media imposed a practical limit on the calculated heated area. In the present study, however, we shall confine the injection interval to a small fraction of the reservoir thickness and assume it has no heat capacity. We therefore direct our attention to regions outside the injection interval into which the conduction of heat is beneficial. In particular, we will endeavor to locate specific isotherms in the media bounding the injection interval. Furthermore, we will construct our model to allow the thermal conductivity to vary arbitrarily with temperature. Thus the model will be applicable to underground retorting of shale where variations in thermal conductivity may be important. THEORY We consider the rock matrix to have no or, at best, a very limited permeability; consequently a

Jan 1, 1965

By G. B. Thomas

Several methods of analyzing pressure build-up data in wells have been presented by various authors. This paper reviews the theory and method of D. R. Horner and presents example calculations performed on data obtained by testing several different types of wells. These calculations include, (1) graphical estimation of final static pressure, (2) determination of the productive capacity of the pay away from the well bore and. (3) the degree to which the formation adjacent to the well bore has been damaged by completion or other causes. The methods of testing and precautions which should be taken to assure the best data possible are discussed. Limitations and reliability of calculated results are also treated. INTRODUCTION Pressure testing of wells is generally limited to the determination of producing and static mean formation pressures. The so-called "static" pressure determination.. along with PVT, electric log and production data, enable the reservoir engineer to determine, within reasonable limits, the drive mechanism of the reservoir and in some cases. the amount of edge water encroachment. Producing pressure tests enable calculation of productivity indices and allow the engineer to plan the systematic production of a pool for optimum conservation of suhsurface energy. The radial flow formula advanced by Muskat' has been based on the assumption of incompressible radial fluid flow. It has been known that reservoir fluids do not behave in an ideally incompressible manner. For example, incompressible flow theory indicates a simple logarithmic relationship between the difference of the instantaneous and static well pressures when plotted against time. 'The latter stages of this type of plot of pressure build-up data generally show a marked devia-, tion from the earlier straight line trend. which deviation may be shown to be due to the compressible flow of fluids toward the well bore. Shut-in times of 24, 48. 72. or at most 96 hours are currently in wide use for determining so-called "static" reservoir pressure. Due to the continuation of compressible flow of fluids into the well bore long after this arbitrarily taken shut-in time, the determination of static pressures has almost invariably resulted in lower than equilibrium values. Materials balance calculations made early in the life of a reservoir often result in a calculated reserve which later observations prove to be too low. Failure to obtain reliable "static" reservoir pressures within the prescribed 24 or 48-hour build-up period has undoubtedly been a major factor in obtaining these low estimates. Comparison of theoretical and actual productivity performance has indicated that formation damage and not being able to attain true static pressures have been partially responsible for the observed discrepancies. Research into the theory of compressible flow behavior has resulted in methods of applying these theories to the testing of wells. Application of the developed theories to pressure build-up performance indicates that in many cases the following can be estimated. 1. The static reservoir pressure. 2. The in-place effective formation permeability away from the well bore. 3. The degree or extent by which the formation has been damaged adjacent to the well bore, either through completion methods or subsequent damage due to fluid entry. REVIEW OF THEORY Horner2 has shown that the pressure build-up within a "point source" well is approximated by the following formula: Equation (1) is the approximate "point source" solution to the radial compressible flow equation advanced by Muskat. The solution assumes the following: 1. A point source well is producing at constant rate from the center of an infinite reservoir with a constant pressure at its external boundary. 2. The fluid flowing is present in one phase only. 3. The compressibility and absolute viscosity of the fluid remain essentially constant over the range of temperature and pressure variation encountered. 4. The well i.; shut-in at the sand face arid there is 110 after production into the well bore. 5. The formation permeability is homogeneous in the direction of flow. From Equation (1) it can be seen that if an ideal well were shut-in while producing from a reservoir under the conditions assumed, the pressure build-up would be a logarithmic func-

Jan 1, 1953

By C. R. Dodson, D. Goodwill, E. H. Mayer

The paper summarizes the historical background of the pressure-volume-temperature analyses of reservoir fluids, the errors involved in both the sampling and testing of reservoir fluids, the type of information required of a PVT determination, and the field conditions that limit the application of any one analysis. Particular emphasis is placed on the necessity for approximating as closely as possible the liberation sequence occurring in the producing formation, flow string, and surface separators. A combined differential and flash or "composite" liberation is suggested as the best means of approximating this liberation sequence. HISTORICAL BACKGROUND Petroleum reservoir engineering commonly is considered to be one of the newest fields of petroleum science, yet much work of a theoretical and intuitive nature was done many years before the modern techniques of reservoir analyses were developed. The period 1910-1924 saw considerable work in the field of reservoir behavior done by the U. S. Bureau of Mines. This work:.' although entirely theoretical, pointed out the importance of gas in the recovery of oil from the reservoir. The statement by J. 0. Lewis that 20 per cent or less of the oil originally in place in a pool was recovered from the ground under uncurtailed production conditions caused some operators to re-examine and modify their production methods. However, though these works showed methods of production that have subsequently resulted in greater oil recoveries, they in general went unnoticed, even though there was a fear that the nation's petroleum resources were being exhausted. The period 1924-1933 saw the industry take considerably more interest in reservoir behavior because of an unfounded fear of its inability to replace oil reserves, the possibility of government regulation, and the energy of one man, Henry L. Doherty. Doherty aroused heated discussions in the industry concerning conservation. To prove or disprove his theories the first experimental work on the reservoir behavior of petroleum was undertaken. The papers of Dow and Calkin,' Beecher and Parkhurst," and Mills and Heitheckere are classics although their experimental procedures were crude. These papers proved that the properties of petroleum in the reservoir are quite dissimilar to the properties measured at the well head. Gas dissolved in the oil phase was recognized as having considerable importance as had been predicted by the earlier theorists. Slightly later Miller and his co-workers assembled the first text correlating all the known data on the reservoir behavior of petroleum and used, these data to show the economic value of conservation. Quantitative application of pressure-volume-temperature analyses of reservoir fluids was given in the paper of Coleman, Wilde, and Moore"; this work showed that with sufficient laboratory and production data prediction of reservoir behavior was possible. The refined sampling and labora-

Jan 1, 1953

By F. W. Jessen, G. F. Sears

Cavities in massive salt for the purpose of storage of liquid hydrocarbons have assumed a prominent position in recent years. This paper describes a program to facilitate leaching operations for the formation of specifically shaped storage cavities. Various forms and sizes of cavities may be possible through use of the techniques developed. INTRODUCTION The creation of large under ground storage facilities for natural gas and liquified petroleum gases has been practiced for many years. Use of cavities obtained through solution of salt for this purpose is a fairly new and novel approach, but has gained increasing importance due mainly to the economics of this type of mining operation opposed to hard rock mining or surface and pit storage.1 The idea of controlling the shape of any cavity dissolved from massive salt has not been a prime consideration of companies engaged in the formation of storage space. This has been due mainly to an insufficient knowledge of the mechanics of the leaching process and a dearth of published information dealing with both the desirability and ease of control possible for this type of operation. Two distinct advantages are readily ascertainable. From a stability standpoint (i.e., the ability of the completed cavity to withstand stress imposed by overburden pressure and lateral tectonic stresses) a controlled cavity may be generated which will yield the most favorable attitude to these external forces and remain in operable use for longer periods of time.2 A second advantage is that for a given volume, a sphere (which is one controlIed shape possible) represents the minimum surface area exposed. This may become more important in the future when dealing with refrigeration and product losses in underground cavities. The basic solution mining process, without regard for controlling the final shape, is quite simple in that the equipment and materials required to dissolve a cavity in a salt dome or layered salt section are a source of fresh water, a circulating pump, several strings of tubing and a means of disposal of the return brine. To add control measures involves the use of an inert blanket material above the position at which solution proceeds. Initially the annulus of the largest wash pipe is filled with this blanket down to the top of the proposed cavity. The bottom of the proposed cavity is determined by the depth of the original drilled hole or a blanking plug. The blanket material is added incrementally in stages and displaces the water in the enlarging cavity downward. Where the blanket material has displaced the water, no further solution takes place. The rate at which to add this controlling blanket material so as to be able to form specifically shaped cavities of any particular size is of considerable importance. THE PROBLEM It is desirable to have as much information as possible as to the behavior of the mining operation since visual observation obviously is impossible. It would be advantageous, for example, to have a step-by-step program showing how much salt is to be removed, at what rate and the time required. This type of program would facilitate procurement of surface equipment and provide for proper management of manpower requirements, as well as yield a host of intangible benefits. As mentioned earlier, only rule-of-thumb estimates were available until recently and made this type of pre-planned operation haphazardous at best. Recent work3-' has done much to clear up the ambiguities and inaccuracies attendant to solution mining of salt and has helped to place the operation on a more scientific basis. This study attempts to correlate much of the information thus far developed into a complete mining program; to extend the idea of controlled solution mining to include not only spheres, but solid conic sections of the ellipsoidal variety; and to refine the computation of the rate of removal of salt during the various stages of mining. Dommers3 and IIemson4 made extensive laboratory

Jan 1, 1967

By W. D. Von Gonten, R. L. Whiting

Regression analysis was used to correlate the physical properties of 478 sandstone and 90 carbonate core samples. Porosity, permeability, electrical formation resistivity factor, capillary pressure and the sonic velocity of the shear and compressional waves were measured. Prediction equations for porosity, permeability and electrical formation resistivity factor were found which should be useful in understanding the relationships between the physical properties of porous media in formation evaluation. INTRODUCTION The physical properties of porous media are important to the petroleum engineer and geologist. To evaluate fully the potential and behavior of a subsurface formation such as a petroleum reservoir, certain physical properties of the porous medium must be known. The problem of accurately determining physical properties of subsurface formations has not been solved because many determinations must be made by indirect measurements and because of the difficulties caused by complex pore structure and presence of clays in most naturally occurring porous media. Due to the difficulty in measuring some of the physical properties of porous media, it would be advantageous to be able to predict a certain physical property of a rock from other physical properties of the rock which could be measured more easily and more accurately. Since most potential reservoir rocks are heterogeneous, relationships between the physical properties are very complex and thus far no satisfactory correlations based on theory or laboratory models have been developed. It appears that empirical relationships obtained by measuring the physical properties of a large number of samples of naturally occurring porous media and applying regression analysis to develop equations for one physical property in terms of other rock properties is the best approach. REVIEW LITERATURE The three physical properties used as dependent variables for correlating purposes were porosity, permeability and formation factor. Since formation factor is more difficult to determine, a brief review of the literature is provided on this property. The first work on determining formation factors was published by Archie1 in 1942. He defined this property of a porous medium as F - Ro/Rw......(1) where Ro is the resistivity of the porous medium when completely saturated with a brine of resistivity Rw. Archie found the best correlation between formation factor and porosity was the following equation, where is the porosity fraction and m is the constant characteristic of the rock. The value of m was 1.3 for unconsolidated sand packs, and ranged from 1.8 to 2.0 for consolidated sandstones. In 1950 Patnode and Wyllie2 and De Witte3 observed that the formation factor as determined by Eq. 1 was valid only when the porous medium contained no conductive solids such as clay or shale. When conductive solids are present the formation factor is also dependent on the resistivity of the saturating fluid. Therefore, in samples containing conductive solids the formation factor decreased as the resistivity of the saturating fluid increased. Because of this. the measured formation factor was called the apparent formation factor and was designated Fa. Patnode and Wyllie proposed the following equation,

By L. A. Rapport, J. R. McCord, W. P. Drews, K. H. Coats

Water movement about an oil or gas reservoir can be predicted provided an aquifer influence function is known. This function is generally determined from an electric analyzer study or by fitting an idealized mathematical model to field data. A prohlem arises in determining that influence function which best reproduces field data and satisfies certain smoothness requirements, without making any idealizations concerning aquifer geometry and heterogeneity. This study shows the problem to be solvable through the linear programming technique and presents applications to one oil reservoir and two gas storage reservoirs. In addition, a theorem is stated and proven which confirms influence junction smoothness requirements heretofore accepted on physical intuition and establishes additiotzal smoothness constraints formerly not recognized. INTRODUCTION The pressure decline accompanying production results in water encroachment into oil and gas reservoirs situated adjacent to aquifers. The importance of this water movement derives from the significant dependence of production rate upon reservoir pressure and of pressure, in turn, upon water encroachment. Pressure and rate of water movement are uniquely related through an "influence" function which may be in the form of reservoir pressure response to a unit water encroachment rate or of water movement caused by a unit reservoir pressure drop. This influence function reflects the heterogeneity and geometry of the aquifer and is therefore particular to each reservoir. If the influence function is known, then reservoir pressure may be predicted from a production rate schedule, or vice-versa, with the aid of the material balance equation which relates water influx to reservoir production and pressure. The problem of determining this influence function for a given reservoir has been approached in three different ways. Van Everdingen and Hurst calculated and tabulated certain influence functions pertinent to idealized mathematical models which satisfy simplifying assumptions such as aquifer homogeneity and elementary reservoir-aquifer geometries.' Other studies developed additional tables for a variety of elementary reservoir-aquifer geometrical con- figurations and demonstrated their validity and utility.:' Hicks, et al. described the use of an electrical network to determine an influence function giving "best" match of field history and allowing prediction of future performance.' Finally, Hutchinson and Sikora5 and Katz, Tek and Jones attempted to develop a calculation method for deriving the influence function directly from field data with no idealizing assumptions concerning aquifer geometry and homogeneity. They encountered difficulties with inaccuracies in field data and with satisfaction of well-known smoothness requirements of the influence function. This study is concerned with the determination of influence functions from field data and, more specifically, treats the following problem: given field pressure-production data of arbitrary accuracy, how does one determine the influence function which, subject to satisfaction of known smoothness requirements, gives best agreement with field history? A rigorous solution to this problem is ohtained through linear programming as described below. Example applications to one oil and two gas storage reservoirs are presented. In addition, a theorem is proved which confirms influence function smoothness requirements heretofore accepted on physical intuition and establishes additional smoothness requirements formerly not recognized.

Jan 1, 1965

By F. M. Stewart

A paper by Hurst and van Ever-dingen in 1949 led to the practical solution of many nonsteady-state flow problems.' Subsequently, applications of this material have been discussed by several authors. The purpose of this note is to outline a simplified direct procedure for water-influx pressure predictions for pressures above the bubble point. The prcposed procedure is unique in that after the "time-conversion factor" and the "influx constant" are obtained, predictions are made without trial and error. The elimination of trial and error is made possible by expressing certain PVT characteristics as explicit functions of pressure. The procedure is applicable to fields producing at pressures above the bubble point and under radial water drive from an infinite aquifer. In the usual applications it is required to predict (1) withdrawals given future pressure history or (2) reservoir pressure, given future withdrawals or production. The theory, assumptions and conditions under which the nonsteady-state solutions hold are covered in the references and are not repeated. The simplified direct procedure outlined here has been developed from published solutions for the so-called "pressure case". A similar direct procedure for the same special conditions has been developed for the "rate case"; however, in this instance the procedure is lengthy and offers no advantage over the trial-and-error solution. Consequently, this paper covers only the procedure based on the pressure case. DISCUSSION Under the pressure case prediction of withdrawals from a water-drive field, given the future pressure history, involves the use of a nonsteady-state water-influx term together with the material balance equation. The solution of a particular problem requires the following steps: (1) division of past pressure and production data into equal time intervals; (2) computation of cumulative water influx vs pressure drop for each time interval of past history using the material balance equation; (3) determination of the best values for the "time-conversion factor" and the "influx constant" by equating the nonsteady-state influx term with the influx computed in Step 2 (this requires a trial-and-error procedure); (4) prediction of influx for a given pressure drop using the influx term and the computed factors from Step 3; and (5) calculation of fluids produced using the predicted influx in the material balance equation. A forecast of the water-oil ratio or water production rate is required to complete the calculation. It is assumed that a reliable estimate of original oil in place is available, as well as the production and pressure history, and data on reservoir fluid characteristics. If the original oil in place is unknown, it may be computed simultaneously with the tinie-conversion influx factors2,3. However, the oil in place so computed may be considerably in error3; consequently, it is usually preferable to use oil in place estimated from pore-volume data unless there is considerable doubt as to its accuracy. The Material Balance Equation Above the Bubble Point Above the bubble point the simpli-

By R. C. McFarlane, T. D. Mueller, J. E. Walstrom

In evaluating uncertainty, experiments are usually performed repeatedly and then conclusions are drawn from the distribution of results. With the advent of high-speed electronic computers, it is possible to perform experiments using mathematical models constructed to simulate complex experiments or operations. Statistical methods are then applied to the results of the simulated experiments. This procedure forms the busis of this paper. Demonstrated is the need for properly accounting for uncertainty in petroleum engineering problems. How uncertainty affects solutions is evaluated in three example illustrations. The method used to evaluate uncertainty in petroleum engineering studies is the Monte Carlo simulation procedure.'-" INTRODUCTION The solution to most technical problems may be derived from interrelationships among several quantities called variables or parameters. There may be only a few variables or several hundred. Interrelationships among parameters may be explicit or implicit, well established or only approximate. Some variables that fully or partially depend on the magnitude of others are called dependent variables. Input variables for most practical problems are not precisely known; there is usually an uncertainty in their value. The degree of uncertainty may vary from one variable to another. Variables that are known accurately are called determinates.' For instance, the gravity of crude obtained from a particular pool may be known precisely, and therefore is a determinate. The degree of precision with which a quantity can be determined increases as data describing the pool are accumulated during the development of the field and the producing life of the pool. The uncertainty of a parameter may result from difficulty in directly and accurately measuring the quantity. This is particularly true of the physical reservoir parameters which, at best, can only be sampled at various points, and which are subject to errors caused by presence of the borehole and borehole fluid or by changes that occur during the transfer of rock and its fluids to laboratory temperature and pressure conditions. Uncertainty may also result in attempting to predict future parameter values. This type of uncertainty is particularly evident in investment analyses involving future costs, prices, sales volumes and product demand. Uncertainty in the solution to investment problems is often called risk, and its study is called risk analysis.' Uncertainty also enters into biological and sociological analyses in which indeterminate factors are often important due to limited control of the experimental material. It is customary, in evaluating uncertainty, to perform repeated experiments and to draw conclusions from the distribution of the results of these experiments. With the advent of the high-speed electronic computer, it is possible to construct mathematical models which simulate complex experiments or operations and to perform the experiments repeatedly, utilizing the models. Statistical methods are then applied to the results of the simulated experiments This method forms the basis of the investigation reported here. PROBABILITY DISTRIBUTIONS FOR VARIABLES The uncertainty in the value of a variable may be indicated by a probabilistic description accomplished by expressing the quantity by a probability distribution. Many recognized probability distributions can be used to describe physical quantities. Recent studies used various types of distributions to describe core analysis data.',' However, for the examples in this paper, the uniform and triangular distributions are believed to reasonably approximate the data used (Fig. 1). The uniform distribution confines the variable between an upper and a lower limit. The variable may lie anywhere between the two limits. This distribution is used when no one range of values for a variable is more probable than any other, but information or intuitive reasoning indicates the variable will lie somewhere between the chosen limits. The triangular distribution is used for a variable when more data are available to indicate a central tendency of distribution. This allows postulating a "most likely" value to the distribution and upper and lower limits. In this case, as for the uniform distribution, the variable is not expected to assume a value less than the lower limit or greater than the upper limit. However, with improved quality of data it can be postulated that the variable will tend to assume a value close to the most likely value, and that there will be a decreasing probability for values away from the most likely value. The area under either of these probability distributions is equal to unity since it is assumed that there is a 100 percent probability that the variable will lie somewhere under the curve. An ordinate erected at any particular value of the variable divides the area under the curve into two parts: the area to the left of the ordinate represents the probability that the value of the variable will be equal to or less than the value of the variable at the position of the ordinate, and vice versa. The probability is zero that the variable will have any specific deterministic value. If two ordinates are drawn for any two values of the variable, the probability that the variables will have a value lying between these ordinates is equal to the area under the curve lying between the ordinates.