Search Documents

By R Saha, D Barry, C Mawdesley, O Bertoli

Kriged estimates are not typically calculated on blocks substantially smaller than the sample spacing (eg the Standard Mining Unit or SMU). Moreover, the sample spacing for quality variables is often much greater than the SMU size. For short to medium term planning decisions, the shorter term variability in the quality variables is a key economic consideration in defining, setting and meeting coal product specifications.Global uncertainties can be calculated using Drill Hole Spacing Analyses (DHSA). However, knowing the precision of a variable for a one or two year production volume is usually not adequate for variables with low spatial continuity that may fluctuate significantly over daily, weekly, SMU, or shipment volumes.Is there a feasible, reasonable and valid geostatistical method to model the uncertainty of coal quality variables on blocks substantially smaller than the drill spacing (say on the SMU)?The validity of using co-kriging variances, calculated on SMUs to characterise local uncertainty, was tested using a suite of conditional simulations generated for several quality variables for BHP’s Blackwater operation. The conditional simulation results were used as a benchmark against which the results derived from the kriging variances could be compared.The agreement between the kriging and simulation methods for thickness is excellent, with much closer correlations than for the quality variables. Typically, the kriging variances give very reasonable approximations to the simulation-based uncertainties, although there could be some systematic discrepancies which may become much more critical if a classification is to be based on a thresholding of the SMUs relative uncertainties.CITATION:Mawdesley, C, Barry, D, Bertoli, O and Saha, R, 2017. Local uncertainty benchmarking – a coal case study, in Proceedings Tenth International Mining Geology Conference 2017, pp 43–50 (The Australasian Institute of Mining and Metallurgy: Melbourne).

Sep 20, 2017

By Stephen D. Schwochow

Aggregate producers nationwide well recognize the numerous problems often faced in securing local approval to open and operate new gravel pits and rock quarries. Similar problems are faced in Colorado, but this region also presents certain political, social, and geological complexities that have become particularly troublesome in the last 10 years. The gravity of the local permitting dilemma can be related directly to the worsening supply-demand imbalance for aggregates in the state's most populated area, the Front Range urban corridor, particularly the Denver metropolitan area. The drive to develop crushed-rock reserves to supplement dwindling alluvial supplies and to meet increasing regional growth demands was somewhat predictable but perhaps more a "last resort" alternative forced onto producers by a hesitant and often apathetic government.

Jan 1, 1982

By I Douglas, L Allen

Localised indicator kriging was developed following a similar approach used by Marat Abzalov (2006) in his paper titled æLocalised uniform conditioning (LUC): A new approach for direct modelling of small blocksÆ. Abzalov introduced the idea of mapping panel distributions from uniform conditioning onto selective mining unit-sized blocks that were created using ordinary kriging. This concept of remapping panel distributions was applied to multiple indicator kriging histograms.At Newmont a program was developed that made the remapping of MIK histograms to SMU sized blocks extremely efficient and was subsequently tested on several gold deposits. The results have been encouraging as the localised indicator kriging (LIK) method has eliminated many issues that have plagued both ordinary kriging and multiple indicator kriging. Use of localised indicator kriging methodology limits the smearing of high-grade and the over-smoothing of grade that is common and difficult to control in OK but retains the ability for localising grades. LIK also simplifies change-of-support corrections, eliminates un-mineable slivers of low or high-grade when dealing with small tonnages, and allows for the mapping of linear and non-linear estimates into the same block.

Sep 26, 2011

By W Assibey-Bonsu

The paper provides a brief review of the multivariate uniform conditioning and the localised multivariate uniform conditioning techniques and presents a case study based on a porphyry copper gold deposit in Peru. Two other important scenarios have been investigated and will be particularly emphasised: The recovered selective mining unit (SMU)s depend on a combination of copper and gold content, which can be expressed by means of ænet smelter returnsÆ (NSR). Two approaches are compared. Firstly, the application of a composite cut-off on NSR on the joint distribution of each grade and the application of uniform conditioning to three variables is investigated, ie for the first approach NSR is used as the main variable. The second approach considers copper as the main variable. In the case of multiple domains the procedure used to apply the technique to process mixed panels is clarified and recommendations are made.

Sep 26, 2011

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 V. Maazallahi, A. Rouleau, A. Saeidi

