Final Report DOE Contract No. DE-FG36-04G014294 ICEKAP 2004: A Collaborative Joint Geophysical Imaging Project at Krafla and IDDP P.E. Malin, S.A. Onacha, E. Shalev Division of Earth and Ocean Sciences Nicholas School of the Environment Duke University Durham, NC 27708

In this final report, we discuss both theoretical and applied research resulting from our DOE project, ICEKAP 2004: A Collaborative Joint Geophysical Imaging Project at Krafla and IDDP. The abstract below begins with a general discussion of the problem we addressed: the location and characterization of “blind” geothermal resources using microearthquake and magnetotelluric measurements. The abstract then describes the scientific results and their application to the Krafla geothermal area in Iceland. The text following this abstract presents the full discussion of this work, in the form of the PhD thesis of Stephen A. Onacha. The work presented here was awarded the “Best Geophysics Paper” at the 2005 Geothermal Resources Council meeting, Reno. This study presents the modeling of buried fault zones using microearthquake and electrical resistivity data based on the assumptions that fluid-filled fractures cause electrical and seismic anisotropy and polarization. In this study, joint imaging of electrical and seismic data is used to characterize the fracture porosity of the fracture zones. P-wave velocity models are generated from resistivity data and used in locating microearthquakes. Fracture porosity controls fluid circulation in the hydrothermal systems and the intersections of fracture zones close to the heat source form important upwelling zones formore » hydrothermal fluids. High fracture porosity sites occur along fault terminations, fault-intersection areas and fault traces. Hydrothermal fault zone imaging using resistivity and microearthquake data combines high-resolution multi-station seismic and electromagnetic data to locate rock fractures and the likely presence fluids in high temperature hydrothermal systems. The depths and locations of structural features and fracture porosity common in both the MT and MEQ data is incorporated into a joint imaging scheme to constrain resistivity, seismic velocities, and locations of fracture systems. The imaging of the fault zones is constrained by geological, drilling, and geothermal production data. The objective is to determine interpretation techniques for evaluating structural controls of fluid circulation in hydrothermal systems. The conclusions are: • directions of MT polarization and anisotropy and MEQ S-splitting correlate. Polarization and anisotropy are caused by fluid filled fractures at the base of the clay cap. •Microearthquakes occur mainly on the boundary of low resistivity within the fracture zone and high resistivity in the host rock. Resistivity is lowest within the core of the fracture zone and increases towards the margins of the fracture zone. The heat source and the clay cap for the hydrothermal have very low resistivity of less than 5Ωm. •Fracture porosity imaged by resistivity indicates that it varies between 45-5% with most between 10-20%, comparable to values from core samples in volcanic areas in Kenya and Iceland. For resistivity values above 60Ωm, the porosity reduces drastically and therefore this might be used as the upper limit for modeling fracture porosity from resistivity. When resistivity is lower than 5Ωm, the modeled fracture porosity increases drastically indicating that this is the low resistivity limit. This is because at very low resistivity in the heat source and the clay cap, the resistivity is dominated by ionic conduction rather than fracture porosity. •Microearthquakes occur mainly above the heat source which is defined by low resistivity at a depth of 3-4.5 km at the Krafla hydrothermal system and 4-7 km in the Longonot hydrothermal system. •Conversions of S to P waves occur for microearthquakes located above the heat source within the hydrothermal system. Shallow microearthquakes occur mainly in areas that show both MT and S-wave anisotropy. •S-wave splitting and MT anisotropy occurs at the base of the clay cap and therefore reflects the variations in fracture porosity on top of the hydrothermal system. •In the Krafla hydrothermal system in Iceland, both MT polarization and MEQ splitting directions align with zones that have high fracture porosity below the clay cap. These zones coincide with fault zones trending in the NNE-SSW and NW-SE directions in otherwise uniform volcanic rocks and laterally continuous geology. The NW-SE orientation is parallel to the regional shear fractures while the NNE-SSW trending polarizations align parallel to rift zone fracture swarms. This suggest that correlations between MT polarizations and MEQ splitting may be related to fluid filled fractures. •In areas of high resistivity (60Ωm), the P-wave velocity approaches that of the rock matrix. •S-wave splitting polarization is determined from measurements of angles of rotation to get the optimum direction of polarization. •The use of MEQ and resistivity for imaging fractures requires that the MEQ data acquisition system be located close to the fracture zone.« less