Egbert Jolie 1, James Faulds 2, Inga Moeck 1.

Similar documents
STRUCTURAL ASSESSMENT AND 3D GEOLOGICAL MODELING OF THE BRADY'S GEOTHERMAL AREA, CHURCHILL COUNTY (NEVADA, USA): A PRELIMINARY REPORT

Targeting of Potential Geothermal Resources in the Great Basin from Regional Relationships between Geodetic Strain and Geological Structures

Shevenell et al. 2004, Geothermal Resources Council Transactions, v. 28: 47-52, Accomplishments at the Great Basin Center for Geothermal Energy

Nevada Bureau of Mines and Geology, University of Nevada, Reno, NV U.S. Geological Survey, Menlo Park, CA 94025

Structural Constraints of Buffalo Valley Hot Springs, North-Central Nevada

3D analysis of geothermal fluid flow favorability: Brady s, Nevada, USA

3D geophysical inversion modeling of gravity data to test the 3D geologic model of the Bradys geothermal area, Nevada, USA

Geothermal Potential Assessment in Northern Rwanda

Geothermal Systems: Geologic Origins of a Vast Energy Resource

Bog Hot Valley. (updated 2012)

NOTICE CONCERNING COPYRIGHT RESTRICTIONS

Ann Moulding and Tom Brikowski University of Texas at Dallas, Department of Geosciences

ANN ROBERTSON TAIT PRESENT POSITION. Business Development Manager / Senior Geologist EXPERTISE

A Revised Interpretation of 3D Seismic Data, Hawthorne Army Depot, Nevada: Faulted Basin Reflections or Sill Intrusions?

Aeromagnetic map of the Death Valley ground-water model area, Nevada and California

Structural deformation across the southwest Mina deflection, California-Nevada: Field studies in the Huntoon Springs area.

Geothermal Resource Characterization and Evaluation at Astor Pass, Nevada

Empirical Green s Function Analysis of the Wells, Nevada, Earthquake Source

Answers: Internal Processes and Structures (Isostasy)

Structural Style and Tectonic Evolution of the Nakhon Basin, Gulf of Thailand

The Nevada Play Fairway Project Phase II: Initial Search for New Viable Geothermal Systems in the Great Basin Region, Western USA

2. Governing Equations. 1. Introduction

Tu D Understanding the Interplay of Fractures, Stresses & Facies in Unconventional Reservoirs - Case Study from Chad Granites

James E. Faulds and Nicholas H. Hinz. Nevada Bureau of Mines and Geology, University of Nevada, Reno, NV

Simulations of Carbon Dioxide Injection, Seismic Monitoring, and Well Logging for. Enhanced Characterization of Faults in Geothermal Systems

Structural Controls of a Blind Geothermal System in the Northern Pyramid Lake Area, Northwestern Nevada

predictive mineral discovery*cooperative Research Centre A legacy for mineral exploration science Mineral Systems Q1 Architecture

Petroleum geology framework, West Coast offshore region

HDR/EGS POTENTIAL OF THE VILKYCIAI AREA, WEST LITHUANIA

DETAILED IMAGE OF FRACTURES ACTIVATED BY A FLUID INJECTION IN A PRODUCING INDONESIAN GEOTHERMAL FIELD

Seismic Reflection Data and Conceptual Models for Geothermal Development in Nevada

GEOTHERMAL WELL SITING USING GIS: A CASE STUDY OF MENENGAI GEOTHERMAL PROSPECT

Summary. Introduction

Kinematic inversion of pre-existing faults by wastewater injection-related induced seismicity: the Val d Agri oil field case study (Italy)

Marine and Land Active Source Seismic Imaging of mid Miocene to Holocene aged Faulting near Geothermal Prospects at Pyramid Lake, Nevada

Structural Controls of the Tuscarora Geothermal Field Elko County, Nevada

Slip- and dilation tendency analysis: Implications for geothermal exploration in the Upper Rhine Graben

Activity Pacific Northwest Tectonic Block Model

Thermal Modeling of the Mountain Home Geothermal Area

Regional Patterns of Geothermal Activity in the Great Basin Region, Western USA: Correlation With Strain Rates

ABSTRACT. * ,

Seismic Reflection Imaging across the Johnson Ranch, Valley County, Idaho

Characterization of Subsurface Permeability of the Olkaria East Geothermal Field

Structural investigations of Great Basin geothermal fields: Applications and implications

Fault seal analysis in Move

Building confidence of CCS using knowledge from natural analogues

