ESTIMATION OF DEEP SUBSURFACE STRUCTURE IN EUROPEAN HOT DRY ROCK TEST SITE , SOULTZ-SOUS-FORÊTS , FRANCE , BY USE OF THE AE REFLECTION METHOD
N. Soma,H. Niitsuma,R. Baria
Abstract:Recently, the further Hot Dry Rock (HDR) development toward great depth has been started at the Soultz-sous-Forêts, France. Therefore the understandings deeper subsurface structure is taking on more important role. In this paper we present an estimation of deep subsurface structure in the Soultz HDR site by applying a type of reflection technique which we call the AE reflection method in timefrequency domain, where induced acoustic emission (AE) signals are used as wavesource. First, we mention the AE reflection method in time-frequency domain briefly. Next, we show the three-dimensional deep imaging which includes some meaningful reflectors. Then we precisely relate detected reflectors along the trajectory of well GPK-1 and deepened GPK-2, and compare the AE reflection method with other borehole observations. The results indicate a good agreement with distribution of fracture zone in GPK-1 and corresponds with possible fracture zone from a synthetic petrographic log derived from microscopic examination of cuttings during deepening of GPK-2. INTRODUCTION The Soultz-sous-Forêts Hot Dry Rock (HDR) site is located in northeast France, Alsace, about 50 km north of Strasbourg (Fig. 1), and has been developed since 1987 with supports by France, Germany, and European Commission (EC) (Baria et al., 1995). An artificial geothermal reservoir was created by a hydraulic fracturing in the depth of 2800 3500 m from well GPK-1, and a second well GPK-2 was drilled and succeeded in developing a circulation system in 1995 (Baumgärtner et al., 1996). In summer and autumn 1997 a 4 months long term circulation test was performed. The result demonstrated a great possibility of HDR system as one of future actual energy resources because it showed a enough production rate, negligible water losses, a minimum input energy, and etc., although the temperature of hot water was not enough for real production of electricity (Baumgärtner et al., 1998). Fig. 1. Location map of the Soultz-sous-Forêts site For a practical and industrial production of electricity by using Soultz HDR system, it is necessary to obtain higher temperature that may be around 200 °C. In order to develop a scientific pilot plant of HDR system, the well GPK-2 was deepened toward 5000 m in spring of 1999. Therefore, it becomes important to investigate deeper area below the developed artificial reservoir. A measurement surrounding pre-existing geothermal reservoir is generally useful not only in an HDR site but also in conventional geothermal fields, because there are reports of an existence of massive deep resources below a conventionally used shallow reservoir (Muraoka et al., 1993). However, there is no easy way to extract useful information under the geothermal condition, such as high temperature, high pressure, and often great deep. For instance, considering seismic reflection survey, the high temperature and high pressure of the geothermal field prevent us from using powerful downhole artificial seismic source such as explosives. Furthermore in the Soultz HDR site, the developed area lies in high velocity basement crystalline rock mass overlaid by sedimentary rock, which prevents wave energy to deep important region. We developed one kind of reflection technique where acoustic emission (AE; that is often called microseismicity) signals are used as wavesource which we call the AE reflection method (Soma and Niitsuma, 1997; Soma et al., 1997). Recently we have enhanced its performance to obtain precise and high-resolution deep image by developing timefrequency domain analysis with the wavelet transform (Soma, 1998). In this paper, we will describe the estimation of deep subsurface structure in the Soultz HDR site by the AE reflection method. First, we introduce the AE reflection method in time-frequency domain briefly. Then we show the application to the Soultz HDR site using AE events which were observed during a massive hydraulic fracturing in 1993. We discuss a comparison between estimates by the AE reflection method with other borehole observations which include new data from the well GPK-2 deepened in 1999. AE REFLECTION METHOD IN TIMEFREQUENCY DOMAIN Detection of reflected waves In the AE application to geothermal field, we usually can not detect reflected waves from the raw observed AE waveforms directly because they are covered by complex coda; randomly scattered wave (Soma and Niitsuma, 1997; Soma et al., 1997). Therefore, we developed a method to detect reflected waves by examining the shape of three-dimensional hodogram which is a trace of the particle motion associated with the wave arrivals. It is known that the shape of the hodogram changes corresponding to the wave condition (Nagano et al., 1986). The shape of the hodogram is like spherical when encounters incoherent signal such as randomly scattered wave, on the other hand, it changes linear at that moment when the coherent signals such as direct Pand Swaves are detected. Therefore, we can detect reflected waves by analyzing the linearity of the hodogram under the assumption that reflected waves are coherent signal. The concept of change of linearity in hodogram is shown in Fig. 2. Fig. 2. Concept of linearity of hodogram For the evaluation of the shape of the hodogram in Time-Frequency domain, we use a covariance matrix method by using the wavelet transform. The wavelet transform can be regarded as one of time-frequency signal representation (Rioul and Vetterli, 1991), that shows smoothed time varied spectral which is similar to short time Fourier transform (STFT). Hence we can define a covariance matrix of the hodogram in time-frequency (shift-scale: terms in wavelet transform) domain as formula (1),