"The geometry of the terrain topography before and after the failure are the essential factors for analyzing a landslide and for the design of remedial works. This paper presents a method to determine these factors, considering that the determination of the slide surface profile is difficult and costly using traditional methods and the pre-slide topography data is often not available. The presented method to locate the 3-dimensional geometry of the slide surface is based upon determining the slide surface in 2- dimensional vertical sections. This method utilizes post-slide topography of the area to determine the pre-slide topography and uses equations to estimate the pre-slide situation using a statistical process applied to points located on the sliding mass. This method is tested on the Maskun landslide in Iran. The results indicate the applicability of the proposed method to obtain the key factors for the back analysis of a landslide in a cost effective manner.INTRODUCTIONProbably the most reliable method of determining the strength of soil and rock mass is to back analyze a failed or failing slope (Caia et al. 2007; Wyllie & Mah, 2004). This method is used in geotechnics for estimating geotechnical parameters of the rock and soil mass (Sakurai et al., 2003; Frayssines and Hantz, 2009; Cepuritis et al. 2010), but it requires that the failure mode is well established and that there is complete and precise information available on the slide surface and the sliding mass (Sjöberg, 1996). Measurement of the slope profile before failure provides the basic information required for a back analysis (Hock & Bray, 1981). The post-slide topography is another key factor. The back analysis process for determining the actual strength parameters of soil or a decomposed rock mass after a landslide requires knowledge of the sliding mass geometry, both before and after the sliding event. Many studies are reported in the literature about the application of back analysis for determining geotechnical and geomechanical parameters of rock and soil mass (Sakurai et al., 2003; Cepuritis et al. 2010; Frayssines and Hantz, 2009, Hussain and Stark, 2011), but they generally present little information about the geometry of the slide surface, and the weight of rock and soil mass that slid, although these factors are the most important parameters for back analysis process. Liang and Xue-song (2012) determined the slip surface having the minimum factor of safety (called critical slip surface in deterministic analysis of slopes) and the slip surface with the minimum reliability index (called critical reliability slip surface), but they have not considered the actual profile and geometry of the slide surface. Rocscience (2013) proposes the back analysis option in the Slide software for slope stability study, but it requires an imposed slide surface as input; alternatively, one can estimate by a finite element analysis the surface with the minimum factor of safety. However, this estimated surface is very likely different from the actual sliding surface, considering the uncertainties in input data of the finite element method."

Jan 1, 2015

By W. H. Pelton

Polymetallic sulphide (PMS) deposits have recently been discovered associated with hydrothermal vents on the sea floor off Vancouver Island, B.C. Visual estimates indicate reserves in excess of 1 million tons. Reliable reserve estimates have not yet been made, however, due to the extreme difficulty of drilling cores in 2 km of water. Few geophysical techniques can aid in the exploration of these deposits. There is little or no difference in radioactivity, magnetic susceptibility or electrical resistivity between the deposits and the surrounding sea floor sediments. This precludes radiometric, magnetic, electromagnetic and resistivity methods. There is a density contrast, but the logistics of gravity surveys at such depths are formidable. The physical property which shows the most contrast with the sea floor sediments, and which is perhaps one of the easiest to measure (with modern equipment) is polarizability. The induced polarization (IP) method is one of the major electrical geophysical techniques for exploring for PMS deposits on land.

Jan 1, 1985

