Search Documents

By Yi Luo, Jianwei Cheng, Mingming Li, Lei Xiao

Concealed thermal events in longwall gobs present safety hazards for underground coal miners and can significantly interrupt mining production. The probable causes for such thermal events range from coal oxidation at slightly above ambient temperature to smoldering fires due to spontaneous combustion of the broken coal left in the gobs. Locating such “hot spots” is the key step for planning and implementing mitigation measures to bring such thermal events under control. If a mine section, or even an entire mine, has to be sealed, understanding the status of the sealed mine atmosphere is important for the decisions guiding subsequent actions. The gas composition data obtained through the tube bundle mine atmospheric monitoring system is firsthand information of great value. Fire ratios proposed by various researchers and derived from gas composition data are useful in determining the fire status and even approximating the source locations (Timko and Derick, 2006). However, the applicability of each of these ratios under particular conditions should be carefully evaluated. The use of properly combined fire ratios can enhance the certainty of analysis. The explosibility of a sealed mine atmosphere deserves primary concern for safety reasons during seal construction, idle time and reopening of the sealed area.Mine atmospheric monitoring data from a longwall mine with flame fire have been analyzed. Efforts have been made to identify the possible causes, locations and statuses of these events using the obtained gas data, as well as other ventilation, geological and mining parameters. The effect of injected inert gas to control the events has also been quantified. The explosibility of the sealed atmosphere has been determined using two methods.

Jan 1, 2012

By T de Wit, S Meyer

The location of seismic events in mines is traditionally calculated based on the arrival times of P- and S-waves for a given velocity model (homogeneous or heterogeneous). This technique is successful when the sources of seismic events are sudden and localised (eg episodes of fracturing in solid rock or dynamic slip along a plane of weakness) and thus the radiated seismic waves have clear arrivals. In massive mining operations such as block caves, strong seismic signals often originate from the movement of loose rocks on the perimeter of the fracture zone. These signals are easily recordable but do not contain easily identifiable direct P- and S-waves, making the estimation of location through traditional methods more complicated. This problem has been approached by using a ‘brightness’ location method that estimates the origin of seismic waves based on travel times of envelopes of seismic energy. This tool is extended to be applicable in a heterogeneous velocity environment that now accounts for the presence of a cave or other strong velocity perturbations. The suggested method is particularly important when caving occurs in low-stress environments (eg close to the surface).CITATION:Meyer, S and de Wit, T, 2016. Location of seismic events associated with caving using brightness method, in Proceedings Seventh International Conference and Exhibition on Mass Mining (MassMin 2016), pp 225–236 (The Australasian Institute of Mining and Metallurgy: Melbourne).

May 9, 2016

By D. B. Rebuli

Seismic networks in deep gold mines of South Africa are often planar owing to the tabular geometry of the orebody and, therefore, the mining infrastructure. Consequently, seismic event locations are poorly constrained in the direction perpendicular to the reef plane, and symmetrical locations above and below the reef plane may fit the observed data equally well. This problem may be resolved in situations where the strata above and below the reef plane have significantly different seismic velocities, which is the case for the Ventersdorp Contact Reef at Kusasalethu mine (previously Elandskraal, the amalgamation of Deelkraal and Elandsrand). A location method using a layered location model that takes into account the velocity contrast between the layers was tested. Observed direct- and head-waves were used to determine the velocity bounds of the P-wave velocities in the two-layer model. The data-set was also used to find the best P- and S-wave velocities for the two layers by grid searching over a range of velocities and minimizing the sum of the travel-time residuals. The results from the grid search and those from the calibration of velocities using direct- and head-waves give similar P-wave velocities for each of the two layers. A second set of seismic events was then relocated with this layered model and locations closer to reef were found. This simple layered location method goes some way to reduce the reefperpendicular location errors for Kusasalethu mine in particular. Keywords mine seismology, velocity model, location methods.

Jun 1, 2013

By D. Yu Pleev, V. G. Krivolapov

Mines of Southern Kuzbass working out, basically, coal formations of flat seams. Seams, as a rule, are inclined to self-ignition. Detection at early stages of the centers of endogenous underground fires has great value for safe conducting mountain works. One of methods of detection of an early stage of underground fires in the worked out space is the method of superficial radon survey. In the report the method is de-scribed and examples of superficial radon survey on concrete layers are given.

Jan 1, 2003

By Michel L. Bilodeau, Brian W. Mackenzie, Dallas W. Davis

"The paper formulates an evaluation procedure for examining the effects of location on the economics of developing mineral deposits. This methodology is applied to assess locational aspects of Nova Scotia's competitive position f or mine development relative to other Canadian regions and provinces. The appraisal concerns two representative mine development opportunities, a medium-size open-pit tungsten deposit, and a small underground gold deposit. The influence of regional economic conditions is initially assessed. Then, the consequences of taxation are isolated. Finally, the combined impact of regional and taxation conditions is examined. SettingThe economic potential of mining may be assessed in a number of ways. One approach, increasingly used for both corporate planning and government policy purposes, is to evaluate the competitive position of an actual or potential mineral producer in the international marketplace relative to other suppliers.Such examinations of competitive position vary from cost league comparisons of breakeven costs for operating mines, through the assessment of medium-term costs associated with developing new mine capacity, to long-term evaluations of the economics of mineral supply from the start of exploration. The complexity of the analysis increase s with the scope and timeframe of interest. Many factors are involved. One component of competitive position is the quality of the mineral resources available. Other key factors include existing technology, wage and salary rates, labour productivity, the cost of capital, infra structure requirements, distances from sources of materials supply an d to mineral markets, management attributes, organizational efficiency, government mineral policies (particularly taxation), and inflation and exchange rates. The paper fits within this competitive position framework.The variable economic conditions which exist within Canada have a direct bearing on the viability of mine development and, thereby, on mineral exploration potential. The paper examines locational aspects of Nova Scotia 's competitive position."

