Underworld-GT applied to Guangdong, a tool to explore the geothermal potential of the crust

Abstract

Geothermal energy potential is usually discussed in the context of conventional or engineered systems and at the scale of an individual reservoir. Whereas exploration for conventional reservoirs has been relatively easy, with expressions of resource found close to or even at the surface, exploration for non-conventional systems relies on temperature inherently increasing with depth and searching for favourable geological environments that maximise this increase. To utilitise the information we do have, we often assimilate available exploration data with models that capture the physics of the dominant underlying processes. Here, we discuss computational modelling approaches to exploration at a regional or crust scale, with application to geothermal reservoirs within basins or systems of basins. Target reservoirs have (at least) appropriate temperature, permeability and are at accessible depths. We discuss the software development approach that leads to effective use of the tool Underworld. We explore its role in the process of modelling, understanding computational error, importing and exporting geological knowledge as applied to the geological system underpinning the Guangdong Province, China.

Keywords

Underworld-GT / geothermal potential / computational modelling approach / Guangdong

References

[1] Balay S, Eijkhout V, Gropp W D, . Arge E, Bruaset A M, Langtan-gen H P, . Efficient Management of Parallelism in Object Oriented Numerical Software Libraries. Modern Software Tools in Scientific Computing, 1997 Boston: Birkhauser, 163-202.

[2] Batory D, Sarvela J N, Rauschmayer A. Scaling Step-Wise Refinement. Proceedings of the 25th international Conference on Software Engineering. International Conference on Software Engineering. IEEE Computer Society, Portland, May 3–10, 2003, 2003, 187-197.

[3] Capitanio F A, Faccenna C, Zlotnik S, . Subduc-tion Dynamics and the Origin of Andean Orogeny and Bolivian Orocline. Nature, 2011, 480: 83-86.

[4] Chopra P, Holgate F. A GIS Analysis of Temperature in the Australian Crust. Proceedings World Geothermal Congress, Antalya, April 24–29, 2005, 2005

[5] Danis C. Use of Groundwater Temperature Data in Geothermal Exploration: The Example of Sydney Basin, Australia. Hydrogeology Journal, 2014, 22: 87-106.

[6] Danis C, O’Neill C, Quenette S. Is It Hot enough down there? Assessing Geothermal Potential in the Sydney-Gunnedah-Bowen Basin. Australian Geothermal Energy Conference, Sydney., 2012, 14-16.

[7] Farr T G, Rosen P A, Caro E, . The Shuttle Radar Topography Mission. Rev. Geophys., 2007, 45 RG2004.

[8] Gropp W D, Lusk E. User’s Guide for mpich, a Portable Implementation of MPI, 1996

[9] Honza E, Fujioka K. Formation of Arcs and Back-Arc Basins Inferred from the Tectonic Evolution of Southeast Asia since the Late Cretaceous. Tectonophysics, 2004, 384: 23-53.

[10] Kiczales G, Lamping J, Mendhekar A, . Aksit M, Matsuoka S, . Aspect-Oriented Programming. Proceedings of the European Conference on Object-Oriented Programming, 1997 Berlin, Heidelberg, and New York: Springer-Verlag, 220-242.

[11] Lafaurie B, Nardone C, Scardovelli R, . Modelling, Merging and Fragmentation in Multiphase Flows with SURFER. Journal of Computing Physics, 1994, 113: 134-147.

[12] J. Geophys. Res., 2010, 115 B7

[13] Liu T Y, Wu Z C, Zhan Y L, . Wavelet Multi-Scale Decompostition of Magnetic Anomaly and Its Application in Searching for Deep-Buried Minerals in Crisis Mines: A Case Study from Daye Iron Mines. Earth Science—Journal of China University of Geosciences, 2007, 32: 31-41.

[14] Ma X H. Outline of Chinese Lithospheric Dynamics, 1987 Beijing: Geological Publishing House

[15] Mallet J L. Turner K. GOcad: A Computer Aided Design Program for Geological Applications. Three-Dimensional Modelling with Geoscientific Information Systems: Nato ASI Series C, 1992 Dordrecht: Kluwer Acad. Publ., 123-141.

[16] Meixner A J, Holgate F L. The Cooper Basin Region 3D MapVersion 1: A Search For Hot Buried Granites, 2009