Fleet management systems have been adopted by many open pit mines throughout the world and the resulting improvement in truck cycle times has been one of the most significant contributors to mine efficiency. Ongoing enhancements to fleet management systems such as interfaces for vital signs monitoring have given incremental efficiency gains. Each incremental improvement is another step towards a grand vision of a totally automated mine. A relatively new fleet management system is being integrated with two other technologies to provide a highly reliable platform for a networked information system which monitors mine machinery and production. The paper reviews the integration project, the features of the new mine site information network platform and the contribution to productivity and safety. The new fleet management system has a predictive truck dispatch algorithm. It also has a distributed database that replicates rapidly so that each machine in the fleet has an identical database, thereby making each machine aware of other machinesÆ locations on the site. High precision systems for navigating and monitoring machinery such as drills, dozers, draglines and shovels provide individual machine operators with cockpit style real-time information to increase their operational efficiency. The integration of the new fleet management systems capabilities with the high precision systems that provide production information, machine vital signs data and operator information means that the mine operations and maintenance managers have real-time access to the whole mine siteÆs production data. Such a mine monitoring and reporting system is dependent on reliable wireless communications and reliable positions. Positions in fleet management and high precision systems for open pit mines are usually provided by global positioning system (GPS). Other global navigation satellite systems (GNSS) such as Glonass, and in the future Galileo, will provide additional satellite cover; however, some parts of open pit mines do not have enough sky visibility for reliable GNSS positioning and so a new centimetre level wireless positioning system called Locata is being integrated with the fleet and high precision systems to augment and backup the GNSS positions. The final integrated system will provide the highest levels of reliable position-based information to all stakeholders in a mineÆs operation, ensuring gains in productivity and safety.

Jan 1, 2008

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

By M. R. J. Wyllie

Eight laboratory-prepared aqueous base drilling muds representing common mud types, and 15 aqueous base drilling muds sampled in the field, have been used in an experimental investigation of the relationship between the streaming potential developed across a mud filter cake and the differential pressure causing filtrate flow. It has been established that the relationship is of the form, ES = kPy, where Es is the streaming potential, P the differential pressure and k and y are constants for any particular mud system at constant temperature. Beyond a certain minimum thickness this relationship is independent of filter cake thickness. Experiment has shown that both k and y are temperature dependent parameters and that in certain mud systems there is a critical value of P above which the simple relationship Es = kPY ceases to hold. The constant k tends to be a function of the mud resistivity at 77ºF, but the exponent y has not been found to be related to any commonly measured mud property. Certain field experiments are described which serve to demonstrate that the total SP (self potential) developed opposite a formation may be quantitatively accounted for as the algebraic sum of an electrochemical potential and an electrokinetic potential. A chart to permit the simple use of mean solution activities in the interpretation of electrochemical SP data is given and its use described. Various field methods for computing the electrokinetic SP component of the total SP are considered and rejected as being unduly complicated and inaccurate. It is shown that the employment of low resistivity drilling muds is essential if optimum accuracy is desired in the use of the S.P. curve for the purpose of computing connate water salinities by the electrochemical SP method. The use of such low resistivity drilling muds is accordingly recommended. INTRODUCTION In a previous paper' a quantitative analysis of the electrochemical compovent of the SP (self potential) curve was made and it was shown that the electrochemical* emf, E, could be expressed as E = ZtRT/F In a,/a,, where t is the transference number of the chloride ion, R is the gas constant, T the absolute temperature of the formation water and mud, a, the mean activity of the connate water (assumed to be sodium chloride) and a, the mean activity of the mud (as sodium chloride). The essential accuracy of this expression was confirmed by an investigation 2 of the SP kicks opposite sands of known water salinity. The logs examined had been routine runs in oil wells in a number of Mid-Continent states. Further confirmation of the original analysis of the electrochemical SP is implicit in the paper by Tixier3 wherein a relationship between the SP and the ratio of connate water and mud resistivities, based only on field data, is given. Allowing for the fact that mud and connate water resistivities were used instead of activities, and that electrokinetic effects were ignored, the constant reported by Tixier is essentially 2tRT/F. The error involved in the universal substitution of resistivities for activities will receive further reference below. Possible reasons for the fact that the electrochemical properties of shales do not appear to be significantly dependent upon either the nature or the quantity of cation exchange materials, such as clays, they contain has been considered at some length by Wyllie and Patnode.4 In the same paper, the probable effects of shale thickness and compaction on electrochemical efficiency, and hence the inherent difficulty of obtaining reliable laboratory-scale data, were also pointed out. However, it is

Jan 1, 1951

By M. R. J. Wyllie