Jan 1, 1986

The Charters Towers goldfield in northeastern Queensland produced over 206 tonnes of gold from gold-quartz veins hosted by Mid-Palaeozoic granodiorite, tonalite, diorite and to a lesser extent amphibolite-grade meta-sedimentary schist and meta-igneous rock. Mineralisation occurs in grouped, moderately dipping shear zones with altered wallrock rafts and a complex network of quartz veins.Over 75 per cent of the production came from three main vein systems in the central goldfield: the Brilliant, Day Dawn and Towers Hill vein systems.The Brilliant vein system is hosted by undeformed tonalite and contains north-plunging oreshoots on a north-dipping main fissure and on shorter east-dripping hangingwall fissures. The Day Dawn vein system is localised along a north-dipping dyke at the contact of two deformed granitoids in a brittly reactivated mylonite, where horizontal oreshoot plunges may be related to the original megascopic undulations in that structure. The Towers Hill vein system is hosted in two granodiorite types and has multiple ore controls.Structural analysis indicates that the distribution of the oreshoots within the vein systems is variously controlled by the primary fabric and igneous contacts of the granitoid host rocks. Splaying, changes in dip, and fissure intersections are the main controls. The largest vein systems are both in reactivated shallow-dipping mylonite, or in fissures in younger undeformed granitoids. Each vein system lies in a separate structural domain, usually where older east-striking structures are preserved.

Jan 1, 1990

By Maribel Guzmán Córdova

Ante la problemática del tratamiento de drenajes ácidos de mina (DAM) se ha previsto solucionarlo a través del uso de tecnologías convencionales existentes como la precipitación química, osmosis inversa, intercambio iónico y separación por membranas. Sin embargo, estas técnicas presentan significativas desventajas como son los requerimientos de reactivos y costos de energía; además de la generación de lodos tóxicos. En este sentido, el presente texto describe el uso de “lodos activados”- en tanto materiales de desecho que procedentes de un proceso industrial, pueden remover metales pesados de un drenaje ácido de mina,

Sep 14, 2009

By A. R. Jr. Amos

LOG washers have been used for many years in the washing of clay iron ores, phosphate rock and manganese ores, but not until the past 15 years have they been employed to any extent in the preparation of aggregate materials and fluxing stone. For a machine of such importance where tenacious clay or soft rock must be eliminated, there are few data generally available. The purpose of this paper is to set forth the results obtained from several operations, together with a few suggestions for planning a specific installation. Where the material to be washed out consists of loose surface soil, silt, or small amounts of quickly soluble clay, there is little competition with other forms of washers, because greater power and maintenance are required for the operation of log washers, but where the foreign matter is a stiff clay, such as occurs with high-calcium limestones and in marginal gravel deposits, the ability of other kinds of washers to eliminate com-pletely the hard clay balls is at least uncertain. Then the log washer comes into its own field, for with the breaking and churning action of the logs it will clean out the toughest clay. Even with material from the hydraulic dredge, if it contains tough clay, the log washer is desirable unless the gravel is subjected to a costly soaking and rehandling process. Recently there has been a considerable development in other types of machines to accomplish these results but it is not within the scope of this paper to make comparisons between competitive types, but to present the best obtainable data on log washers only. From the primitive design that gave it its name, a long log hewn from a tree, with paddles forged to shape from strap iron and attached by lag screws, the log washer has gradually been developed to the present all-steel type. As now manufactured, the log washer consists essentially of two built-up logs on which are mounted a series of paddles, carried by shafts fitted at the ends of the logs on bearings fastened outside the enclosing steel box, or trough, which is set at an angle of about 5° to 10° and into which the material is fed. Stuffing boxes at the places where the shafts pass through the box at the lower end are equipped with water-

Jan 1, 1936

By A. R. Jr. Amos

LOG washers have been used for many years in the washing of clay iron ores, phosphate rock and manganese ores, but not until the past 15 years have they been employed to any extent in the preparation of aggregate materials and fluxing stone. For a machine of such importance where tenacious clay or soft rock must be eliminated, there are few data generally available. The purpose of this paper is to set forth the results obtained from several operations, together with a few suggestions for planning a specific installation. Where the material to be washed out consists of loose surface soil, silt, or small amounts of quickly soluble clay, there is little competition with other forms of washers, because greater power and maintenance are required for the operation of log washers, but where the foreign matter is a stiff clay, such as occurs with high-calcium limestones and in marginal gravel deposits, the ability of other kinds of washers to eliminate com-pletely the hard clay balls is at least uncertain. Then the log washer comes into its own field, for with the breaking and churning action of the logs it will clean out the toughest clay. Even with material from the hydraulic dredge, if it contains tough clay, the log washer is desirable unless the gravel is subjected to a costly soaking and rehandling process. Recently there has been a considerable development in other types of machines to accomplish these results but it is not within the scope of this paper to make comparisons between competitive types, but to present the best obtainable data on log washers only. From the primitive design that gave it its name, a long log hewn from a tree, with paddles forged to shape from strap iron and attached by lag screws, the log washer has gradually been developed to the present all-steel type. As now manufactured, the log washer consists essentially of two built-up logs on which are mounted a series of paddles, carried by shafts fitted at the ends of the logs on bearings fastened outside the enclosing steel box, or trough, which is set at an angle of about 5° to 10° and into which the material is fed. Stuffing boxes at the places where the shafts pass through the box at the lower end are equipped with water-

