Three-dimensional gravity models are a useful part of improving the geological understanding of large areas in various geological settings. Such models can assist seismic interpretation, particularly in areas of poor seismic coverage. In general, forward modelling and inversion are conducted until a single model is derived that fits well to the observed gravity field. However, the value of such a model is limited because it shows only one possible solution that depends on a fixed set of underlying assumptions. These underlying assumptions are not always clear to the interpreter and an arguably more useful approach is to prepare multiple models that test various scenarios under a range of different assumptions. The misfit between observed and calculated gravity for these various models helps to highlight flaws in the assumptions behind a particular choice of physical parameters or model geometry. Identifying these flaws helps to guide improvements in the geological understanding of the area. We present case studies for sedimentary basins off western Africa and western Australia. The flawed models have been used to rethink assumptions related to the geology, crustal structure and isostatic state associated with the basins, and also to identify areas where seismic interpretation might need to be revised. The result is a more reliable interpretation in which key uncertainties are more clearly evident.

Abstract # : 1479734 Paper # : GP43B-1142 Session : GP43B Potential-field and EM methods for geologic problems of the mid and upper crust Developments for 3D gravity and magnetic modeling in spherical coordinates Richard Lane - Geoscience Australia - rjllane@gmail.com Qing Liang - China University of Geosciences (Wuhan) - qingliang.cug@gmail.com Chao Chen - China University of Geosciences (Wuhan) - chenchao@cug.edu.cn Yaoguo Li - Colorado School of Mines - ygli@mines.edu At Geoscience Australia (GA), Australia's Commonwealth Government geoscientific agency, we perform gravity and magnetic modeling at a range of scales, from broad regional crustal studies with thousands of kilometer lateral extent and tens of kilometer vertical extent, to detailed local studies with kilometer or less lateral extent and meters to hundreds of meters vertical extent. To achieve greater integration and coherence, and to better understand the geological significance of this work, we are investing in a number of development projects; * Spherical coordinate gravity and magnetic modeling, * Modeling using High Performance Computing facilities, * Utilizing rock property data as an input into the modeling and interpretation of gravity and magnetic data, * Better management of geoscience data and models, and * Visualization of spatial data in a Virtual Globe format. In collaboration with the Colorado School of Mines (CSM) and the China University of Geosciences (CUG), we are developing a capability to model gravity and magnetic data in a spherical coordinate framework. This will provide more accurate calculations and permit us to integrate the results into a single framework that more realistically reflects the shape of the Earth. Modeling gravity and magnetic data in a spherical coordinate framework is far more compute intensive than is the case when performing the corresponding calculations in a Cartesian (rectangular) coordinate framework. To reduce the time required to perform the calculations in a spherical coordinate framework, we will be deploying the modeling software on the National Computational Infrastructure (NCI) High Performance Computing (HPC) facility at the Australian National University (ANU). This will also streamline the management of these software relative to the other main option of establishing and maintaining HPC facilities in-house. We are a participant in the Deep Exploration Technologies Cooperative Research Centre (DET CRC). In combination with this involvement, we are expanding our support for systematic management of rock property data, and developing a better understanding of how these data can be used to provide constraints for the modeling work. We are also using the opportunities afforded through the DET CRC to make progress with documentation and standardization of data storage and transfer formats so that the tasks of management, discovery and delivery of this information to users are simplified and made more efficient. To provide the foundations of integration and analysis of information in a spatial context, we are utilizing and customizing 3D visualization software using a Virtual Globe application, NASA World Wind. This will permit us to view the full range of information types at global to local scales in a realistic coordinate framework. Together, these various development activities will play an important role in the on-going effort by Geoscience Australia to add value to the potential field, rock property, and geological information that we possess. We will then be better able to understand the geology of the Australian region and use this knowledge in a range of applications, including mineral and energy exploration, natural hazard mitigation, and groundwater management.

Abstract for SGA Townsville 2009

The Paterson AEM survey was flown over the Paterson Orogen, the eastern Pilbara Craton and the on-lapping Officer and Canning Basins in NW Western Australia between September 2007 and October 2008 as part of the Commonwealth Government's Onshore Energy Security Program. The survey was designed to provide pre-competitive data for enhancing uranium and other mineral exploration. Flight lines were at a variety of spacings from 6, 2 and 1 km to 200 m targeting known deposits and other covered highly prospective rocks for a total area of 45,330 km2. The survey data has afforded new insights into the Paleozoic paleotopography of the region which is blanketed by regolith including Phanerozoic sediments including Permian glaciogene, Mesozoic and Cenozoic sediments. These insights have major implications for mineral prospectivity.