Overview of the Seismic Source Characterization for the Palo Verde Nuclear Generating Station

A New Method of Evaluation of Chemical Geothermometers for Calculating Reservoir Temperatures from Thermal Springs in Nevada

Exploration of Geothermal High Enthalpy Resources using Magnetotellurics an Example from Chile

EVALUATING HEAT FLOW AS A TOOL FOR ASSESSING GEOTHERMAL RESOURCES

Today: Basic regional framework. Western U.S. setting Eastern California Shear Zone (ECSZ) 1992 Landers EQ 1999 Hector Mine EQ Fault structure

Overview of Indonesian Geothermal System

SUPPLEMENTARY INFORMATION

Summary. Introduction

Mapping the Preferential Flow Paths within a Fractured Reservoir

Conceptual model for non-volcanic geothermal resources - examples from Tohoku Japan

Great Basin Center for Geothermal Energy

Chapter 16. Mountain Building. Mountain Building. Mountains and Plate Tectonics. what s the connection?

Applicability of GEOFRAC to model a geothermal reservoir: a case study

NUMERICAL MODELING STUDY OF SIBAYAK GEOTHERMAL RESERVOIR, NORTH SUMATRA, INDONESIA

The Geologic Framework of the Fallon FORGE Site

Earth Science, (Tarbuck/Lutgens) Chapter 10: Mountain Building

The Geologic and Structural Framework of the Fallon FORGE site

FLUID STRATIGRAPHY OF THE COSO GEOTHERMAL RESERVOIR

Keywords: Play fairway, Great Basin, permeability, geothermal potential map, structural setting, geostatistics

GeothermEx, Inc. GEOTHERMAL RESERVOIR ASSESSMENT METHODOLOGY FOR THE SCIENTIFIC OBSERVATION HOLE PROGRAM, KILAUEA EAST RIFT ZONE, HAWAII TASK 1 REPORT

ABSTRACT. Introduction

Project Overview. Justyna Ellis, Ernst Huenges, Simona Regenspurg

Integration of Data in a Play Fairway Analysis of Geothermal Potential Across the State of Hawaii

Worldwide Power Density Review

Structural Geology tectonics, volcanology and geothermal activity. Kristján Saemundsson ÍSOR Iceland GeoSurvey

Geophysical exploration methods at European sites

Preface and Overview. Folded strata in the mountains of Italy (ca AD), Leonardo da Vinci

Copyright McGraw-Hill Education. All rights reserved. No reproduction or distribution without the prior written consent of McGraw-Hill Education

Determination of Calcite Scaling Potential in OW-903 and OW-914 of the Olkaria Domes field, Kenya

Using Fuzzy Logic to Identify Geothermal Resources and Quantify Exploration Risk through Play Fairway Analysis

Identification of natural fractures and in situ stress at Rantau Dedap geothermal field

Seismic Reflection and Magnetotelluric Imaging of Southwestern Dixie Valley Basin, Nevada

3D Geometry of the Xade Complex inferred from Gravity and Magnetic Data

HIGH TEMPERATURE HYDROTHERMAL ALTERATION IN ACTIVE GEOTHERMAL SYSTEMS A CASE STUDY OF OLKARIA DOMES

Hijiori HDR Reservoir Evaluation by Micro-Earthquake Observation

Composition of the earth, Geologic Time, and Plate Tectonics

Chapter 15 Structures

STRATIGRAPHIC RELATIONSHIPS IN MESOZOIC BASEMENT ROCKS AT THE DESERT PEAK EAST EGS AREA, NEVADA

Estimation of fracture porosity and permeability using radon as a tracer

GIS INTEGRATION METHIOD FOR GEOTHERMAL POWER PLANT SITEING IN SABALAN AREA, NW IRAN

Beneath our Feet: The 4 Layers of the Earty by Kelly Hashway

Lateral extrusion and tectonic escape in Ilan Plain of northeastern Taiwan

caribbean basins, tectonics and hydrocarbons university of texas institute for geophysics

SEISMIC MONITORING OF EGS STIMULATION TESTS AT THE COSO GEOTHERMAL FIELD, CALIFORNIA, USING MICROEARTHQUAKE LOCATIONS AND MOMENT TENSORS

Focal Mechanism Analysis of a Multi-lateral Completion in the Horn River Basin

Modelling the surface heat flow distribution in the area of Brandenburg (northern Germany)

Japan Engineering Consultants, Inc., Energy and Industrial Technology Development Organization,Tokyo, Japan