Jan 1, 1936

By Robert E. Bush, E. S. Mardock

The objective of this paper is to present practical methods of applying radioactivity logs to problems of interest both to those engaged in evaluating fundamental reservoir data as well as to those engaged in subsurface structural work. The basic theories underlying the operation of radioactivity logging methods are briefly discussed along with some of the problems arising in the interpretation of the log. Relative methods of quantitative interpretation are described, with a discussion of relative reference lines and specific examples are given showing the application of these relative quantitative interpretation methods. To show the variation of the methods described, specimen wells penetrating formations widely divergent in point of geologic age were selected to demonstrate the application. A relative reference line is determined for the Neutron Curve and by comparison to known petrophysical data such as porosities from diamond cores, a direct quantitative relationship is developed. The application of this relative calibration curve to logs of unknown wells will permit a reasonably accurate determination of porosity in limestone and dolomite formations. Specific examples are given showing the application of these relative methods of interpretation to logs of the Smackover limestone in Louisiana, the Edwards limestone and dolomites of South Texas, and Permian dolomites of West Texas. Radioactivity logs have been used extensively by the petroleum industry for a number of years and recognition of their value has grown steadily, not only as tools to aid in the solution of specific stratigraphic problems, but also as potential sources of the basic data necessary for the estimation of net pay and reserves. The value of radioactivity logs for stratigraphic work has been pointed out in the literature,1,2,4,5,6 and now has long been recognized by the industry, but the quantitative application of the logs has advanced but little beyond the stage of a strong conviction that radioactivity logging curves will some day provide the information on such important quantities as the porosity and the liquid saturation of the formations. During the past few years the staffs of the Lane-Wells Company and Well Surveys, Incorporated, have devoted a great deal of time, thought, and experimentation to the relationships existing between radioactivity curves and the lithologic conditions encountered, and this paper is, in a manner of' speaking, intended as a progress report of this work and includes suggestions which the authors have found useful in the interpretation of radioactivity logs. The section on theory is not intended as a complete exposition of the physical processes involved, but rather as a brief introduction to the reader of these processes so that he may acquire a better understanding of the behavior of the radioactivity curves under the varied conditions encountered in normal well logging. Three relative methods for using neutron logs for the quantitative estimation of porosity in limestone or dolomite are presented. These methods are based on the comparison of logs and porosity data from wells in Louisiana, South Texas, West Texas, and New Mexico. These areas were selected because not only were there more core data available from these areas but also the majority of these data were from limestone or dolomitic zones. The correlation between core data and logs is simpler for limestones and dolomites than it is for sandstones due to the fact that sandstones are more complex than limestones and differ more widely among themselves as to composition of cementing materials, siltiness, etc. When data applicable to the various types of sandstones are available, the authors hope to present a discussion similar to this paper covering sandstones or possibly a more general discussion covering all types of reservoir rocks. The relative methods for estimating porosity which are presented in this paper fall short of the absolute method sought. However, these methods possess the advantage that they can be applied to almost any neutron log as long as the range of application is held within the limits imposed by the local lithologic conditions which prevailed during calibration, i.e., these methods are for short range application only. THE THEORETICAL BASIS FOR RADIOACTIVITY LOGGING The two redioactivity logging curves commercially available today are the

Jan 1, 1950

By M. H. Waxman, M. Gondouin, H. J. Hill

Streaming-potential experiments were conducted within the Muddy- and Dakota-sandstone interval of a Denver basin well. Analysis of the data shows that, for this case, streaming potentials opposite sands were higher than those opposite shales when the mud was fresh. Streaming potentials opposite shales, however, were of sufficient magniiIlde to be important in SP interpretation. They were linearly dependent upon pressure differential, and they increased with mud filtrate resistivity. The magnitude of streaming potentials opposite sands was influenced by the characteristics of the original mud cake formed at the time of drilling. Use of mud samples taken at the time of logging may not always be satisfactory for determining streaming-potential corrections. Both the streaming and electrochemical potentials were found to be affected by the previous history of the borehole system. Following a change in mud characteristics, several days were required for SP stabilization. This was, in part, due to the mud cake acting as an imperfect shale membrane. INTRODUCTION Laboratory investigations by Wyllie,1 Salisch,2 Moore," Sen Gupta and Bannerjee,0 and Hill and Anderson' have shown that streaming potentials developed across mud cakes deposited on permeable filter beds are of sufficient magnitude to be critically important in SP interpretation. Streaming potentials across shales have been reported by Schenck,' Gondouin and Scala; and Hill and Anderson." The magnitude of the shale streaming potentials obtained by Schenck and by Gondouin et a1 led these authors to conclude that in many cases the shale streaming potential would approximately equal the mud-cake streaming potential. Hill and Anderson, however, obtained somewhat lower magnitudes of shale streaming potential. The field experiment described in the present paper was a cooperative project carried out by the respective operations and research departments of Schlumberger Well Surveying Corp., Shell Development Co. and Shell Oil Co. The primary objective of the experiment was to determine in the borehole the streaming-potential behavior of both shales and mud cakes under conditions approaching, as nearly as possible, those normally encountered in logging practice. METHODS AND EQUIPMENT WELL A Muddy- and Dakota-sandstone vie11 in the Denver basin was chosen for the experiment. Both sandstones were water-bearing and, since reservoir pressures in these sandstones are abnormally low (gradient of approximately 0.25 psi/ft), it was expected that studies of streaming potential could be conducted from low-pressure differentials (approximately 200 psi) to about 27000 psi by bailing mud from the hole, filling up and then applying wellhead pressure. The selected well encountered the Muddy sand-

