Search Documents

By K George

The Dugald River lead-zinc-silver deposit is located approximately 60 km north-west of Cloncurry and 87 km north-east of Mount Isa. The development project is expected to turn operational in 2018. The resource is estimated to contain 61.1 Mt at 12.6 per cent zinc, 1.9 per cent lead and 34 g/t silver, with a projected mine life of 25 years. The evolution from exploration project to development project has presented multiple challenges for the technical team. Understanding and modelling the structural complexity of the deposit has been one of these challenges.Historically fault zones were logged as a lithology, with structural measurements recorded in a separate database table. When underground drill programs and development commenced in late 2012, it was recognised the deposit was more structurally complex than had been modelled. The extensive number of fault zone lithologies in core logging, made it difficult to identify fault trends and model faults in 3D software packages.In 2014, a structural project was conducted to produce a fault model for the deposit. The project highlighted that by logging the faults as a lithology, a shortage of fault characteristics was being recorded. Additionally, structural information could not be filtered to identify different fault sets due to the lack of information available to filter on. One of the major outcomes of the study was to change the methodology for the structural logging of drill core.Faults are now logged in a separate structural log. The structural log collects two types of information: interval data and point data. The interval data represents the damage zone of the fault in the surrounding rock and the point data the plane of movement for the fault. The faults are further characterised by their intensity on a scale of one to five, the type of fault structure (fault breccia, fault vein etc), the minerals present on the fault surface and structural measurements when available.The new data set allows fault data to be filtered based on fault characteristics, making implicit and explicit modelling of the faults more timely and accurate. The data set can be used to produce fault wireframes as single surfaces from the point data or hanging wall and footwall surfaces from the interval data. The logging method and how its applied in fault modelling form the focus of this paper.CITATION:George, K, 2017. A new approach to logging and modelling faults at the Dugald River Pb-Zn-Ag deposit, North West Queensland, in Proceedings Tenth International Mining Geology Conference 2017, pp 315–320 (The Australasian Institute of Mining and Metallurgy: Melbourne).

Sep 20, 2017

During the first five years of operation, the Lihir gold mine has steadily increased plant throughput rate to over 4 Mtpa, well above the original design capacity of 2.8 Mtpa. Understanding how future variations in ore type impact on plant performance has been an integral part of LihirÆs development program. This paper discusses the key drivers for significant improvements to the plant, including a second oxygen plant, a pilot flotation plant, an autoclave heat recovery circuit, a pebble crushing circuit and additional autoclave feed pumping capacity. The result is a process plant with the flexibility to accommodate variations in ore type without significant loss of throughput.

Jan 1, 2003

By R Talebi

A research project to better understand the rock mass behaviour and the seismic response has been conducted during the design and construction of the 9215 level at the E-Block sector of the Northern Star Kanowna Belle Mine. The research program was undertaken at a tunnel located 1100 m below the surface with the results applicable to feasibility studies on underground mining projects in highstress conditions, where failure can occur very soon after tunnel construction. Semi-circular excavation shapes coupled with the use of de-stress blasting techniques were designed to monitor the excavation in conditions of very high stress. The ground support was designed considering high energy dissipation ground support schemes (HED) as described by Arcaro et al (2021). The results indicate a decreased rock mass demand due to the change of excavation shape. Also, localised reduction in rock mass stiffness was achieved due to the implementation of de-stress blasting. Finally, increased ground support capacity was gained by the implementation of the high energy dissipation ground support schemes. The Kanowna Belle Mine has significant ore reserves at depths exceeding 1000 m below the surface. The results of this research will be used to optimise the extraction of the remaining E Stoping Block at depths exceeding 1350 m.

Mar 29, 2023

By T N. Little

This paper is concerned with the influence of engineering geology on design of Wall Control Blasts. In the Introduction, the Blast Management Framework and the Blast Design Requirements classification system (previously developed by the author) is drawn upon to set the context. Blasting objective for all wall control blasts involves both profile and damage control. Four classes of approaches for blast design are outlined and briefly discussed. Five distinct wall control blasting techniques are summarised. The general relationship between geology and blast design is explored in the paper. Influence of geology on pre-split and post-split design and blasting direction is highlighted. Wall damage modelling from the closest row of blastholes is also introduced using the Analytical Damage Model (ADM). The relationship of rock properties and explosive selection and in particular the Mach wave generation is introduced to assist blast designer in this regard. Finally a Conceptual Model that relates vibration, damage and rock mass condition has been presented along with some applications that have been undertaken by the developer. Eight conclusions have been formulated based on extensive work in this area and the current study.