11.1 Rock Deformation

Seismicity in Texas in Relation to Active Class I Underground Injection Control Wells: Preliminary Observations

PROCEEDINGS, INDONESIAN PETROLEUM ASSOCIATION Thirty-Ninth Annual Convention and Exhibition, May 2015

Status of geothermal energy exploration at Buranga geothermal prospect, Western Uganda

Tectonic Seismogenic Index of Geothermal Reservoirs

Taller de Geotermica en Mexico Geothermal Energy Current Technologies

Depositional History and Petroleum Potential of Ombilin Basin, West Sumatra - Indonesia, Based on Surface Geological Data*

Transcription:

PROCEEDINGS, Thirty-Seventh Workshop on Geothermal Reservoir Engineering Stanford University, Stanford, California, January 30 - February 1, 2012 SGP-TR-194 THE DEVELOPMENT OF A 3D STRUCTURAL-GEOLOGICAL MODEL AS PART OF THE GEOTHERMAL EXPLORATION STRATEGY A CASE STUDY FROM THE BRADY S GEOTHERMAL SYSTEM, NEVADA, USA Egbert Jolie 1, James Faulds 2, Inga Moeck 1 1 GFZ German Research Centre for Geosciences, Telegrafenberg, 14473 Potsdam, Germany 2 Nevada Bureau of Mines and Geology, University of Nevada, Reno, Nevada, 89557, USA e-mail: jolie@gfz-potsdam.de ABSTRACT In the framework of geothermal exploration campaigns 3D structural-geological modeling plays an important role in the understanding of geothermal systems. The focus for the Brady s geothermal system located in the Basin and Range province is on the identification of structural controls on fluid flow and permeability anisotropy. In addition to 3D structural-geological modeling, the applied exploration strategy also includes stress field analysis and surface geochemical surveys. We have used 1) detailed geological maps, 2) borehole data, 3) processed 2D seismic and gravity data, and 4) a digital elevation model as input parameters of the 3D model. Based on these data, four representative cross sections have been developed as a major input for a preliminary 3D geological model. Well logs are used to verify the stratigraphic structure between the cross sections. The major strike direction of the faults is NNE. Normal faulting is the dominant stress regime. Dip angles range from 45 to 80. The 3D model consists of eight geological units. The Mesozoic basement consists of granites and metamorphic rocks. Above, a sequence of Tertiary ash-flow tuffs, lacustrine sediments, and lava flows of different composition has been encountered. 3D structural models populated with geomechanical and stress data can help to delineate between dilational and shear zone both being prone for channeling fluids. In a later stage, stress data derived from fault plane analysis shall be integrated into the 3D structuralgeological model applying the slip and dilation tendency technique to estimate hydraulically active fault zones. These results shall be verified by surface gas measurements to understand the impact of individual faults on fluid flow. GEOLOGICAL SETTING The Basin and Range-province is characterized by E- W to WNW-ESE extensional motion of the lithospheric crust, which began in the Paleogene about 29-30 Ma ago and reached main extension in the middle Miocene (Eaton, 1982; Eaton, 1984; Fosdick and Colgan, 2010; McQuarrie and Wernicke, 2005). Due to the extension, an alternating pattern of parallel elongated mountain ridges and desert basins formed, which increased its original width since the Oligocene between 100 % and 200 % and resulted in crustal thinning (Proffett, 1977, Velasco et al., 2010, Eaton, 1984). The extension resulted in approximately N to NNE-striking fault zones. Various geothermal surface manifestations along normal faults give evidence to the existence of hydrothermal systems. In this study we focus on the Brady s geothermal field, which occurs along a normal fault and appears to be amagmatic. It is located at the in the Basin and Range province (Nevada, USA) approximately 80 km northeast of Reno (Faulds et al, 2010). FIELD MAPPING The entire area around the Desert Peak and Brady s geothermal systems has been mapped in detail at 1:24,000s scale (Faulds and Garside, 2003; Faulds et al., 2010, unpublished data). Together with well logs and the results of geophysical surveys, the geological map provides a major input dataset for the development of geological cross sections from which the 3D structural-geological model can be derived. DEVELOPMENT OF CROSS SECTIONS Four cross sections have been developed, based on mapped surface faults and rock units, 2D seismic reflection data, gravity data, and re-interpreted drill cuttings and core from wells in the study area (Faulds, unpublished data). Faults and horizons have been attributed in each cross section, which allows