By R. H. Winn

A method for the selection of the most suitable corrosion inhibitor for a particular system is given. The method involves the evaluation of surface passivity by means of copper ion displacement after the metal surface has been subjected to the corrosive atmosphere. The application of the method to the selection of suitable inhibitors for use in producing wells is described. Correlation with recognized laboratory tests has beet1 found to be good within the limitations of the test. The use of the test to evaluate inhibitor persistencc is described. INTRODUCTION Corrosion as it is associated with the industry shows wide variations in its qualitative as well as its quantitative aspects. This complicates the work of the corrosion engineer. In the past few years the efforts of the corrosion engineer have been further taxed by the large numbers of products being marketed as corrosion inhibitors by various manufacturers. These products are submitted to him together with a limited amount of data to support their use as inhibitors. The data are usually authentic and reliable yet all too often when these products are applied to an actual corrosive system the results are disappointing due to the ineffectiveness of the inhibitor and the corrosive system. The corrosion engineer must select a comparatively limited number of inhibitors for field trial from the large number available. The selection is usually made on the basis of a standardized laboratory test such as that adopted by the National Association of Corrosion Engineers or by field experience in an area; however, testing of this type is discouraging on large numbers of compounds, particularly in view of the specificity of environment. It has been an object of this work to devise a method which will allow the corrosion engineer to select the best inhibitor for a specific environment under field conditions. The test would of necessity be rapid and inexpensive. It would preferably eliminate the need for expensive and delicate instruments. STATEMENT OF THEORY The test as developed makes use of the fact that metal ions will be displaced from solution by metals having a more positive electrode potential. Thus, if an active iron surface is exposed to a solution of a copper salt, free copper will plate spontaneously on the iron surface as it is displaced from the salt solution by iron ions. If the iron surface is passified, as by coating with an insulating-type organic corrosion inhibitor, it is no longer capable of displacing the copper ions due to its own limited ionization. The result will be a change in the amount of copper plated upon exposure to a neutral copper salt solution. In practice it has been found possible to measure the magnitude of the alteration of the cell potential by determining the copper deposited on a standard iron surface under controlled conditions. This means that, as the film builds up, the resulting shift in the electrode potential tends to cause the iron surface to act more like a noble metal, resulting in a reduction in the cell potential to a point that copper will no longer plate spontaneously. Obviously, the salts of other metals such as silver, gold or platinum might be used for this operation; however, copper was selected because of the color difference of the plate, the comparatively low cost of copper salts and the ease of evaluation of the plated copper. Furthermore, its high cell potential provides a substantial driving energy. The determination of the copper is readily done colorimetrically by dissolving it in 5 per cent ammonium hydroxide solution containing 2 per cent potassium per-sulfate. The blue color of the copper ammonic complex ion so developed may be compared with standards prepared by dissolving known quantities of copper powder in the same reagent. In routine screening it is seldom desirable to determine the amount of copper plated. It is usually more important to determine the concentration of inhibitor required to prevent the plating of copper.

By G. R. Pickett

Examples are presented which show that the velocity~ amplitude, attenuation and apparent frequency of several acoustic waves can be recorded in the borehole. Examination of such recordings, termed "character" logs, indicates that the wave types observed include a refracted compressional wave and a wave which travels with formation shear velocity. Laboratory data are used to show that compressional and shear wave velocities are dependent on porosity, effective stress and lithology; but that the change in reciprocal velocity per unit change in porosity is larger for shear waves than for compressional waves. We, therefore, conclude that. the accuracy of porosity determinations can sometimes be improved by use of shear wave velocities, provided that the shear wave amplitudes are large enough to delineate the shear arrival from the preceding compressional arrival on the character log. Borehole data are presented which show that the difference between shear wave and compressional wave reciprocal velocities can be used to predict porosities. This is a refinement which may allow the prediction of porosities from single-receiver acoustic logs without introduction of errors from borehole fluid traveltimes. Laboratory and field data are presented to show that the relationship between compressional and shear wave velocities can be used to indicate lithology. An example is presented to show that fractures usually cause a greater reduction in borehole shear wave amplitudes than in compressional wave amplitudes, an effect which may offer a more reliable means of detecting fractures. The complexity of the borehole acoustic wave train can rake presently available cement bond logs highly sensitive to the gate and bias settings used. The character log offers a means to circumvent possible misinterpretations by recording all amplitudes, from which the interpreter can select the appropriate data for evaluating the cement bond. Character logs may also be used as a quality control for open-hole transit-time logs when existence of small compressional wave amplitudes interferes with the proper functioning of bias-controlled timing devices. Evaluation of the potential uses of character log data is not complete; but a character log presented in a form convenient for routine use would be a desirable addition to currently available logs. To summarize, possible applications for such a log in formation evaluation include the following (1) quality control of transit-time logs, (2) refinement of porosity predictions, (3) determination of lithology, (4) improvement of fracture detection and (5) improvement of cement bond evaluation. Suggestions are made regarding the requirements for a suficient but practical character log for routine use. INTRODUCTION Acoustic logs have become a widely used porosity tool in formation evaluation. In addition, there is a growing application of acoustic logs in cement bond evaluation and fracture detection. These applications have mainly involved the use of logs of first-arrival transit times and amplitudes and have not included detailed studies of the complete signal. The purpose of this paper is to show that significant benefits in formation evaluation can be gained by a more complete use of the acoustic wave train generated in the borehole by an acoustic logging tool. We hope that this discussion will also stimulate further development of logs suitable for routine use so that these benefits may be realized. Examples of acoustic wave train logs, termed "character" logs, are presented to show that several identifiable acoustic waves are present in the borehole. The measurable characteristics of these acoustic waves and some of their relations to formation properties of interest are also discussed. The more obvious potential uses of character logs are listed, and some suggestions are made regarding the requirements for a sufficient but practical character log for routine use. CHARACTER LOGS Some 10 years ago, Vogel' and Summers and Broding' noted that the signals received uphole from an acoustic logging tool located in a borehole had a number of interesting characteristics. The logging tool consisted of two or more pressure transducers spaced on an acoustically insulated body (Fig. la). One of the pressure transducers was used as a transmitter to generate pressure waves in the borehole fluid. The other transducer served as a receiver to detect any pressure waves reaching it in the borehole. The receiver then converted these pressure waves to electrical signals which were transmitted to the surface and displayed on an oscilloscope as a record of time vs receiver-signal amplitude. Fig. lb is a schematic representation of a typical record. The interesting characteristics seen in the earlier' and subsequent experiments were (1)