Jun 22, 2022

By H S. Pescini

A material aim when feeding ore blends into a processing facility is to reduce the variability of the properties of the feed material and to maximise the output, or throughput. At the Mining Geology conference of 2019, we presented a paper on optimal blending of coking coal given geological variability. The paper introduced integer-programming techniques to both maximise value and minimise feed or product variability. We discussed an application where the geological uncertainty is characterised using a normal distribution and exact solutions are readily obtainable using second order conic constraints in the mathematical model. Geological uncertainty is not always normally distributed, many properties will be characterised as log normal or some other skewed distribution. For example, copper can be some form of lognormal and Fe in Iron ore has a left skewed distribution. In this paper, we discuss some distribution free techniques to optimise the blends of plant feed or output given a quantified geological uncertainty. We discuss the sample average approximation method and how to set-up and apply the method to an ore blending problem. Evaluating a case study demonstrates that the explicit use of non-normally distributed geological uncertainty can better optimise feed into an ore processing facility and potentially add significant value. The approach is equally applicable to all mining and processing operations that blend material into a processing facility or directly into a saleable product.

Mar 22, 2022

By Y Lanwin

The Australian mining industry, a cornerstone of the nation’s economy, is on a journey towards greater equity and inclusion. Th paper introduces a comprehensive exploration of how universal design principles, stemming from architectural roots, have evolved through automation and digital technology into the more encompassing frameworks of inclusive design and equitable design. The paper proposes a practical framework for those involved in work, facilities and equipment design in the mining industry to identify and assess for equity and inclusion. Universal Design emerged in the 1960s with a primary focus on architecture, aiming to create spaces accessible to people of all abilities. Its foundational principles sought to ensure equitable use, flexibility, simplicity, perceptibility, tolerance for error, low physical effort and optimal spatial considerations. Over time, these principles evolved and expanded beyond physical spaces, finding applications in various fields, including technology, education and transportation. In response to advancements in technology and a growing understanding of the multifaceted nature of diversity, Universal Design evolved into Inclusive Design. Inclusive Design goes further than mere accessibility, considering a broader spectrum of users with varying cognitive, sensory, cultural and socio-economic backgrounds. It emphasises user centred design, customisation options, diverse user research and continuous iteration. Equitable Design emerged as a framework that addresses not only the diversity of needs but also systemic disparities. It emphasises fairness, justice and the elimination of structural inequities. Equitable Design encompasses accessibility and affordability, representation and diversity, ethical considerations and community engagement. This paper will explore the applications of these design principles through a series of case studies including village and on-site facilities, remote operations and office based work design and technology design. By applying these frameworks in a mining context, the paper proposes that mining companies and suppliers can improve diversity, equity and inclusion leading to safer, more engaged and productive workplaces.

Apr 16, 2024

By C Artica

Two quantitative methods for Mineral Resource classification have been applied to a copper skarn deposit beneath a large open pit that is mining a world-class porphyry complex. A drill hole spacing study (DHSS) using the single block kriging method (SBK) and a conditional simulation method (CS) were applied to the skarn deposit and the results compared. Inputs to both processes included the proposed mining rate and stope designs to estimate the likelihood of achieving grade and metal production targets. Mineral Resource classification of the deposit considers the relative precision of tons, grade and metal predictions over quarterly and annual production periods. These relative precision estimates are calculated at 90 per cent confidence interval (CI). The SBK method represented the production period as a single large block. Kriging estimates were performed in the large block for various synthetic drill hole grids, and the estimation errors derived. The error was assumed to be normally distributed due to the large block size and quantity of samples, allowing CIs to be deduced. The CS method followed a conventional sequential Gaussian simulation workflow where grades were simulated with a 5 ft × 5 ft × 5 ft grid node spacing inside mineralisation domains. Each realisation was then averaged up to monthly production volumes represented by a combination of designed stopes, and the coefficient of variation (CV) for the averaged grade was calculated. CI on the relative precision of grade predictions were then deduced from the CVs. The two methods produced similar results; despite being calculated independently. CS results were used to report uncertainty in the proposed life-of-mine plan and highlight areas for additional drilling. SBK results were used as the basis for Mineral Resource classification as this method has been previously applied to the skarn deposit.