[17] Monaghan J J. Smoothed Particle Hydrodynamics. Annual Review of Astronomy and Astrophysics, 1992, 30: 534-574.

[18] Moresi L, Dufour F, Mulhaus H B. A Langrangian Integration Point Finite Element Method for Large Deformation Modeling of Viscoelastic Geomaterials. Journal of Computing Physics, 2003, 184: 476-497.

[19] Moresi L, Quenette S, Lemiale V, . Computational Approaches to Studying Non-Linear Dynamics of the Crust and Mantle. Physics of the Earth and Planetary Interiors, 2007, 163: 69-82.

[20] Musson A, Harrison B, Gordon K, . Budd A R, Gurgenchi H, . Thermal Thinking: Optimal Targeting for Australian Geothermal Explorers. Proceedings of the 2009 Australian Geothermal Energy Conference, 2009

[21] Pavlis N K, Holmes S A, Kenyon S C, . An Earth Gravitational Model to Degree 2160: EGM 2008, 2008 Vienna: General Assembly of the European Geosciences Union

[22] Pohl K, Böckle G, van der Linden F. Software Product Line Engineering—Foundations, Principles, and Techniques, 2005

[23] Quenette S. Software Evolution and Computational Performance—The Semantics of Software as an Inverse Problem, 2014 Melbourne: Monash University

[24] Quenette S, Hodkinson L. StgDomain—Scalable Parallel Domain Software Components for Particle-in-Cell Finite Element Methods. Concurrency and Computation: Practice and Experience, 2009, 22(12): 1593-1603.

[25] Quenette S, Moresi L. Models Based Experimentation: Numerical Modelling of 3D Basin Scale Architecture Heat & Fluid Flow. AGU Fall Meeting Abstracts, December 10, 2010, San Francisco, 2010

[26] Quenette S, Moresi L, Sunter P D, . Explaining StGermain: An Aspect Orientated Rnvironment for Building Extensible Computational Mechanics Modelling Software. HIPS Workshop, Parallel and Distributed Processing Symposium, Proceedings of the 19th IEEE International. 210, 2007

[27] Somerville M, Wyborn D, Chopra P N, . Hot Dry Rock Feasibility Study, 1994 Canberra: ERDC

[28] Spector A, Grant F S. Statistical Models for Interpreting Aeromagnetic Data. Geophysics, 1970, 35: 293-302.

[29] Stegman D R, Moresi L, Turnbull R, . Next Generation Visualization Framework for High-Performance Computational Geodynamics. Visual Geosciences, 2008, 13(1): 71-84.

[30] Sulsky D, Zhou S J, Schreyer H L. Application of a Particle-in-Cell Method to Solid Mechanics. Comput. Phys. Commun., 1995, 87: 236-252.

[31] USGS Shuttle Radar Topography Mission. Global Land Cover Facility, 2004 College Park, Maryland: University of Maryland

[32] Wan T. The Tectonics of China: Data, Maps and Evolution, 2011 Beijing: Higher Education Press, 173.

[33] Wand M, Kiczales G, Dutchyn C. “A Semantics for Advice and Dynamic Join Points in Aspect-Oriented Programming,” ACM Trans. Programming Languages and Systems, 2004, 26(5): 890-910.

[34] Yang W C, Shi Z Q, Hou Z Z, . Discrete Wavelet Transform for Multiple Decompostion of Gravity Anomalies. Chinese Journal of Geophysics, 2001, 44: 534-541.

[35] Zhang B T, Wu J Q, Ling H F, . Estimate of Influence of U-Th-K Radiogenic Heat on Cooling Process of Granitic Melt and Its Geological Implications. Science in China Series D: Earth Sciences, 2007, 37: 155-159.

[36] Zhang Z, Xu T, Zhao B, . Systematic Variations in Seismic Velocity and Reflection in the Crust of Cathaysia: New Constraints on Intraplate Orogeny in the South China Continent. Gondwana Research, 2013, 24(3): 902-917.

[37] Zhou L, Xie J, Shen W, . The Structure of the Crust and Uppermost Mantle beneath South China from Ambient Noise and Earthquake Tomography. Geophysical Journal International, 2012, 189(3): 1565-1583.

[38] Zhou X, Sun T, Shen W, . Petrogenesis of Mesozoic Granitoids and Volcanic Rocks in South China: A Response to Tectonic Evolution. Episodes, 2006, 29 1 26.