Eight laboratory-prepared aqueous base drilling muds representing common mud types, and 15 aqueous base drilling muds sampled in the field, have been used in an experimental investigation of the relationship between the streaming potential developed across a mud filter cake and the differential pressure causing filtrate flow. It has been established that the relationship is of the form, ES = kPy, where Es is the streaming potential, P the differential pressure and k and y are constants for any particular mud system at constant temperature. Beyond a certain minimum thickness this relationship is independent of filter cake thickness. Experiment has shown that both k and y are temperature dependent parameters and that in certain mud systems there is a critical value of P above which the simple relationship Es = kPY ceases to hold. The constant k tends to be a function of the mud resistivity at 77ºF, but the exponent y has not been found to be related to any commonly measured mud property. Certain field experiments are described which serve to demonstrate that the total SP (self potential) developed opposite a formation may be quantitatively accounted for as the algebraic sum of an electrochemical potential and an electrokinetic potential. A chart to permit the simple use of mean solution activities in the interpretation of electrochemical SP data is given and its use described. Various field methods for computing the electrokinetic SP component of the total SP are considered and rejected as being unduly complicated and inaccurate. It is shown that the employment of low resistivity drilling muds is essential if optimum accuracy is desired in the use of the S.P. curve for the purpose of computing connate water salinities by the electrochemical SP method. The use of such low resistivity drilling muds is accordingly recommended. INTRODUCTION In a previous paper' a quantitative analysis of the electrochemical compovent of the SP (self potential) curve was made and it was shown that the electrochemical* emf, E, could be expressed as E = ZtRT/F In a,/a,, where t is the transference number of the chloride ion, R is the gas constant, T the absolute temperature of the formation water and mud, a, the mean activity of the connate water (assumed to be sodium chloride) and a, the mean activity of the mud (as sodium chloride). The essential accuracy of this expression was confirmed by an investigation 2 of the SP kicks opposite sands of known water salinity. The logs examined had been routine runs in oil wells in a number of Mid-Continent states. Further confirmation of the original analysis of the electrochemical SP is implicit in the paper by Tixier3 wherein a relationship between the SP and the ratio of connate water and mud resistivities, based only on field data, is given. Allowing for the fact that mud and connate water resistivities were used instead of activities, and that electrokinetic effects were ignored, the constant reported by Tixier is essentially 2tRT/F. The error involved in the universal substitution of resistivities for activities will receive further reference below. Possible reasons for the fact that the electrochemical properties of shales do not appear to be significantly dependent upon either the nature or the quantity of cation exchange materials, such as clays, they contain has been considered at some length by Wyllie and Patnode.4 In the same paper, the probable effects of shale thickness and compaction on electrochemical efficiency, and hence the inherent difficulty of obtaining reliable laboratory-scale data, were also pointed out. However, it is

Jan 1, 1951

By Milton Williams

A comparison is made between the amount of interstitial water found by analysis in cores and that estimated from the resistivity curve of the electric log for the corresponding strata of producing reservoirs. It is shown that while the determined and the estimated water saturations are in good agreement in some cases, this agreement may be due to a fortuitous compensation of errors. Estimates of water saturations made from the electric 1og by use of existing relations should, therefore, be employed with caution. Also presented are data on the variation of salinity of interstitial waters in reservoir sands, and on the effect of variation of salinity and porosity on the validity of estimates of water content made from the electric log. INTRODUCTION Throughout the history of electric logging there have been recurrent cycles of optimism and pessimism as to the value of the electric log in the quantitative estimation of reservoir fluid content. The very early days of electric logging were an optimistic period, in which there was widespread faith in the ability of the resistivity curve to reflect fluid content, and in the ability of the potential curve to indicate textural characteristics. In time the accumulation of exceptions to this idealistic behavior brought about a general disillusionment, and caused the more pessimistic of those dealing with electric logs to believe that the sole importance of the log lay in stratigraphic correlation. Within the last several years optimism as to the quantitative use of the electric log has been revived, primarily as the result of the stimulating work of C. E. Archie' on the estimation of interstitial water in sands from resistivity data. Even though the general trend of thought on the quantitative value of the electric log reflects optimism, it is inevitable that there are some widely divergent opinions of the validity of quantitative estimates based on the electric log. A valid method of estimation of fluid content by means of the electric log and a general knowledge of the textural character of the reservoir rock would be of signal importance, not only in that it would save time and money in securing information on reservoirs now being drilled, but also because it would permit the evaluation of reserves in many fields drilled without adequate core data subsequent to the advent of electric logging. Because of this importance, Humble's Production Research Division, in cooperation with various other company divisions, initiated several years ago series of tests designed to evaluate the reliability of estimates of fluid content made from the electric log. The investigation was limited at the outset to sand reservoirs, since it was felt that the inclusion of limestone reservoirs would entail needless complications. In devising these tests it was immediately apparent that the chief difficulty in evaluating the worth of any method of estimation lay in establishing the true value of fluid content. The results of direct analyses of cores cut with ordinary mud were not sufficient, because of the probability of contamination of the fluids by filtrate from the drilling mud; the results of restored state tests on cores were deemed unsuitable because these results can hardly be said to furnish a primary standard. The method adopted was the direct analysis of cores cut with an oil-base drilling fluid. Such cores, when properly handled, are believed to afford the most precise determination of water content available at the present, and, in addition, lend themselves readily to the determination of various other properties which are pertinent to the estimation of fluid content from the electric log. Since it was desired to compare the fluid contents estimated from the electric log run in ordinary mud with cores cut in the same hole, the use of oil-base fluid necessitated the changing over, after all cores were cut, to a conventional mud. This procedure presented several advantages, one of which was that the use of oil-base fluid permitted