May 24, 2023

By P Dowd, A Adeli, C Xu, X Emery

Geometallurgical modelling is increasingly being incorporated into mineral resource modelling and estimation as a means of increasing efficiency, decreasing operating costs and reducing risk in mining operations. Traditional geostatistical workflows for geometallurgical modelling face several difficulties, in particular, related to the non-additivity of process response variables, to highly heterotopic data sets and/or to preferential sampling designs, all of which may explain why most attempts to estimate geometallurgical response variables from assay data and mineralogy are still based on statistics and machine learning. The aim of the work presented here is to demonstrate the applicability of geostatistics to estimate copper recovery from copper sulfide ores at the Prominent Hill Iron Oxide Copper-Gold (IOCG) deposit in Australia. As recovery is non-additive, the preferred workflow is the joint estimation of the masses of metal in the feed and in the concentrate, both of which are additive variables. While the mass of metal in the feed can be estimated from abundant online assay analyses, the mass of metal in the concentrate can only be estimated from a very limited number of laboratory-scale batch flotation tests. Traditional approaches to estimating the masses of metal in the feed and in the concentrate, such as simple and ordinary co-kriging, are described and discussed. A modified version of co-kriging is introduced, which incorporates linear relationships between the mean values of the input variables, in addition to a spatial correlation model. Co-kriging with related means is shown to outperform traditional simple and ordinary co-kriging in terms of precision, accuracy and consistency of the estimates, and offers a practical and efficient tool to deal with highly heterotopic and preferential sampling designs.

May 24, 2023

By A Batnasan, K Haga, A Shibayama, K Higuchi, H Takeuchi

High silica (SiO2), alumina (Al2O3), and phosphorus (P) contents in iron ore adversely affect the quality of iron products. The removal of those impurities from iron ore is an essential point in further iron ore processing. In this paper, the removal of P, SiO2 and Al2O3 impurities from iron ore by alkaline drying followed by reverse cationic flotation is discussed. The reverse cationic flotation experiments are performed using the iron ore samples without and with alkaline drying pretreatment. All other flotation conditions except pulp pH (4–12) and particle size distributions (-45 to -250 μm) are held constant. A typical collector and frother, such as dodecylamine acetate (DAA) and methylisobutylcarbinol (MIBC), are used in the flotation experiments. Alkaline treatment is conducted before flotation using sodium hydroxide at fixed conditions. The impact of alkaline drying on the performance of iron ore flotation is evaluated based on the removal rates of P and SiO2 + Al2O3, iron grade, and iron recovery compared with the iron ore flotation without alkaline drying (direct flotation). X-ray diffraction (XRD) analysis indicates hematite and goethite are the main minerals in the iron ore, which assays 62.8 mass per cent Fe, 0.232 mass per cent P, 7.1 mass per cent SiO2 + Al2O3 (total) determined by X-ray fluorescence analysis (XRF). Results indicate that a combined approach is more effective for removing P and SiO2 + Al2O3 impurities from iron ore than the direct processing technique because the maximum removal rates of P and SiO2 + Al2O3 are 48.2 per cent and 66.5 per cent from the alkaline dried iron ore when 8.4 per cent and 13.1 per cent from the non-alkaline dried iron ore.

Sep 18, 2023

By Y Li, D J. Williams, E Topal

A mine waste rock dump does not have any intrinsic value, yet is an inevitable structure in surface mining operations. Since the purpose of mining is to maximise profit, possible expenditure reduction in hauling waste rock to the dump, which can represent up to 50 per cent of the total operating cost, is highly desirable. Apart from cost, potential acid mine drainage (AMD), due to seepage from surface waste rock dumps containing reactive waste rock, could threaten the ecological balance, drinking water supplies and irrigation systems downstream from a mine site. In order to minimise these potential risks, selective placement of benign and reactive waste rock, and the encapsulation of reactive waste rock by benign waste rock, need to be carried out during the construction of the rock dump.In this paper, a new methodology, based on mixed integer programming (MIP) is developed to manage waste rock dumping. The model aims to minimise waste rock haulage costs by locating the optimum destination(s) for each waste rock block, while at the same time controlling the placement of reactive waste rock to minimise the potential for AMD. Tests on small-scale virtual data demonstrated that the model is capable of dealing with selective placement of different waste rock types, encapsulating reactive, potentially acid forming (PAF) waste rock, and stockpiling benign waste rock for encapsulation and rehabilitation purposes. The solution time for the MIP model is far less than the traditional manual approach, and the resulting schedule is up to 46 per cent more cost effective.CITATION:Li, Y, Topal, E and Williams, D J, 2012. Mathematical approach for better mine waste rock dumping management, in Proceedings Life-of-Mine 2012 , pp 207-214 (The Australasian Institute of Mining and Metallurgy: Melbourne).