Like many of the basins along Australia's eastern seaboard, there is currently only a limited understanding of the geothermal energy potential of the New South Wales extent of the Clarence-Moreton Basin. To date, no study has examined the existing geological information available to produce an estimate of subsurface temperatures throughout the region. Forward modelling of basin structure using its expected thermal properties is the process generally used in geothermal studies to estimate temperatures at depth in the Earth's crust. The process has seen increasing use in complex three-dimensional (3D) models, including in areas of sparse data. The overall uncertainties of 3D models, including the influence of the broad assumptions required to undertake them, are generally only poorly examined by their authors and sometimes completely ignored. New methods are presented in this study which will allow estimates and uncertainties to be addressed in a quantitative and justifiable way. Specifically, this study applies Monte Carlo Analysis to constrain uncertainties through random sampling of statistically congruent populations. Particular focus has been placed on the uncertainty in assigning thermal conductivity values to complex and spatially extensive geological formations using only limited data. As a case study these new methods are then applied to the New South Wales extent of the Clarence-Moreton Basin. The geological structure of the basin has been modelled using data from existing petroleum drill holes, surface mapping and information derived from previous studies. A range of possible lithological compositions was determined for each of the major geological layers through application of compositional data analysis. In turn, a range of possible thermal conductivity values was determined for the major lithology groups using rock samples held by the NSW Department of Primary Industries (DPI). These two populations of values were then randomly sampled to establish 120 different forward models, the results of which have been interpreted to present the best estimate of expected subsurface temperatures, and their uncertainties. These results suggest that the Clarence-Moreton Basin has a moderate geothermal energy potential within an economic drilling depth. This potential however, displays significant variability between different modelling runs, which is likely due to the limited data available for the region. While further work could improve these methods, it can be seen from this study that uncertainties can provide a means by which to add confidence to results, rather than undermine it.

Exploration models for Rot Rock geothermal energy plays in Australia are based primarily on high-heat producing granites (HHPG) in combination with overlying low-conductivity sedimentary rocks providing the insulator necessary to accumulate elevated temperatures at unusually shallow (therefore accessible) depths. Unknowns in this style of geothermal play include the composition and geometry of the HHPG and thermal properties, and the thickness of the overlying sediments. A series of 3D geological models have been constructed to investigate the range of geometries and compositions that may give rise to prospective Hot Rock geothermal energy plays. A 3D geological map of the Cooper Basin region which contains known HHPG beneath thick sedimentary sequences, has been constructed from gravity inversions and constrained by geological data. The inversion models delineate regions of low density within the basement that are inferred to be granitic bodies. Thermal forward modelling was carried out by incorporating measured and estimated thermal properties to the mapped lithologies. An enhancement of the GeoModeller software is to allow the input thermal properties to be specified as distribution functions. Multiple thermal simulations using Monte-Carlo methods would be carried out from the supplied distributions. Statistical methods will be used to yield the probability estimates of the in-situ heat resource, reducing the risk of exploring for heat. The two thermal modelling techniques can be used as a predictive tool in regions where little or no temperature and geological data are available.

Broken Hill Managed Aquifer recharge Projects 3D models and Fly-through

The northern Perth Basin is an under-explored part of the southwest continental margin of Australia. Parts of this basin have proven hydrocarbon potential. The basin is extensively covered by mostly 2D seismic reflection data and marine gravity and magnetic data. The seismic data helps to resolve the structural framework of the basin, but in deepwater regions, the basement-cover contact and deeper basement structure are generally not well imaged. To help overcome this limitation, integrated 3D gravity modelling was used to investigate crustal structure in onshore and offshore parts of the basin. Such modelling also relies on knowledge of crustal thickness variations, but these variations too are poorly constrained in this area. Multiple models were constructed in which the seismic data were used to fix the geometry of sedimentary layers and the fit to observed gravity was examined for various different scenarios of Moho geometry. These scenarios included: 1) a Moho defined by Airy isostatic balance, 2) a Moho based on independently-published Australia-wide gravity inversion, and 3) attempts to remove the Moho gravity effect by subtracting a long-wavelength regional trend defined by GRACE/GOCE satellite data. The modelling results suggest that the best fit to observed gravity is achieved for a model in which the thickness of the crystalline crust remains roughly constant (i.e. deeper Moho under sediment depocentres) for all but the outermost parts of the basin. This finding has implications for understanding the evolution of the Perth Basin, but remains susceptible to uncertainties in sediment thickness.

The GA Animator Software is a tool used to create high quality fly-through animations of geoscience data for internal and external stakeholders. It is build using the NASA World Wind Java SDK, as a companion tool to the publicly available 3D Data Viewer.

Predictive 3D geological models of the subsurface can be developed using a range of available tools. Each tool is suited to slightly different problems and datasets. The method described here, using the UBC-GIF inversions algorithms, allows rapid development of models using an objective, automated procedure. It has flexibility to include as little or as much geological information as is available, making it ideal for greenfields exploration or mapping programs. The steps involved are: 1) develop a solid understanding of the expected physical properties; 2) convert geological observations into physical property constraints; 3) perform geologically-constrained inversions; 4) apply geological classifier to recovered 3D physical property models. The procedure is demonstrated for the southern Agnew-Wiluna greenstone belt in WA, a highly mineralised region with a high proportion of surface cover. The predictive 3D lithology models developed for the area are particular effective at mapping the extent of dense mafic and magnetic ultramafic rocks, and provide new insights about their distribution at depth.