Jan 1, 1950

By E. S. Mardock, Robert E. Bush

The object of this paper is to report recent developments in the quantitative interpretation of radioactivity logs. The use of reference lines is described in the application of the new zero radioactivity reference line for the determination of neutron derived porosities is discussed. Illustrated examples are shown of logs from the reef limestone of the North Snyder Field and from the Lansing-Kansas City limestone of central Kansas. A method of using an interrelation factor calculated from the neutron and gamma ray curves, called the neutron Productivity number, for the estimation of productivity from sandstone formations is discussed. When used in a relative manner the neutron productivity number has aided in Predicting whether a given zone will produce gas or liquid or whether it is too shaly to produce. Examples of logs from wells in the Freites District, Venezuela, and in East Texas are shown to illustrate the application of the neutron productivitv number. INTRODUCTION The petroleum industry was introduced to some of the preliminary results of investigations into the quantitative interpretation of radioactivity logs in October, 1949." It was stated at that time that the methods presented were entirely relative because they involved the use of reference points obtained from the radioactivity log, i.e., shale reference lines, dense zone reference lines, or a combination of both shale and dense zone reference lines. These relative methods offered many advantages to the industry in interpreting radioactivity logs in terms of porosity with particular application to limestone and dolomite reservoirs. However, these relative methods did not offer the wide range of applicability desired because the esta1)lishment of reliable reference points from shale in some areas were difficult. It was recognized that a more reliable mean; of establishing reference lines was vitally needed. Such a method is now available through the introduction of an instrumeutal zero device. It is the purpose of this paper to discuss the application of this new zero reference to the quantitative interpretation of radioactivity logs. Previously all applications of quantitative interpretations of radioactivity logs have been made in limestone reservoirs with little or no emphasis placed on sandstone reservoirs. A brief resumé of investigations into the quantitative interpretation of radioactivity logs in sandstone reservoirs is presented here. This discussion deals primarily with the interrelation of the gamma ray and neutron curves and the response of the latter to the density of the fluids filling the pore space in relatively homogeneous sandstones of uniform porosity. The authors believe that the methods presented herein are only applicable locally to specific fields and reservoirs, and that they probably do not offer a wide range of applicability. However, tile authors believe that the wider use of the new instrumental zero references in sand areas will reveal more information about the specific response of the gamma ray and neutron curves to changes in lithology and saturation. Two applications of quantitative interpretation using the new instrumental zero references are presented. Examples are given which show the comparison of logs and porosity in the limestones and dolomites of West Texas and Kansas. Logs of wells in 'Venezuela and East Texas are compared with the drill stem test data, production data, etc.. to illurtrate the interrelation of the gamma ray and neutron curve: and the rise of methods involving these interrelations to determine certain saturation variables. THE QUANTITATIVE APPLICATION OF RADIOACTIVITY LOGS TO LIMESTONE AND DOLOMITE RESERVOIRS Tile theoretical basis for the gamma ray and neutron curves has been discussed elsewhere.9 he method of recording the logs and apparatus involved has also been thoroughly covered and will not be discussed here.'."".'" It has been shown that the neutron curve responds exponentially to liquid filled porosity over the 1 1/2 to 35 per cent