By M. Felsenthal

Engineers are frequently faced with the problem of having to predict oil recovery from a solution gas drive reservoir in the early life of a field. This is often the time when actual laboratory or field-derived gas-oil permeability (kg/k0) data are not yet available. Although excellent collections of k0/k0 data for various types of rocks may be found in the literature,1,2,3 a need exists for a more detailed correlation of k0/k0 data with core properties such as porosity, permeability, clay content or pore-size distribution. It is the purpose of this article to present this type of correlation for sandstone material. The correlation is based on laboratory tests performed on over 300 core samples from 19 reservoirs in the United States and Canada. EXPERIMENTAL Cores were analyzed by routine procedures, using samples of conventional size (7/8-in. diameter, 7/8-in long) cut parallel to the bedding planes. Air permeability was measured at 1.67 atm mean pressure in a Hassler-type apparatus. Porosity was determined with the aid of a Ruska bulk volume meter and a Boyle's Law rock volume meter, An evaporation method' was used to evaluate irreducible water saturations. Core samples were saturated with 25,000 ppm NaCl brine. The brine was reduced to the irreducible saturation, and the remainder of the pore space was filled with kerosene. This was followed by multiphase flow tests using the so-called "dispersed feed" technique.5 n this technique, gas and oil (kerosene and kerosene-saturated air) were dispersed with the aid of a porous Alun-durn feed head, the face of which contained a waffle-shaped system of Lucite-lined grooves. Oil entered the back side of the feed head and emerged at the raised areas of the front face. Gas was introduced directly into the system of grooves at the front face. The gas-oil mixture then entered the test core which was kept in capillary contact with the feed head. Permeabilities were determined from rate and pressure measurements and saturations from weights. All determinations were performed at equilibrium or steady-state conditions. Pore-size distributions were determined on selected samples by a mercury injection technique.6 GENERAL CORRELATIVE TRENDS When examining data from a given reservoir, it was noted that k0/k0 vs saturation curves became generally less steep as specific air permeability increased, a development that has a favorable effect on oil displacement efficiency by solution gas drive. This trend has also been reported by others in the literature.' An effect of porosity on k0/k0 data was also noted. This effect was not generally discernible in a study of relative permeability data for a given reservoir but became apparent when data for sandstone reservoirs of similar lithology but differing average porosity were compared. For instance, in a comparison of argillaceous and/or calcareous sandstones from 11 reservoirs ranging in average porosities from 14 to 28 per cent, a definite trend was noted, indicating that for a given permeability k0/k0 curves became less favorable (i.e., steeper) as porosity increased. A similar trend was observed for a group of comparatively clean sandstones from five reservoirs ranging in average porosities from 15 to 21 per cent. For a given permeability and porosity, comparatively clean sandstones gave more favorable k0/k0 curves than argillaceous and/or calcareous sandstone or chert reservoirs. The laboratory studies also showed that least favorable k0/k0 curves were exhibited by shaly sandstone, conglomerate, and sandstone containing carbon inclusions. Table 1 lists the various source reservoirs classified by lithology types. To facilitate a correlation of k0/k0 data with core characteristcis, it was found expedient to work with

By L. S. Allen, R. L. Caldwell, W. R. Mills