the later merging of data points from the same fault plane or horizon. The cross sections are constructed to a depth of 1,000 m bsl, which forms the basis of the model. the modeling workflow started. The structural model includes only fault planes and no horizons. The location of geothermal surface manifestations helps to identify the structural controls of fluid circulation. Control points between the cross sections help to constrain the fault plane computation. The workflow consists of four different stages: 1) definition of model information, 2) stepwise implementation of fault sets, 3) fault tree building (hierarchy model), and 4) fault modeling. Figure 1: Study area with fault traces (red), well locations, and developed cross sections (B, C, D, E). METHODS The development of the 3D structural-geological model has been completed in two steps. First, a fault model has been calculated, which is followed by the development of a horizon model. The model has a size of 7.2x8.3x3.0 km. Its base lies within the Mesozoic basement. Figure 3: Development of a fault plane I: Input data from four geological cross sections and the fault trace map. Figure 4: Development of a fault plane II: Computed fault plane, based on the input data from Fig. 4. Figure 2: Process of data preparation for the development of a 3D structural model. The entire model is based on four geological cross sections and a geological map. Fault modeling In total, 61 fault planes have been implemented into the model. This required the consistent attribution of each fault in each cross section. The cross sections have been digitized, georeferenced, and edited before The implementation of faults has been accomplished in a stepwise approach. Faults have been separated into three different classes and priorities according to their length, fault plane surface, and hierarchy. Each fault has been classified in the hierarchy system. The hierarchy has been created according to the developed cross sections. Faults with the highest priority have been implemented first, as they will have the strongest effects on the model design. The second fault set was modeled according to the given hierarchy of set 1. Set 3 is the final fault set, which contains minor fault planes associated with faults of set 1&2. The separation of huge fault sets with complex interactions into different sets appeared to be a very successful approach. The majority of the faults have been implemented as dying faults ; only

a few faults cross the entire model. The selected gridding method for the fault planes is a 2D minimum tension gridding algorithm with trend control (DGI, 2009). The fault model will provide the basis for the stress field analysis. Each fault plane can be assessed according to their orientation within the present stress field. The stress field analysis includes the stress inversion from surface measurements of fault slip data, and the slip and dilation tendency analysis after Morris et al. (1996). Figure 5: Final fault model with 61 fault planes. The dominant strike direction of the fault system is NNE (~N30E). reduced by the calibration with borehole data. However, the most accurate models with highest resolution can only be achieved by using additional information such as high definition 3D seismic reflection data. For the Brady s model a simplification of the geological succession was necessary due to heterogeneous volcanic units. The horizons of the cross sections are digitized in the same manner as the fault planes. Each horizon was consistently attributed throughout all cross sections. For the correct display of the geological succession, horizons are imported stepwise bottom-up. This allows maximum control on the interpolation of each horizon. The calibration of all horizons was carried out for each fault block. It is suggested to conduct the editing perpendicular to the major strike direction. This method allows the maximum control on the correct modeling of structural displacements. Once the interpolation of a horizon has been successfully performed, the data points of the overlying horizon are imported and interpolated in the same manner. At this stage eight different geological units have been integrated into the presented model. Figure 7: Horizon model of the Brady s geothermal system. The top horizon represents the early to late Miocene. Figure 6: Results of the dilation tendency analysis on selected fault planes within the Brady s geothermal system. Yellow to red colors indicate an increased tendency of faults to dilate. Horizon modeling The process of 3D structural-geological modeling, also referred to as 3D mapping (Moeck et al., 2010), includes the import of stratigraphic horizons. Their assembly is arranged according to the developed fault model. 3D structural-geological modeling provides information at each point of the model within the defined z-range. The uncertainty of models can be Figure 8: 3D geological model with the proposed stimulation well 15-12. RESULTS The 3D structural-geological model of the Brady s geothermal system comprises 61 fault planes and