Jan 1, 1951

By Robert E. Bush, E. S. Mardock

The object of this paper is to report recent developments in the quantitative interpretation of radioactivity logs. The use of reference lines is described in the application of the new zero radioactivity reference line for the determination of neutron derived porosities is discussed. Illustrated examples are shown of logs from the reef limestone of the North Snyder Field and from the Lansing-Kansas City limestone of central Kansas. A method of using an interrelation factor calculated from the neutron and gamma ray curves, called the neutron Productivity number, for the estimation of productivity from sandstone formations is discussed. When used in a relative manner the neutron productivity number has aided in Predicting whether a given zone will produce gas or liquid or whether it is too shaly to produce. Examples of logs from wells in the Freites District, Venezuela, and in East Texas are shown to illustrate the application of the neutron productivitv number. INTRODUCTION The petroleum industry was introduced to some of the preliminary results of investigations into the quantitative interpretation of radioactivity logs in October, 1949." It was stated at that time that the methods presented were entirely relative because they involved the use of reference points obtained from the radioactivity log, i.e., shale reference lines, dense zone reference lines, or a combination of both shale and dense zone reference lines. These relative methods offered many advantages to the industry in interpreting radioactivity logs in terms of porosity with particular application to limestone and dolomite reservoirs. However, these relative methods did not offer the wide range of applicability desired because the esta1)lishment of reliable reference points from shale in some areas were difficult. It was recognized that a more reliable mean; of establishing reference lines was vitally needed. Such a method is now available through the introduction of an instrumeutal zero device. It is the purpose of this paper to discuss the application of this new zero reference to the quantitative interpretation of radioactivity logs. Previously all applications of quantitative interpretations of radioactivity logs have been made in limestone reservoirs with little or no emphasis placed on sandstone reservoirs. A brief resumé of investigations into the quantitative interpretation of radioactivity logs in sandstone reservoirs is presented here. This discussion deals primarily with the interrelation of the gamma ray and neutron curves and the response of the latter to the density of the fluids filling the pore space in relatively homogeneous sandstones of uniform porosity. The authors believe that the methods presented herein are only applicable locally to specific fields and reservoirs, and that they probably do not offer a wide range of applicability. However, tile authors believe that the wider use of the new instrumental zero references in sand areas will reveal more information about the specific response of the gamma ray and neutron curves to changes in lithology and saturation. Two applications of quantitative interpretation using the new instrumental zero references are presented. Examples are given which show the comparison of logs and porosity in the limestones and dolomites of West Texas and Kansas. Logs of wells in 'Venezuela and East Texas are compared with the drill stem test data, production data, etc.. to illurtrate the interrelation of the gamma ray and neutron curve: and the rise of methods involving these interrelations to determine certain saturation variables. THE QUANTITATIVE APPLICATION OF RADIOACTIVITY LOGS TO LIMESTONE AND DOLOMITE RESERVOIRS Tile theoretical basis for the gamma ray and neutron curves has been discussed elsewhere.9 he method of recording the logs and apparatus involved has also been thoroughly covered and will not be discussed here.'."".'" It has been shown that the neutron curve responds exponentially to liquid filled porosity over the 1 1/2 to 35 per cent

Jan 1, 1951

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 P. Majani, F. P. Kokesh, M. Grosmangin