Borehole models used in the study of nuclear well logging are reviewed and the merit of heterogeneous us homogeneous formation simulation is discussed. .A heterogeneous model for simulating sandstone formations of variable porosity is described. Simulation adequacy was checked by steady-state and pulsed-neutron measurements. INTRODUCTION In the course of developing new methods of nuclear well logging, the problem of testing new logging instruments under controlled conditions always arises. This testing is usually done in wells of known lithology, porosity and fluid saturation, or in laboratory test facilities which simulate the borehole environment. Although the former approach is possibly more satisfying, it is frequently difficult or inconvenient to find wells with enough variation in the cited parameters for a truly comprehensive study, On the other hand. laboratow- test facilities can simulate a wide variety of formation and borehole conditions if care is exercised to maintain adequate simulation under all conditions. Workers concerned with the development of nuclear well logging methods have reported several types of research facilities in the past. Tittle, et a1 used cubical tanks filled with Ottawa silica sand to study the steady-state behavior of thermal and epithemal neutrons in simulated sandstone logging environments.1 In addition to the effects of fresh vs salt water, a wide variety of casing sizes and compositions were investigated. Tittman also used a large tank of Ottawa sand to measure the slowing down length of sio2;2 for making the same measurement in CaCO3, however, an assembly of solid Vermont marble blocks was used. More recently a very useful and carefully planned test pit for nuclear logs was constructed under the auspices of the API.3 This facility consists of three massive limestone sections, each 6-ft thick and more than 5 ft in diameter. The blocks are buried, one on top of the other, in a concrete- lined pit having a total depth of 24 ft. Each section is made of a different limestone to achieve porosity variation, the selection being Carthage (2 per cent), Indiana (19 per cent) and Austin (26 per cent). The blocks are traversed by a 7 7/8-in. diameter borehole and are saturated with fresh water. Great care was taken to insure total saturation of the stone. Homogeneous mixtures and materials have not been used exclusively for simulating earth formations. Tittle, et a1 4 showed that a heterogeneous lattice of water and graphite could be used to simulate Liquid furfural with very little effecton measured thermal and epi thermal neutron distributions. Subsequently, a lattice-type model for simulating limestone earth formations was designed and built by Berry This limestone model was constructed of 7/8- x 11/4- x 36411. marble rods glued together with a plastic resin to form a 36-x 42- x 42in. cube. The cube was honeycombed vertically with 7/8- x 7/8- x 3641-1. holes for porosity and had a 111/2-in. diameter hole in the center to serve as the borehole. Lower porosities were simulated by placing marble rods in the square holes. The finished cube was used with the borehole vertical and set in a tank to hold liquid for filling the holes in the rock matrix. A second type of heterogeneous borehole model has been proposed by Kukharenko, et a1.6= In this concept the model is composed entirely of solid material. The hydrogen equivalent of formation fluids is simulated by making the model of alternating layers of dense rock and celluloid. Additional hydrogen is included by drilling vertical holes in the rock sections and filling them with solid forms of hydrogeneous material. The volume of the holes is equal to the volume of the celluloid layers, so an axial and radial dispersion of the hydrogeneous material is thus obtained. One such model was composed of alternate 50 mm-thick sand layers and 2 mm-thick celluloid sheets. Approximately 200 appropriately spaced plastic tubes were used in each sand section. The plastic tubes were filled with paraffin for additional hydrogen concentration. This arrangement yielded a model with a water-equivalent hydrogen concentration of 10 per cent. Experimental data taken in this model and in a model consisting of the same material but homogeneously dispersed showed that the heterogeneous system adequately

Jan 1, 1966

By C. S. Matthews, H. M. Girner, C. D. Williams, M. J. Edwardson, H. R. Parkison

Quantitative interpretation of electric logs requires knowledge of formation temperature. In this paper, methods are developed for computing changes in formation temperature caused by circulation of mud during drilling operations. The basis of the method is the mathematical solution of the differential equation of heat conduction. The solution of this equation is presented in a series of graphs. These graphs are used to determine formation temperature disturbance at various radii for arbitrary mud circulation histories. Example comparisons with field results show reasonable agreement. It is concluded that, in general, the temperature disturbances caused by circulating mud are small beyond 10 ft from the wellbore but are quite significant near the wellbore. INTRODUCTION Quantitative interpretation of electric logs requires knowledge of the formation temperature in order to establish the resistivity of the formation water with accuracy. To determine the formation temperature, the temperature disturbances produced by circulating drilling mud must be evaluated. -The objective of this investigation was to develop a method for the numerical determination of these temperature disturbances at any distance from the wellbore as a function of time. The first step was to solve the differential equation describing the temperature behavior in the formation during and after mud circulation. This solution was used to calculate a series of curves relating temperature disturbance to shut-in time for various circulating periods. An "exact" method for computing formation temperatures with the use of these temperature-disturbance plots is described, and the results of an application of this method to a well in Montana are presented. Procedures for approximating formation temperatures are also discussed. BASIC EQUATIONS AND ASSUMPTIONS With the assumptions that (1) cylindrical symmetry exists, with the borehole as the axis, (2) heat flow is 1. The formation can be treated as though it is radially infinite and homogeneous in extent, as regards heat flow. 2. The presence of mud cake can be disregarded. 3. The effect of heat generated by bit action is negligible. 4. After mud circulation ceases, the rate of radial heat flow at the wellbore is negligible and is assumed zero for purposes of calculating temperature change with time. These assumptions are believed to be reasonable because of the agreement achieved between calculated and measured borehole temperature in a well, to be discussed later. Further, the agreement at the wellbore indicates that temperatures in the formation were also approximately correct. since it is the formation temperatures that determine wellbore temperature after ceasing circulations SOLUTION OF THE DIFFERENTIAL EQUATION AND ITS APPLICATION The solution of Eq. 2 for the temperature behavior during and after a single period of mud circulation is

By E. A. Breitenbach