Jul 10, 2012

By T Horberry

This paper describes work conducted as part of an overall research project that aims to identify current and emerging trends in the development of automation and remote mining technologies, and to provide an operator-centred perspective to the design and deployment of such systems. More specifically, this paper outlines a research and development process being constructed to examine the issues of the human skills and capabilities required to operate or maintain automated/remote controlled mobile mining equipment. A human factors approach presented here, where the focus is upon designing equipment interfaces to fit operators, rather than purely designing and deploying the equipment based on a technology-centred approach. However, the importance of training to help develop the required skills and capabilities is also emphasised.The work built on a full review of the literature about the effects of automation in mining and elsewhere, extensive interviews with different stakeholders in mining automation and recognised best practice regarding designing for the human element in automated systems. The human element research and development process proposed is in four stages, ranging from understanding the tasks involved in both manual and automated equipment operation through to user trials to ensure that the automated systems are acceptable to operators and are fit for purpose.It is outside the scope of the current paper to present the full data collection described by this four stage process, but an overview of what will be collected is presented. As such, this paper focuses on the broad methods being developed, rather than the specific results to be obtained. The ultimate objective is to apply a broad ‘technology assessment’ approach to the design, deployment and evaluation of automated mining technologies, in particular by considering the user-centred design implications of these technologies.

Nov 22, 2011

By K Eglinton

The extractive industries response to community engagement due to social and political pressures in the 1980s and 1990s was the starting point of development leading to the modern function of corporate social responsibility.Challenged by communities who have greater access to information, demand greater transparency and responsibility of corporate organisations and who possess greater awareness of the communities they participate in, the mining sector faces another step change in how it approaches and integrates corporate social responsibility, transparency and engagement.Mining is seeking reserves of lower grades, in more sensitive environments and in countries with minimal transparency and governance. Not only are the locations changing but so to the stakeholders to engage with and issues to contend with. An increasingly savvy opposition, including not just the ‘green hippies’ but the Gen-Ys interested in the wider sustainability debate and Baby Boomers moving to life style environments, are now questioning mining’s behaviour and encroachment into ‘their space’. Issues including effect on lifestyle amenity, property value, energy efficiency, social provision; concepts such as ‘doing more with less’, participatory engagement, increasing transparency and balancing economic, environment, cultural and social aspects are becoming the norm.The traditional model of supporting local communities, communicating plans and providing employment, community investment and indigenous support programs is not enough. Our industry faces challenges on multiple fronts. It must now be able to address these from a personal/community perspective if it is to remain relevant within communities and play a role in addressing the global sustainability issues.This paper draws on corporate social responsibility trends, integrated functions and implementation approaches through the application of activities within Newmont operations across Asia Pacific.The approach presented works on the philosophy that communities do not necessarily have to like what we do, but should be able to identify and recognise the way in which we involve them in decision making, communicate our intentions and respond to issues.

Nov 22, 2011

By R J. Thompson

In truck-based hauling systems, the mine haul road network is a critical and vital component of the production process. As such, underperformance of a haul road will impact immediately on mine productivity and costs. Operations safety, productivity and equipment longevity are all dependent on well-designed, constructed and maintained haul roads. With the advent of autonomous haul trucks, the haul road itself becomes all the more critical to the success of these type of operations; not only in relation to mine operators requirements for safer and more efficient and predictable haulage systems, but also in response to autonomous truck manufacturers’ requirements for a more predictable and controlled operating environment. This paper presents a brief summary of the state-of-the-art in surface mine road design and then proceeds to examine the design and technological challenges associated providing a safe, predictable and affordable road for autonomous mine haul trucks. The paper serves as a basis for evaluating the contributions that enhancements in road design and management can deliver to autonomous mining operations, an initiative that also has scope for adoption within the operating life of many mines, haul roads and haul trucks.

Nov 22, 2011