Determination of the quality of cementation of casing in oil wells in the past has involved inflow and circulation tests to insure that the producing zones were adequately sealed off from the adjacent zones. Existing logging methods, such as temperature and radioactivity surveys, may detect the presence of cement behind the casing. However, the qualities of the cement (i.e., its hardness and particularly its bond to the casing) are not indicated. The new logging method described in this paper operates on the principle that the attenuation of a sonic pulse transmitted by a casing is greatly increased when that casing is bonded to an outer annulus of hard material (such as set cement) which has an appreciably smaller sonic-wave velocity than that of the casing. The down-hole tool contains a source of recurrent sound pulses which are detected by a receiver spaced a few feet from the source The amplitude of the detected casing-borne pulse is measured, and the resulting signal is transmitted to the surface where it is recorded vs depth. Since amplitude is a function of attenuation, the log is readily interpreted. Laboratory studies have shown straightforward relationship between attenuation and such variables as source-detector spacing and per cent of circumference bonded. It is shown that cement not set or not bonded to the casing has compara- tively little attenuating effect. Field examples show not only the cement top, but also the variation in cementation quality below the top. Further, the increase of bonding with time and after squeeze cementation is depicted. The detection of poor cement jobs is confirmed by production tests and by formation-test results. It is anticipated that the method will have wide application in evaluating cementation quality prior to formation testing in completions and recompletions. The analysis it affords may aid in further improving cementation techniques. INTRODUCTION The main purpose of oilwell cementing is to isolate a production zone from other undesirable zones. To investigate whether this purpose has been accomplished, several logging methods have been used such as temperature logs, radioactive tracer logs, etc. While these logs all respond to the presence of cement behind the casing, they do not indicate the degree of bonding of the cement to the casing. Early in the application of sonic logging, it was noticed that considerable attenuation of sound signals takes place in cemented pipe and is often made evident on the standard Sonic log by cycle skipping.' The development of a circuit capable of continuously recording the amplitude of the casing-borne sound signal has made possible an extensive series of laboratory and field tests, which gave the following results. The amplitude of a sound signal after it has traveled in a firmly cemented pipe is only a small fraction of that recorded by the same device in free pipe. This provides a wide spectrum of energy levels; for given local conditions, empirical values of the amplitude can be correlated with the quality of cementation. Interpretations made in this manner generally have been confirmed by production tests, circulation tests and squeeze cementation. The purposes of this paper are to give a general description of the new logging method and to present some laboratory and field results. CYCLE SKIPPING Early attempts to study the quality of the cement behind the casing were performed with a standard Sonic log, which measures the transit time At, and were based on the well known phenomenon of cycle skipping.' Cycle skipping normally is interpreted as being a manifestation of weak signals at the receivers. The log of Fig. 1 was run with the recording instrumentation adjusted to enhance cycle skipping. A mirror-image presentation was used for better visual interpretation. The transit time of sound in steel is about 58 microseconds/ft (corresponding to a velocity of 17,000 ft/ sec), and the portions of the log where this value is recorded are interpreted as zones with no cement bond. Where cycle skipping produces a higher value of At, weak signals (or a high rate of attenuation) are indicated, and a good cement bond may be present. The occurrence of cycle skipping, however, depends too much on instrument adjustment to give a uni-

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 J. P. Timmons, J. H. Moran, G. K. Miller, M. A. Coufleau

A prototype equipment has been designed and built for the digital recording of well logs on magnetic tape at the same time that the regular film recording is made. The format of the digital tape produced is such that it can be used directly at the input of the ZBM 704, 7090 or other models of ZBM computers which accept digital magnetic tape. This apparatus has been used for the experimental field recording of dipmeter tape logs which were subsequently computed by means of an ZBM 704 or 7090. In this paper the equipment and the digital tape are described briefly, and their application to the computer-interpretation of dipmeter data is discussed. A principal element in the interpretation of the dipmeter log is the correlation of the three microresirtivity dipmeter curves to determine the depth displacements between them. Several correlation methods for computer use are considered, with particular attention to their sensitivity to error and their consumption of computer time. The tape data were used to compute information content of the dipmeter microresistivity curves in terms of their frequency spectra. The results show that the sampling rate used in recording the digital information is quite adequate and illustrate a use of the digital tape in evaluating the characteristics of new tools. Some examples of field results are shown. It can be foreseen that, when digital tape recording becomes available for general field use, a whole new realm of possibilities will be opened up for the processing of other well logs through computations, which hitherto were not feasible because they were too laborious and time-con.sunzing. INTRODUCTION The last few years have seen a revolution in the design and production of data-processing equipment. Stored-pro-gram digital computers have progressed from a research curiosity to the basis of a major industry. There are now hundreds of such machines in daily use in the United States. With the acceptance of a technique that was, in fact, already clearly described by John von Neumann in 1945, the last decade has seen great strides in the development'of components, reliability, programming systems and, most spectacularly, in the sheer number of machines built and in use. In 1957 there were enough digital computers available to the oil industry to justify the suggestion that it would be worthwhile to investigate the possibility of using these machines in processing well log data.' The first result of this investigation was the appearance of what may be referred to as the input-output bottleneck. Well logs are customarily recorded on film. To get these data into a machine required then (and still does): a time-consuming semi-automatic reading of the film; conversion of the log data to digital form; and recording these digital data in some medium acceptable for computer input, such as cards, magnetic tape, or punched paper tape. However, the recording, reading, and re-recording could only result in deterioration of the data. Therefore, it was concluded that the fist step should be the development of a new, more direct recording technique supplemental to the film recording, which would provide easy access to the digital computer. There are many solutions to the problem of recording log data in an easily recoverable form. After careful consideration it was decided to adopt the boldest solution which, it was felt, was also the most elegant. It was decided to record well logs directly, in the field, on magnetic tape in such a way that this tape could be used without further modification as an input to the IBM 704 or 7090 computer. To realize practical field recording of magnetic tape logs, it became necessary to develop in a rather small package, an analog-to-digital converter, a tape recorder, and the necessary multiplexing and control circuits to allow the simultaneous recording of a multiplicity of logging signals. The magnetic tape recording was to be made simultaneously with the conventional logging operation in such a way as not to interfere with it. Along with the development of hardware, it was necessary to begin development of interpretation techniques and machine programs that would exploit the power of the digital computer. Here, again, there is a long list of possible applications. After much consideration it was decided to concentrate on the interpretation of the dipmeter log as a first application. It is the object of this paper to describe in some detail the developments sketched in the last three paragraphs.

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 G. H. Pardue, L. H. Gollwitzer, J. H. Moran, R. L. Morris