A computer program has been developed to afford rapid and complete quantitative log analysis for exploration and production decisions. The computation consists of automatic selection of tops and bottoms of porous intervals from the digitized data, and then point-by-point calculations within each selected interval. Nearly all log types can be analyzed. This paper presents the calculation techniques found to be appticable to machine evaluation and gives examples of their use. INTRODUCTION As the quantitative interpretation of well logs entails repetitive use of charts and equations, it is natural that digital computer programs would be written. A Universal Log Interpretation Computer Program (ULICP) has been developed to afiord rapid and complete analysis for exploration and production decisions. Application of this program enables the log analyst to: (I) apply rapid quantitative log analysis for exploration and production decisions on either a field, a well or on a single zone basis; (2) apply concepts of interpretation requiring detailed numerical analysis; (3) analyze all porous intervals on each log, rather than a few selected zones, and complete the analysis in less time than previously required. Computations are reported for every digitized point in each zone; (4) use empirical techniques applicable to a given area as an integral part of the computation; (5) experiment with the empirical coefficients and exponents in the interpretation equations to find the best possible solution; and (6) automatically plot both the original and computed data to scales analogous to the field prints. The primary utility of a log interpretation program stems from its ability to do an overwhelming amount of work with very little man-power. If a calculation procedure or thought process can be formalized to the extent that step-by-step logic can be written, a computer program can be developed that follows this logic. The major question is one of economics. A feasibility study of the costs for such a program indicates that digital processing is economical primarily as a means of increasing the productivity of the log analyst. In effect, the probability of missing productive intervals in any well, because of a lack of time to do detailed calculations, is greatly reduced. The nucleus of ULICP is programmed to compute large sections of a suite of logs by selecting zones automatically and then performing all pertinent computations on a data point by data point basis within each zone.'.' An entire suite of logs can be processed in this manner with very little manual intervention. Sufficient programming logic is available so that each log analyst can request computations pertinent to his area. These requests are made by simple additions or deletions of information on the input header cards. Hence, the log analyst is always in complete control of the computation process. The evaluation of a suite of logs requires pre-editing, digitization, computation, presentation of results and interpretation. Work by the log analyst is reqilired only in pre-editing and interpretation. Thus, he is allowed more time for comprehensive interpretation, rather than calculation. For continuity, the discussion of ULICP is organized sequentially: pre-editing, digitization, computation and presentation of results. DISCUSSION PRE-EDITING The extent of pre-editing prior to computation is dependent on the format of the original data. For analog prints, it requires inspection, correlation and editing of the logs, plus the entering of required data an special forms. For digitized data such as magnetic tape field recordings, only the special forms are necessary. The process for analog prints will be given here. The inspection and correlation process involves the selection of sections for digitization and the correlation of the traces to a common depth. To decrease the cost of digitization, traces that exhibit large variations in shale formations can be redrawn to a non-zero baseline. The next step is to enter all pertinent data on the input header cards. The header cards presently used for ULICP are presented as Figs. 1 through 6. Cards 1, 2, 3 and 4 (Fig. 1) are defined as the Main Header Cards. They describe the particular well for output identification and give basic information. Cards 5 through 15 (Figs. 2 through 6) are defined as the Block Header Cards. As such, they define the log types and the interpretation parameters for the block of data immediately following. Cards 1 through 10 are required for every computer run. Cards 11 through 14 pertain only to nuclear logs. Card 11 is required if any nuclear log is supplied, and cards 12 through 14 are required only when gamma-ray, density or neutron log data, respectively, are supplied. Card 15 (Fig. 6) must be supplied when cross-plot calculations are required. Either two- or three-component

Jan 1, 1967

Jan 1, 1969

By S. E. Szasz, E. J. Moore, B. F. Whitney