eight geological units. Due to the complex fault pattern, the model consists of 71 large fault blocks. This compartmentalization has effects on the hydraulic conductivity, as it can isolate or connect production or re-injection wells from each other. The dip angles of the faults vary from 45 to 80. The major strike direction is NNE. The known geothermal reservoir is located at a depth of 600-1,500 m below surface and lies within the Oligocene ash-flow tuffs and Miocene lava flows. Production wells target the NNE-striking Brady s fault, which has a dip range from 70-80. The Brady s fault occurs on the western limb of an extensional syncline. Some of the encountered stratigraphic units could also contribute to the reservoir in addition to the faults. The 3D model for the Brady s area achieved in this study will help to guide the proposed stimulation in Well 15-12 in terms of enhancing understanding of the subsurface faults and major stratigraphic packages that lie between the well and the known geothermal reservoir to the north. OUTLOOK 3D mapping and stress field analysis should be a standard approach in geothermal exploration. It improves the interpretability of already existing data sets. Much of the data can be derived through noninvasive surface methods, which is beneficial especially in the first stages of exploration campaigns, where minimal information is available. The results can be used as a basis for any further exploration work (e.g. dynamic modeling). Geochemical surveys Other helpful methods to verify and improve the accuracy of the 3D model are, for example, geochemical surveys. For that reason, diffuse degassing measurements are planned to be applied at the Brady s geothermal system to actually confirm where potential pathways of geothermal fluids are present at depth and if results match with the modeled fault system. This approach seems to be promising as the Brady System is an active geothermal system, where various gas emanations (e.g. CO 2 ) are likely to occur with high effluxes. Integration with the results of 3D modeling and stress field analysis would improve this exploration technique. The method of diffuse degassing measurements has already been used for various purposes worldwide, such as volcanic hazard analysis and volcano monitoring (Hernandez, 2001; Fridriksson, 2006). REFERENCES Dynamic Graphics, Inc. (2009), earthvision 8.0 User Guide. Eaton, G.P. (1982), The Basin and range province: Origin and tectonic significance. Ann. Rev. Earth Planet. Sci. 1982. 10:409-40. Eaton, G.P. (1984), The Miocene Great Basin of Western North America as an extending Back-Arc Region, In: R.L. Carlson and K. Kobayashi (Editors), Geodynamics of Back-arc Regions. Tectonophysics, 102: 275-295. Faulds, J.E., and Garside, L.J. (2003), Preliminary geologic map of the Desert Peak Brady geothermal fields, Churchill County, Nevada. Nevada Bureau of Mines and Geology Open-File Report 03-27. Faulds, J.E., Coolbaugh, M.F., Benoit, D., Oppliger, G., Perkins, M., Moeck, I., and Drakos, P. (2010), Structural controls of geothermal activity in the northern Hot Springs Mountains, western Nevada: The tale of three geothermal systems (Brady s, Desert Perk, and Desert Queen). Geothermal Resources Council Transactions, 34, 675-683. Figure 9: Work stages of the applied exploration approach. However, both methods cannot qualify and quantify the presence of geothermal fluids in the inferred reservoir. At that stage the degree of uncertainty still has to be reduced. This can only be reached by additional surveys, which will assist in the verification of the results. Fosdick, J.C., Colgan, J.P. (2010), Miocene extension in the East Range, Nevada: A two-stage history of normal faulting in the northern Basin and Range, Geological Society of America Bulletin 2008;120;1198-1213, doi: 10.1130/B26201.1. Fridriksson, T., Kristjánsson, B.R., Ármannsson, H., Margrétardóttir, E., Ólafsdóttir, S., Chiodini, G. (2006), CO2 emissions and heat flow through soil,

fumaroles, and steam heated mud pools at the Reykjanes geothermal area, SW Iceland. Applied Geochemistry, 21, p.1551 1569. Hernandez, P.A., Notsu, K. et al. (2001), Carbon dioxide degassing by advective flow from Usu volcano, Japan. Science, 292, p.83-86. McQuarrie, N. and Wernicke, B.P. (2005), An animated tectonic reconstruction of southwestern North America since 36 Ma, Geosphere; December 2005; v.1; no.3; p.147 172; doi: 10.1130/GES00016.1. Moeck, I., Hinz, N., Faulds, J., Bell, J., Kell-Hills, A., and Louie, J. (2010), 3D geological modeling as a new method in geothermal exploration: A case study from central Nevada. Geothermal Resources Council Transactions, 34, p.807-811. Morris A., Ferrill D.A., Henderson D.B. (1996). Slip tendency analysis and fault reactivation, Geology 24(3), 275-278. Proffett, J.M. (1977), Cenozoic geology of the Yerington district, Nevada, and implications for the nature and origin of Basin and Range faulting, Geological Society of America Bulletin, v.88, p. 247-266. Velasco, M.S., Bennett, R.A., Johnson, R.A., Hreinsdóttir, S. (2010), Subsurface fault geometries and crustal extension in the eastern Basin and Range Province, western U.S., Tectonophysics 488 (2010) p.131 142.