Experimental and theoretical studies of casing-borne signals have provided data on the attenuation efect, compressive strength and shear modulus for different types of cement. These data resulted in the establishment of a definite relationship between the cement compressive strength and the Cement Bond log reading. The relationship is influenced by the casing dimensions, but is independent of the type of cement and the curing conditions. The theoretical analysis provides a mathematical statement of the phenomenon and shows that the attenuation rate is basically a function of the shear modulus of the cement and the thickness of the casing. The lab tests provided a connection between the shear modulus and the compressive strength, and permitted a comparison of the experimental and theoretical work. The results of the experimental and theoretical study have provided means of quantitative interpretation of the Cement Bond log. This interpretation is made more convenient with the use of a nomograph. INTRODUCTION Past interpretation of Cement Bond logs has been essenlially based upon a qualitative comparison between the log response in the zone of interest and the responses in the "fully cemented" and the "uncemented" intervals. Most of the variables influencing the magnitude of the signals were recognized in an earlier paper'; however, insufficient laboratory data were available at that time for quantitative evaluation of these variables. Early attempts to use the numerical amplitudes of the signals without full knowledge of their significance resulted in erroneous comparisons between various types of cement. An extensive program was undertaken to evaluate, both experimentally and theoretically, the attenuation phenomenon and to study the effects of the more common variables. A major objective of the experimental program was to obtain data on the relationship between cement-bond-tool response and physical property standards used by the industry so that practical interpretation methods could be devised. The studies included a comparison of attenuation rate with these properties as a function of time, temperature and pressure, using common sizes of casing. The main result was that a definite relationship was established between the compressive strength of the cement in inti- mate contact with the casing and the response of the cement bond tool. For a given casing thickness, regardless of the type of cement or the curing conditions used, the attenuation rate is a function of the compressive strength. A theoretical analysis showed that the attenuation rate for casing-borne signals should be a function of (among other things) the shear modulus of the cement and the thickness of the casing. The relationship between the compressive strength and shear modulus for curing times similar to those encountered in the field was determined from laboratory tests. Using this relationship, the computed attenuation rates showed reasonable agreement with the measured attenuation rates. The purpose of this paper is to report the results of the experimental and theoretical investigations and to present the first practical method for a quantitative interpretation of Cement Bond logs based on known variables. This method is applicable where the casing-borne signal is not affected by high-velocity formation arrivals. DESCRIPTION OF LOGGING TOOL The Cement Bond logging tool shown schematically in Fig. 1 has been described previously.' The tool consists of two sections: acoustic and electronic. The acoustic section contains a transmitter and a receiver, usually spaced 3-ft apart. Sonic energy emittted by the transmitter travels through the casing and is detected at the receiver. The electronic section measures the amplitude of the