By J J. Grobler

The George Fisher Underground Mine at Mount Isa applies the long hole open stoping methodology. During the reconciliation process of transverse stopes in the George Fisher North mine it was identified that excessive underbreak was regularly experienced, with the major contribution coming from the toe area inside the stope.While over break and the resultant dilution was well investigated and well understood, limited work was done on underbreak improvement. It was decided that a re-engineering process would be followed in an attempt to improve results. The prevailing stoping practices were scrutinised and the possible causes were listed.Some of the identified causes were design, drilling accuracy, equipment condition, mucking practices and direction of firing. It was identified that in most instances the underbreak was initiated in the slot area of the stope during the first firings. Once this happened, it was impossible to get good recoveries on the rest of the stope. The slot area was carefully scrutinised and the blasting practices were changed individually.During the following two years several different strategies were trialled to increase the toe recoveries. The results indicated an overall improvement of six per cent reduction in underbreak on average, by ensuring drilling accuracy and changing the direction of firing. The economic significance of these results outweighed the time and additional work load expended to establish the cause of the problem. This paper outlines the continuous improvement process undertaken to identify and rectify the major contributing factors to the problem. It also presents a set of recommendations to be followed in order to achieve positive results.

Mar 21, 2011

By R Winterberg

Concrete is a well established and reliable construction material for the permanent lining of tunnels, but shows some weak performance characteristics, which can be overcome by the addition of specialised fibre reinforcement. This paper addresses the particular performance design requirements of segmental linings, discussing the improvements of fi bres on concrete properties. Focus is made on the beneficial effects of fibre reinforcement on the cracking behaviour and the toughness performance, providing a safe and durable tunnel lining.Not only the mechanical performance of fibre reinforced concrete (FRC) for segmental linings, but also the durability of the linings has gained growing attention among the designers and consultants. Especially tunnel linings with exposure to media with corrosion potential, such as wastewater sewer tunnels or desalination plant outfall tunnels, have high requirements on their durability. The design of these structures must provide a comprehensive study of durability aspects. Fibre reinforcement can substantially improve the durability of concrete.Research has been carried out on the cracking behaviour of concrete. The beneficial effects of steel fibres on the cracking behaviour reduce the permeability by reducing the flow paths within the concrete. Crack bridging steel fibres create necking effects in bending cracks, leading to earlier closure of the cracks. Further, steel fibres create particular crack formations, which lead to a substantial reduction of the fl ow rates through separation cracks. The results also show that considerable cost savings can be achieved using FRC, as the thickness of the linings can be reduced.The durability design issues will be addressed and underlined by a case history on the Hobson Bay sewer tunnel in Auckland, New Zealand.

Mar 8, 2011

By H Benzer

Grinding clinker to finished cement consumes more energy than any other comminution process. World production of cement in 2009 was 2840 million tonnes and the energy required for comminution was 40-60 kWh/tonne depending on the clinker characteristics and the fineness of the finished cement. Machines used include single and multi compartment ball mills, HPGRs, vertical roller mills, centrifugal air separators, V separators, screens and crushers. The circuit most commonly used consists of a closed circuit HPGR followed by a single compartment ball mill and an air separator, and in operation the variables are adjusted to obtain maximum throughput with no loss of quality of the separator fines.A common objective in cement plants is to reduce energy consumption without affecting the particle size in the finished cement and simulation is used as one procedure to achieve this. Models are now available for all the machines commonly used in clinker circuits and circuit optimisation studies using simulation have been carried out in many plants leading to decrease in energy consumption per tonne of cement. In this paper the models used will be discussed briefly and several case studies will be reviewed which will highlight the effect of change in ball size in a two compartment mill, the role of the HPGR in an HPGR – ball mill – air separator circuit, and the use of a high intensity vibrating mill for producing finished cement with high additive ratio.

Aug 8, 2011

By D J Maconochie

The use of rock bolts as permanent support has been prohibited by client specifications in several recent large tunnel projects in Australia. The justification for this appears to lie in concerns with proper and verifiable installation procedures rather than concerns with the intrinsic durability of the corrosion protection systems. However, in civil tunnels, preventing the use of permanent rock bolts for long-term support results in the use of cast in place concrete lining which in some instances would otherwise be unnecessary.

Mar 8, 2011

By P Moser, R Ivanova, F Ouchterlony

