Monitoring the Groningen gas reservoir by noise interferometry
Since 2013, seismicity in the Groningen gas field in the Netherlands is monitored by geophone strings in two boreholes at reservoir level (3 km). Here we present ambient noise cross- correlation results for one month of data for one of the boreholes.
The cross-correlations show that the noise comes mainly from above. The dominance of cultural noise from the surface is inferred from the asymmetry of the cross-correlations as well as from the observed diurnal and weekly variations. The P and S wave velocity structure of the reservoir is accurately determined from inter-geophone travel time data. The P velocity profile shows good agreement with well log data and has an uncertainty of less than 5%.
Groningen gas field and borehole data
Fig.2. P velocity model and borehole monitoring system (NAM).
DC_T 3040m RO_T 269m
SDM-1
1 2 3 4 5 6 7 8 9 10
Fig.1. Left: The Groningen gas field in the Netherlands. Right: Two years of seismicity (ML 0.1-3.5), subsidence (cm) and two borehole locations (NAM).
The Groningen gas field in the Netherlands is one of the world’s largest onshore gas fields and has been producing from 1963. As a result of gas extraction, subsidence and induced seismicity occur which cause damage and concern in the area.
Since 2013, the year with the highest level of induced seismicity, two geophone strings have been used to monitor the reservoir. For borehole SDM-1, 10 15-Hz geophones are positioned vertically from the top to the bottom of the reservoir with a geophone spacing of 30 m.
We used 1-bit normalization and spectral whitening, together with a bandpass filter from 3 to 400 Hz. After that, for each station pair, the cross-correlation function was calculated for 6 seconds segments with 2/3 overlap. These segmented cross-correlations were stacked over 1 hour.
It is concluded that noise interferometry can be used to determine the seismic velocity structure from deep borehole data. In the future we will try to monitor seismic velocity variations in the reservoir.
Finally, 24(hours)*33(days) segments were obtained from each station pair. In total, for the 10 geophones, 55 independent cross-correlation pairs were calculated for each of the 3 components.
Wen Zhou (w.zhou@uu.nl) & Hanneke Paulssen
Department of Earth Sciences, Utrecht University, The Netherlands
Because the distributions of these (24*33) apparent travel times is skewed, we used kernel density estimations to obtain probability density functions (PDFs) of the travel times.
Data processing and cross-correlation functions
Noise variations
Fig.5. Positive part of the CCFs between the top and bottom geophone for the vertical component (left) and the east component (right). Black dots indicate estimated travel times (maximum of CCFs).
P and S wave velocity estimation
Fig.6. Kernel density estimations for P wave (left) and S wave travel times (right) from the top to the bottom geophone. Histograms are shown in blue. The red curves represent the estimated PDFs.
Both diurnal and weekly variations of the CCFs are observed for the direct downgoing wave. This indicates that during the day the noise is dominated by cultural noise from the surface.
1-bit Normalization Spectral whitening Filter 3-400Hz
Cross- correlation
Stacking over 1 hour
data Fig.3. Processing schedule.
The vertical component cross-correlation functions (CCFs) for the top geophone show a clear downgoing P wave. The east component CCFs do not only show the downgoing S wave but also the downgoing P wave.
Acknowledgements
We thank the NAM and Remco Romijn for providing us with the data. This project is funded from the European Union's Horizon 2020 research and innovation programme under the Marie Sklodowska-Curie grant agreement No 642029 – ITN CREEP.
P wave velocity (m/s)
3000 4000 5000 6000
Depth (m)
2700 2750 2800 2850 2900 2950 3000 3050 3100
logging VP geophone location VP estimation
S wave velocity (m/s)
1500 2000 2500 3000 3500 4000 4500
Depth (m)
2700 2750 2800 2850 2900 2950 3000 3050 3100
NAM VS average model geophone location VS estimation
Fig.4. Left: Vertical component cross-correlations with top geophone. Right: East component cross-correlations. Blue and red lines indicate P and S wave arrivals.
Fig.7. Left: Estimated P wave velocity profile in reservoir (blue) together with well log data. Right: S wave velocity profile (red) with average shear velocity of reservoir. Well log data (green) were provided by the NAM.
Delay time (s)
-0.2 -0.15 -0.1 -0.05 0 0.05 0.1 0.15 0.2
Geophone depth (m)
2751 2780 2809 2838
2868 2898
2927 2957
2987 3017
Delay time (s)
-0.2 -0.15 -0.1 -0.05 0 0.05 0.1 0.15 0.2
Geophone depth (m)
2751 2780 2809
2838
2868 2898
2927 2957
2987
3017
P S
The amplitude as well as the timing of the direct wave varies. The apparent P wave travel time for each geophone pair is measured from the maximum of the vertical component cross-correlation for each of the hourly stacks. The apparent S wave travel time is measured from the east component CCFs.
Delay time (s)
0 0.05 0.1 0.15 0.2
Number of samples
0 10 20 30 40 50
0 0.05 0.1 0.15 0.2
Probability density function
0 10 20 30 40 50 60
t=0.111s
Delay time (s)
0 0.02 0.04 0.06 0.08 0.1
Number of samples
0 20 40 60 80 100 120 140
0 0.02 0.04 0.06 0.08 0.1
Probability density function
0 50 100 150
t=0.074s
The maximum likelihood travel times of all geophone pairs were used to estimate the inter-geophone P wave and S wave velocities. The obtained P wave velocity profile has errors of less than 5% and shows good agreement with the P velocity profile obtained from well logging. The S wave velocity profile also agrees with the inferred average S wave velocity of the reservoir.
date-time of records
24/11/13 01/12 08/12 15/12 22/12
Delay time (s)
0 0.01 0.02 0.03 0.04 0.05 0.06 0.07 0.08 0.09
0.1 -0.05
-0.04 -0.03 -0.02 -0.01 0 0.01 0.02 0.03 0.04 0.05
date-time of records
24/11/13 01/12 08/12 15/12 22/12
Delay time (s)
0 0.02 0.04 0.06 0.08 0.1 0.12 0.14 0.16 0.18
0.2 -0.05
-0.04 -0.03 -0.02 -0.01 0 0.01 0.02 0.03 0.04 0.05