Proceedings 24* Geothermal Workshop 2002

PUTAUAKI MT RESISTIVITY SURVEY

C.J.BROMLEY

Institute of Geological Nuclear Sciences, Taupo, N.Z.

SUMMARY - New magnetotelluric resistivity soundings, supplemented by Schlumberger traversing measurements, have been collected in the vicinity of (Mt ) to investigate the southern and eastern extent of the Geothermal Field. They show that this field is associated with a large area of low resistivity totalling about Approximately 50% of this area is administered by the Putauaki and Te Tahuna Putauaki Trusts, and only a small portion has been tested with exploration drill holes. The volcanic peak of Putauaki itself yrs old) is on the southern boundary of the field. An association is inferred between the magma body feeding this vent and the heat source for the geothermal system yrs old). An area of particularly low resistivities ohm-m at lkm depth) and cold geothermal gas vents geothermometry is located near the centre of the revised low resistivity anomaly, about 1.4 east of the hottest existing Kawerau geothermal wells From this data, a geothermal hydrological model is inferred: a deep upflow, with associated intense hydrothermal alteration, occurs on faults centered beneath the gas vents, and supplies lateral flows to the northwest into the area of present production at Kawerau, as well as to the northeast.

1. INTRODUCTION normal step faults cause this drop in basement depth towards the centre of the Whakatane Fluids from the Kawerau Geothermal Field have Graben (Nairn, 1986). However, Wood and been utilised for process steam (up to 300 Brathwaite (1999) show that the most likely and electricity production (up to by the orientation of local basement faults, based on Tasman Pulp and Paper Mill since 1957 contours of basement elevations in Kawerau (Bloomer, 1997). Two binary plants also generate drillholes, is about 80". Indeed, about 60" about 6 from brine. The existing was the dominant trend of the most accurately steam field is managed by Century Resources on measured faults in the area - an east-west zone of behalf of the New Zealand Government. The bulk en-echelon fractures that ruptured just north of of the drilling and subsequent fluid extraction Kawerau geothermal field during the 1987 took place in a 2 area at the north-western Edgecumbe earthquake. Furthermore, the gravity corner of the field near the largest area of thermal anomaly contours in the south-eastern sector of features (Onepu Springs) and close to the Kawerau (MacDonald et al, 1970) imply a Tarawera River (Figure 1). Temperature predominantly east-west trend to the down- contours, pressure contours and chloride faulted,high-density basement. concentrations of the discharged fluid from the existing boreholes all imply that the fluid source The location of basement faults has important (hot, high pressure, saline fluid upflow) is to the implications for further exploration drilling, and southeast of the drilled area, in the direction of for management of production and injection Putauaki (Mt Edgecumbe) (Allis et al., 1997). strategies at Kawerau, because of the lack of Two exploration wells drilled to the southeast of extensive permeable stratigraphic aquifers. Where the and were found these faults pass through basement they can, in to be very hot (up to but impermeable. places, be highly permeable (eg. but at Subsequent efforts concentrated on the more other places the structures are impermeable (eg. permeable outflow area of cooler temperatures to Because of large uncertainties in the northwest. determining strike and dip of fault zones, and in predicting their productivity, targeting of inferred The known subsurface geology consists of a faults at Kawerau has had variable success (Wood sequence of volcanics (rhyolites, ignimbrites, and Brathwaite, 1999). tuffs and sediments) overlying a basement of greywacke at about 1 depth near the mill, but Some permeable aquifers do occur in the deepening to the northwest, and shallowing to overlying volcanics to the northwest, particularly about 700 m depth to the southeast. It was at depths above about 400 m and at about 800 initially presumed that north-east trending (White et al., 1997). The shallower aquifer has

135 suffered from pressure drawdown of several bars and Te Tahuna Putauaki Trusts. The objective of causing inflows of cooler water from above. this work was to determine the resistivity- Calcite deposition in the well bore was also a structure of this sector, to resolve uncertainties problem, particularly in wells tapping the regarding the eastern boundary of the Kawerau shallower aquifer. The deeper zone has shown Geothermal Field, and to provide information to little, if any, evidence of pressure drawdown after assist with planning future exploration drillholes. more than 25 years of exploitation. This implies The scope of the work included DC and MT that the deeper zone is recharged from a highly resistivity measurements and gas analyses. permeable source of hot water, and that the entire field could probably sustain a much higher rate of DC resistivity measurements were made using the extraction if this source were to be tapped. Schlumberger array with spacings of 500 m and 1000m (giving nominal probing depths of up The natural surface heat flow from thermal to 500 m). This helped resolve uncertainty, activity is estimated (Allis et al., 1993) to be caused by insufficient data, regarding the eastern about 100 Matahina Ignimbrite, boundary of the Kawerau geothermal system. (which occurs between about 100 m and 200 m of Magneto-telluric soundings were then sited using the surface) appears to act as a “blanket” or the traversing data as a guide. These soundings aquiclude, preventing upflows of hot fluid to the provided resistivity information as a function of surface over most of the geothermal system (east depth (fiom near surface to about two of Tarawera River), but is absent to the west kilometres), and as a function of orientation. where the thermal features are located. An additional subsurface outflow of about 50 Gases discharging through a patch of bare ground possibly extends into the (sometimes a shallow lake), at NZMG co- groundwaters of the Rangitaiki plains to the north ordinates 2838433 6340600 about 1.4 and northeast, contributing to anomalously high east of borehole were sampled and chloride and boron concentrations in these waters. analysed. The purpose was to determine whether these gases are of geothermal origin, and to Putauaki (Mt Edgecumbe) is a relatively young obtain the deep temperature of origin from gas cone (mostly years old), compared to the geothermometry. age of the geothermal system years, Browne, 1979). It may be a recent extrusion from Several magneto-telluric soundings were a shallow magma body that has provided a heat conducted around Putauaki (Mt Edgecumbe) in source for the geothermal system (Christenson, order to investigate the deep resistivity structure 1997). Ejecta fiom three hydrothermal of the volcano and its relationship to the vents, located near Tarawera Road, north of geothermal system. Soundings were also located Putauaki, have been dated at 14,500 and 9,000 in the vicinity of several hydrothermal eruption years old. craters near one of the hottest Kawerau boreholes, KA26 (Figure 1). Early resistivity surveys (before 1970) of Kawerau field (MacDonald et al., 1970; Studt, 3. RESISTIVITY AND GAS ANALYSIS 1958) were shallow penetrating (to about 250 m RESULTS depth) and limited in the area surveyed. They suggested a field area of 10 centred on the The DC Schlumberger array measurements Tasman mill. This was the basis for the conducted in the vicinity of Lake near subsequent exploration drilling. Later the eastern boundary of Kawerau Geothermal measurements, using a deeper penetration Field, reveal apparent resistivities that Schlumberger traversing array = 1000 m), consistently decrease with depth spacings showed a much greater field area of at least 19 500 m to 1000 m). They also increase west with a major eastward extension of the to east (about 15 ohm-mto45ohm-m, at = low resistivity ground (Allis et al., 1000 m). In the vicinity of Te Teko, further 1993). However, the location and cause of the northeast, apparent resistivities are greater than north-eastem and southern resistivity boundaries, 50 ohm-m. Figure 1shows the revised resistivity in particular, remained uncertain. The new boundary of the Kawerau Geothermal Field, resistivity work, reported here, allows for an incorporating these new data. The north-eastern improved resistivity model and an inferred “boundary” is actually a zone of gradual increase hydrological model of the Kawerau field to be in apparent resistivity from 20 ohm-m to 70 ohm- constructed. m consistent with a model of outflowing, and gradually cooling, geothermal fluids. The total 2. GEOSCIENTIFIC PROGRAMME area of the Kawerau resistivity anomaly (defined as ohm-mat = 1000 m) is Additional geoscientific exploration on the south- approximately eastern side of Kawerau Geothermal Field was undertaken with the collaboration of the Putauaki The locations of all the magnetotelluric soundings

136 around Putauaki are shown in Figure 1. intense smectite alteration at temperatures of Wherever possible, the orthogonal sounding about occurring within porous volcanic arrays were oriented NE (Ex) and using units. The underlying moderate resistivity layer 50 m electric field dipoles and perpendicular (about 30 ohm-m) coincides with greywacke at magnetic coil antennas. A measurement high temperature but relatively low frequency range of 8 to 0.02 Hz provided permeability, and moderate type penetration depth range from near surface to alteration. In the central area of the anomaly, low about 3 Notch filtering at 50 Hz was resistivities ohm-m) are found at depths of required because of nearby high voltage 1.2 to 2.5 km. This is beyond the depth of the powerlines. Natural signal strengths in the nearest well and the cause is unknown, but frequency bands near 1 Hz and 2 were, as may be related to a local increase in alteration usual, very low, resulting in some noisy portions intensity associated with basement fracturing. of the sounding curves containing poorly correlateddata. The sampled gas vent (within the former Lake Waimangau), lies at the centre of this low Layered interpretation of the filtered and resistivity anomaly. It is aligned with a probable smoothed apparent resistivity curves was fault structure, oriented at that links other undertaken using an automated iterative process. nearby patches of hydrothermally altered ground. The best fitting models, although subject to some The results of the compositional analysis of this uncertainty, have been combined into cross- gas show that it is a typical geothermal gas sections (see Figure 2) in order to assist with sample. Concentrations are given below in geothermal interpretation. Several soundings (dry gas): show evidence of distortions, particularly at low frequencies, that are probably caused by 3- = 853, = 6.52, He = 0.031, = 6.2, dimensional structure. Several others = 0, =58.8, = 75.2, Ar = 0.52 evidence of static effects (frequency independent displacement of the sounding curves The gas geothermometer produces a in the Ex and Ey directions). These are minimum subsurface temperature of for caused by near surface inhomogeneities. To the deep source of this gas (Christenson, pers. resolve some of these uncertainties would require a more detailed MT survey of the area backed up by co-located soundings to account for 4. CONCLUSIONS static shift effects. The resistivity and gas chemistry results reported The cross-sections presented in Figure 2, here strongly support the conclusion that most of however, are reasonable interpretations of the Putauaki Trust block (approximately 10 smoothed subsurface resistivity structure. They and the Te Tahuna Putauaki blocks (1.3 are reveal several new features related to the gross contained within the boundaries of the Kawerau structure of the geothermal system. Sections A- Geothermal Field, which has a revised estimated and B-B' (Figure 2) show that the peak of total area of 22 A subsurface outflow of Putauaki is located above the boundary of the low cooling geothermal fluid is inferred to extend resistivity anomaly associated with the several kilometres to the north and north-east, geothermal system. On the northern side of the causing a more diffuse resistivity boundary in this volcano, the low resistivity layer extends from at direction. The low resistivity anomaly also least 2 depth to within about 100 m depth. extends to the south, as far as Putauaki itself. In Further northeast, (beyond the 20 ohm-m contour the vicinity of the four existing south-eastern of Figure 1) the low resistivity is confined to a Kawerau boreholes layer between about 200 m depth and about 1 where the greywacke basement has high depth, with an underlying resistive basement temperatures (over we found an (presumed to be relatively unaltered, and cool, intermediate basement resistivity (30 ohm-m). greywacke). This resistivity structure is But about 5 km to the north-east, the underlying consistent with a hydrological model of upwelling greywacke has a much higher resistivity (several geothermal fluids feeding a fluid hundred ohm-m). Here it presumably forms a outflow through permeable volcanic formations. cooler, unaltered basement to the northern lateral These fluids have probably generated intense outflow. hydrothermal clay alteration. The most interesting part of the resistivity To the north-west, it is possible to compare anomaly lies at its centre, near an area of cold gas interpreted resistivities along sections B-B and discharge (former Lake Waimangau). Several with nearby borehole geology and soundings in this area reveal very low resistivities temperatures and These (less than 5 ohm-m) at drillable depth range (1 comparisons suggest that the shallow low to 2.5 km). Based on simple extrapolation of resistivity anomaly, in this area, is caused by lateral temperature gradients, very high

137 **

Figure 1. Kawerau Geothermal Field: resistivity boundary, magneto-telluric resistivity (*) and cross-section locations, boreholes and gas vent

NE MT layered resistivity interpretation

0 lkm

SW NE C‘

- -1 30

-2km - 0 lkm 0 lkm

Figure 2. Putauaki magneto-telluric resistivity interpretation. Cross-sections A-A’, B-B and ,showing layered resistivity models and inferred hydrological flows. temperatures can also be expected in this area east G.F., White, S.P. 1993. The natural state of of The low resistivities possibly indicate Kawerau Geothermal Field. Proc. more favourable basement permeability (associated Geothermal Workshop., with highly fractured and intensely altered zones) than at KA26 and Allis, R.G., 1997. The natural state and response to development of Kawerau Geothermal Field, New A sounding located adjacent to one of the large Zealand. Geothermal Resources Council hydrothermal eruption craters reveals a localised Transactions Vol. 21 15. high resistivity layer that may indicate the presence of a relatively young and unaltered intrusive dyke Bloomer, A., 1997. Kawerau Geothermal or sill under the eruption vent. Development: A Case Study. GRC Transactions 21 Analysis of the gas from the Putauaki vent confirms that it is geothermal in origin and has a source Browne, P.R.L. 1979. Minimum age of Kawerau temperature of at least Geothermal Field, North Island, New Zealand. Jnl. Volcanology and Geotheml Research, Vol. Based on the area of the resistivity anomaly, and the production history of the existing development, the potential sustainable production capacity of the Christenson, B.W., 1997. Kawerau Geothermal Kawerau field is likely to be in excess of 100 Field: Geochemical Structure of the reservoir and (e). The new results, when combined with the its response to exploitation. Geothermal knowledge existing Kawerau wells, which Resources Council TransactionsVol. 21 p 17-25. indicate a south-eastern direction for the origin of the deep recharging fluids (high temperature, high MacDonald, W.J.P., Muffler, L.J.P., Dawson, G.B. pressure and high chloride), confirm that a deep 1970. Geophysical investigations at Kawerau. exploration borehole should be drilled in the area of DSIR Geophysics Division Report 62. the gas vents, where high temperatures and permeability can be anticipated at depths of 1 to 2 Nairn, I.A., 1986. Geology of the Kawerau Geothermal Field -results of drilling present. In “Kawerau Geothermal Field” DSIR 5. ACKNOWLEDGEMENTS Geothermal Report Number 10 Compiled by

Trustees of Putauaki and Te Tahuna Putauaki Trusts are gratefully thanked for their Studt, F.E. 1958. Geophysical reconnaissance at encouragement and cooperation with this project. Kawerau, New Zealand. Jnl. Geol. & Stewart Bennie of IGNS contributed significantly Geophysics, 219-246. to the DC and MT resistivity fieldwork. Lew Klyen collected the gas sample and Bruce Christenson White, S.P., M.J., M.J., 1997. assisted with interpretation of the gas analysis. This Models of the Kawerau Geothermal Reservoir work has been supported by a research grant G.R.C. TransactionsVol21. the Foundation of Research Science and Technology. Wood C.P., Brathwaite 1999, The basement at Kawerau Geothermal Field. Proc N.Z. 6. REFERENCES Geotheml Workshop,

Allis, R.G., Christenson, B.W., Nairn, Risk,

139