Wells for Engineered Geothermal Systems (EGS) typically occur in conditions presenting significant challenges for conventional rotary and percussive drilling technologies: granitic rocks that reduce drilling speeds and cause substantial equipment wear. Thermal spallation drilling, in which rock is fragmented by high temperature rather than mechanical means, offers a potential solution to these problems. However, much of the knowledge surrounding this drilling technique is empirical - based on laboratory experiments that may or may not represent field conditions. This paper outlines a new numerical modeling effort investigating the grain-scale processes governing thermal spallation drilling. Several factors affect spall production at the mesoscale, including grain size and size distribution, surface temperatures and material heterogeneity. To investigate the relative influence of these factors, we have conducted a series of simulations using GEODYN - a parallel Eulerian solid and fluid dynamics code. In this paper, we describe a two-dimensional model used to simulate the grain-scale processes and present preliminary results from this modeling effort.

Widespread adoption of geothermal energy will require access to deeply buried resources in granitic basement rocks at high temperatures and pressures. Exploiting these resources necessitates novel methods for drilling, stimulation, and maintenance, under operating conditions that are often difficult or impossible to reproduce in laboratory settings. Physically rigorous numerical modeling tools are vital to highlight potential risks, guide process optimization and reduce the uncertainties involved in developing new technologies for these environments. Lawrence Livermore National Laboratory has developed, and is constantly improving, several multi-physics solid/structural mechanics, fluid dynamics, chemistry, and discrete element codes. Integration of the LLNL simulation tools into a coherent simulation environment will provide a predictive capability for the thermomechanical response - in particular the spall and fracture - of basement rocks at high temperatures and pressures useful for drilling and other geothermal applications. This paper outlines a modeling effort investigating the processes involved in hydrothermal spallation drilling. These include interconnected phenomena on several length and time scales: from system-scale fluid dynamics and heat transfer of the high temperature jet to the rock face to the grain-scale thermomechanics of spallation. Three models are described to capture these different scale processes: a grain-scale model to investigate the onset of spallation; a particulate fluids model to simulate the transport of the produced spalls; and a borehole-scale model to represent the integrated system behavior.

The Engineered Geothermal System (EGS) Exploration Methodology Project is developing an exploration approach for EGS through the integration of geoscientific data. The overall project area is 2500km2 with the Calibration Area (Dixie Valley Geothermal Wellfield) being about 170km2. The Final Scientific Report (FSR) is submitted in two parts (I and II). FSR part I presents (1) an assessment of the readily available public domain data and some proprietary data provided by terra-gen power, llc, (2) a re-interpretation of these data as required, (3) an exploratory geostatistical data analysis, (4) the baseline geothermal conceptual model, and (5) the EGS favorability/trust mapping. The conceptual model presented applies to both the hydrothermal system and EGS in the Dixie Valley region. FSR Part II presents (1) 278 new gravity stations; (2) enhanced gravity-magnetic modeling; (3) 42 new ambient seismic noise survey stations; (4) an integration of the new seismic noise data with a regional seismic network; (5) a new methodology and approach to interpret this data; (5) a novel method to predict rock type and temperature based on the newly interpreted data; (6) 70 new magnetotelluric (MT) stations; (7) an integrated interpretation of the enhanced MT data set; (8) the results of a 308 station soil CO2 gas survey; (9) new conductive thermal modeling in the project area; (10) new convective modeling in the Calibration Area; (11) pseudo-convective modeling in the Calibration Area; (12) enhanced data implications and qualitative geoscience correlations at three scales (a) Regional, (b) Project, and (c) Calibration Area; (13) quantitative geostatistical exploratory data analysis; and (14) responses to nine questions posed in the proposal for this investigation. Enhanced favorability/trust maps were not generated because there was not a sufficient amount of new, fully-vetted (see below) rock type, temperature, and stress data. The enhanced seismic data did generate a new method to infer rock type and temperature (However, in the opinion of the Principal Investigator for this project, this new methodology needs to be tested and evaluated at other sites in the Basin and Range before it is used to generate the referenced maps. As in the baseline conceptual model, the enhanced findings can be applied to both the hydrothermal system and EGS in the Dixie Valley region).

Proposed drill sites for intermediate depth temperature gradient holes and/or deep resource confirmation wells. Temperature gradient contours based on shallow TG program and faults interpreted from seismic reflection survey are shown, as are two faults interpreted by seismic contractor Optim but not by Oski Energy, LLC.