An accurate value of formation water resistivity R, is essential in calculating formation porosity and fluid saturation from electrical well logs. In the cases where R, has not been measured directly, it must be obtained from other data, e.g., the SP curve. This paper deals with another approach: how to calculate R, from the chemical analysis of the formation water. INTRODUCTION It is known that the resistivity of aqueous solutions of pure salts depends on their concentration and on the temperature; the concentrations are given in MPL (mg of solute per liter of solution), or sometimes in ppm (mg of solute per kg of solution): MPL = ppm X specific gravity. Values for different pure salts are available in the literature, but not for solutions of mixtures which are of practical interest. The major component of the dissolved material in almost all formation waters being sodium chloride, it is customary to express the resistivity of formation waters in terms of equivalent sodium chloride concentration, i.e., the concentration of a solution of pure NaCl which has the same resistivity at a given temperature as that of the formation water under consideration. Thus, the problem of calculating R, from the chemical anaylsis can also be stated as how to convert the other constituents of the solute into equivalent NaCl concentration. Salts dissolved in water are at least partly dissociated into ions, and do not conserve their identity. If known amounts of several salts are dissolved in water, the solution does not necessarily contain the same salts in the original proportion, but perhaps some other combination of the ions, along with free ions in solution. This is why the chemical analysis of formation waters is often given in terms of ions, as if all dissolved salts were completely dissociated. Our problem then boils down to how to convert the concentrations of the various ions to equivalent concentrations of Na' and C1-. Dunlap and Hawthorne' have proposed to convert the concentration of all other ions to equivalent Na' and C1-concentrations by means of constant multipliers; e.g., 0.95 for Ca"; 2.0 for Mg"; 0.27 for HCO 3-; 0.5 for SO, -, etc. Their factors were based on measurements made at 68F on 26 formation water samples from the Texas Gulf Coast, ranging in concentration from 1,500 to 75,000 ppm. The Dunlap method is widely used in electric log interpretation, and is often extrapolated beyond its original concentration range. A comparison of R, values calculated by this method and values actually measured on formation water samples has shown large discrepancies, especially at higher concentrations. Therefore, two new methods were developed at Sinclair Oil Corp.'s Tulsa Research Center to calculate equivalent sodium chloride concentration from the chemical analysis of formation water samples. FUNDAMENTAL CONSIDERATIONS The resistivity of a solution, or its reciprocal the conductivity, at a given temperature is determined by the charge, concentration and mobility of the ions actually present. Monovalent ions such as Na' or C1- always carry the same charge. Compounds of polyvalent ions, however. may show incomplete dissociation, e.g., NaSO; + Na' instead of SO,-- + 2Na'. This happens especially in more concentrated solutions. Only very dilute solutions are completely dissociated, as assumed in the chemical analysis report. At higher concentration, the degree of dissociation depends not only on the nature and concentration of the particular salt under consideration but also on the nature and concentrations of the other solutes. Mobility of the ions depends on the viscosity of the solution. It also depends on the degree of hydration of the ions, which in turn is a function of the nature and the charge of the ions and also of the amount of free water available per ion, i.e., the total ionic concentration. The net effect is that the conductivity increases slower than proportional to the concentration, even if a solution contains only one salt such as NaC1, and is different for different salts (Fig. 1). Conductivity can even decline with a further increase in concentration, e.g., if additional salt is little dissociated but ties up some of the free water and/or causes an increase in viscosity. In solutions containing more than one salt, the contribution of one salt to the total conductivity depends not only on the fractional concentration of this same salt, but also on the concentration of all other solutes. A perfect method would give the conductivity or resistivity of a solution as a function of the concentrations of all solutes present. This is so complicated as to be impractical, and a simpler method must be found which is of acceptable accuracy. The Dunlap method, on the other hand. is too simple because it askmes that at any concentration the contrih-

Jan 1, 1967

By G. H. F. Gardner

The velocity and attenuation of elastic waves in sandstones were measured as a function of both pressure and fluid saturation. A large change occurs in these quantities if water is added and the rock is not compressed, but the change is small if the rock is subjected to a large overburden pressure. Measurements were made by vibrating cylindrical samples in both the extensional and torsional modes at frequencies up to 30,000 cycles/sec. Formulas were derived which enable the attenuation of dilatational waves in dry rocks to be deduced from the data. Similar experimental methods were used to investigate the properties of unconsolidated sands. Velocities were found to vary with the 1/4 power of the overburden pressure and attenuations to decrease with the 1/6 power. The effects of grain size, amplitude and fluid saturation were studied. Formulas by which the effects produced by a jacket around the sample may be calculated were derived. The practical application of these results to formation valuation is discussed. INTRODUCTION The attenuation of elastic waves in the earth has been of interest to the seismologist and geophysicist for many years, but only recently to the petroleum engineer. Engineering interest has been brought about by the success of velocity logging devices, for it is possible by modification of these instruments to measure the attenuation of sound waves in addition to their velocity and, hence, deduce the mobility of formation fluids as well as the porosities of the rocks which contain them. The main problem is to decide whether field measurements can be made with sufficient accuracy to be of practical use. This problem can only be solved after we know the magnitude of the attenuations which are typical of the earth at various depths. The logarithmic decrement of a fluid-saturated rock is the sum of a "sloshing" decrement and a "jostling" decrement, the former caused by the mobility of the fluid contained within the rock and the latter by the granular framework of the rock. Sloshing decrements can be calculated' using Biot's theory, but the jostling losses are less well understood. The present paper reports an experimental investigation of jostling losses in consolidated and uncon- solidated sands, particularly with respect to the effect of overburden pressure and fluid saturation. Born' showed that the decrement of a sandstone may increase dramatically when only a few per cent by weight of distilled water is added, and that the additional loss is proportional to the frequency of vibration. His measurements were made with no compressive stress on the framework of the rock. M. Gondouin3 investigated similar phenomena for fluid-saturated plasters but also did not compress the samples. In the present paper it is shown that compression of the framework reduces this effect, so that at depth the jostling decrement of a sandstone may be expected to be almost independent of fluid saturation and frequency. Decrements for many sedimentary rocks have been given by Volarovich,4 but all for the state of zero overburden pressure. Anomalously low velocities have been logged in shallow unconsolidated gas sands. Results of the present investigation confirm that these velocities are not caused by correspondingly high attenuations, because the jostling decrement in a packing of sand grains is small and much less than in a consolidated sandstone at the same depth. Velocities in sands have been measured by Tsareva5 and by Hardin6 as a function of pressure, but the corresponding decrements do not appear to have been measured previously. The widely used "resonant bar method" of measuring velocities and decrements was employed. Comments on variations of this technique have recently been published by McSkimmin.7 The main novelty of the present technique was the application of pressure to the samples. It was found possible to do this by placing the apparatus inside a pressure vessel, provided the conditions leading to large additional losses were avoided. These conditions are discussed below. EXPERIMENTAL TECHNIQUE Cylindrical samples were caused to vibrate in both the extensional and torsional mode of vibration and the amplitude of vibration was measured as a function of frequency in the neighborhood of a resonant frequency. The resonant frequency, fr, is related to the corresponding elastic modulus by the formulas where E and N are Young's modulus and the modulus of rigidity, p is the density of the sample, and A the wavelength of the vibration.

Jan 1, 1965