Abstract:
The geothermal system at Lihir Island, Papua New Guinea, is difficult to model as the top surface of the model evolves with time as the pit for the gold mine is excavated. A method is described for automatically including the evolving topography of the Lihir mine in numerical simulations of the geothermal system. The method makes extensive use of Python scripts to control the simulation process and in particular exploits the abilities of the PyTOUGH libraries to operate directly on TOUGH2 input and output files.