This paper describes model-scale bench blasts on mortar blocks with distorted drill hole patterns, done in 2013–2014 at the Montanuniversitaet Leoben. The main aim of the project is to see if the fragmentation from and the crack development behind blasts with and without drill hole deviations differ significantly.The dimensions of the test blocks were 660 × 280 × 210 mm (L × W × H). They were mounted inside a yoke that allows the blast waves to escape (see Figure 1). Three rows in each block were shot row by row with a nominal pattern of B × S = 70 × 95 mm. The same nominal delay 73 µs (ie 1 ms per m burden) between holes was used throughout. Blocks with holes of stochastic collar position errors (variations in burden or in burden and spacing but constant row volume) and systematically shifted collar positions (blasts with staggered pattern) were shot. The results were compared to reference blocks without collar deviations.Apart from the sieving analysis, the remaining bench faces were measured with a stereo-photogrammetric method and horizontal roughness profiles for each row were constructed. Sawing slabs from block remains after the third row and use of dye penetrant allowed the study of internal blast damage.The analysis showed that none of the changes in borehole collaring had a significant effect on the sieving curve (median). The sieving results confirm earlier findings though that the fragmentation gets finer with the number of rows shot, implying that blast damage from earlier rows has an influence on the blasting results. The tests demonstrated that the nequiv (equivalent n-value) has increased ten per cent by shooting staggered pattern and decreased by seven per cent by pattern with burden and spacing errors. For the rest of the patterns, the nequiv was roughly the same.The surface roughness investigation showed that a smooth surface was achieved by the staggered pattern. The other deviations patterns did not show a clear trend.The internal blast damage, which was measured in the remaining part of blocks, showed that out of eleven, four crack families seemed to be most affected by blasting. The comparison of the reference and distorted patterns showed that less cracks were apparent for the staggered pattern, ie less damage was introduced.CITATION:Ivanova, R, Ouchterlony, F and Moser, P, 2015. Influence of distorted blasthole patterns on fragmentation as well as roughness of, and blast damage behind, remaining bench face in model-scale blasting, in Proceedings 11th International Symposium on Rock Fragmentation by Blasting, pp 693–706 (The Australasian Institute of Mining and Metallurgy: Melbourne).

Aug 24, 2015

Research cruises dedicated to seafloor hydrothermal activity along the Kermadec arc first began with the Sonne-135 expedition in 1998. This was followed by the worldÆs first systematic survey for venting along an arc during the 1999 NZAPLUME I expedition. Since then, cruises dedicated to mapping plumes along the Kermadec arc include NZAPLUME II (2002), NZAPLUME III (2004) and ROVARK (2007). In addition, the manned submersible cruises SWEEP Vents (2004) and NZASRoF (2005) utilised the Shinaki 6500 and Pisces V, respectively, to sample in detail select hydrothermalsystems for mineralised samples, rocks and vent fluids. The ROVARK (2007) cruise was the first to deploy an autonomous underwater vehicle (ABE) along the arc which surveyed Brothers volcano at a far greater resolution than had been done previously by ships on the sea surface. While these expeditions have largely been led by GNS Science, together with its collaborators from National Oceanic and Atmospheric Administration, Pacific Marine Environmental Laboratory (NOAA/PMEL), significant contributions have been made by other groups from various nations, exemplifying the global interest in seafloor hydrothermal activity associated with intraoceanic arcs. This is an ASTRACT only. A full-length paper was not prepared for this presentation.

Jan 1, 2008

This paper presents a review of computer techniques for the automatic contouring of geological and geophysical data with output as interactive display on graph plotters or graphics terminals and passive display on automatic drafting tables. The advantages of the method of fitting a triangulated surface to randomly distributed data using an advancing front algorithm are described, as well as the technique for determining interpolated values for the production of contour line segments and a contour smoothing algorithm. This contouring method is employed by the Geological Survey of Queensland and is implemented as an interactive system. An alternative method ideally suited to randomly distributed data involves the use of piecewise parametric cubic splines. This method which is ideal for sharp variations in magnitude and areal extent typical of geophysical anomalies is discussed. The mathematics of surface interpolation using non-parametric splines and an explicit piecewise surface fitting technique are discussed. Curve tracing techniques which are suitable for automatic contouring are also presented and the inherent advantages of interactive contouring using on-line plotters and graphics terminals for previewing of displays are described.

Jan 1, 1977