Fengjiao Zhang

Application of seismic full waveform inversion to monitor CO2 injection: modelling and a real data example from the Ketzin site, Germany

Seismic full waveform inversion (waveform tomography) is a method to reconstruct the underground velocity field in high resolution using seismic data. The method was first introduced during the 1980’s and became computationally feasible during the late 1990’s when the method was implemented in the frequency domain. This work presents three case studies and one synthetic benchmark of full waveform inversion applications. Two of the case studies are focused on time-lapse cross-well and 2D reflection seismic data sets acquired at the Ketzin CO2 geological storage site. These studies are parts of the CO2SINK and CO2MAN projects. The results show that waveform tomography is more effective than traveltime tomography for the CO2 injection monitoring at the Ketzin site for the cross-well geometry. For the surface data sets we find it is difficult to recover the true value of the velocity anomaly due to the injection using the waveform inversion method, but it is possible to qualitatively locate the distribution of the injected CO2. The results agree well with expectations based upon conventional 2D CDP processing methods and more extensive 3D CDP processing methods in the area. A further investigation was done to study the feasibility and efficiency of seismic full waveform inversion for time-lapse monitoring of onshore CO2 geological storage sites using a reflection seismic geometry with synthetic data sets. The results show that waveform inversion may be a good complement to standard CDP processing when monitoring CO2 injection. The choice of method and strategy for waveform inversion is quite dependent on the goals of the time-lapse monitoring of the CO2 injection. The last case study is an application of the full waveform inversion method to two crooked profiles at the Forsmark site in eastern central Sweden. The main goal of this study was to help determine if the observed reflections are mainly due to fluid filled fracture zones or mafic sills. One main difficulty here is that the profiles have a crooked line geometry which corresponds to 3D seismic geometry, but a 2D based inversion method is being used. This is partly handled by a 3D to 2D coordinate projection method from traveltime inversion. The results show that these reflections are primarily due to zones of lower velocity, consistent with them being generated at water filled fracture zones.

Feasibility of utilizing wavelet phase to map the CO2 plume at the Ketzin pilot site, Germany

ABSTRACT Spectral decomposition is a powerful tool that can provide geological details dependent upon discrete frequencies. Complex spectral decomposition using inversion strategies differs from conventional spectral decomposition methods in that it produces not only frequency information but also wavelet phase information. This method was applied to a time‐lapse three‐dimensional seismic dataset in order to test the feasibility of using wavelet phase changes to detect and map injected carbon dioxide within the reservoir at the Ketzin carbon dioxide storage site, Germany. Simplified zero‐offset forward modelling was used to help verify the effectiveness of this technique and to better understand the wavelet phase response from the highly heterogeneous storage reservoir and carbon dioxide plume. Ambient noise and signal‐to‐noise ratios were calculated from the raw data to determine the extracted wavelet phase. Strong noise caused by rainfall and the assumed spatial distribution of sandstone channels in the reservoir could be correlated with phase anomalies. Qualitative and quantitative results indicate that the wavelet phase extracted by the complex spectral decomposition technique has great potential as a practical and feasible tool for carbon dioxide detection at the Ketzin pilot site.

The First Post-injection Seismic Monitor Survey at the Ketzin Pilot CO2 Storage Site: Results from Time-lapse Analysis

The injection of CO2 at the Ketzin pilot CO2 storage site started in June 2008 and ended in August 2013. During the 62 months of injection, a total amount of about 67 kt of CO2 was injected into a ...

High resolution seismic imaging at the planned tunnel entrance to the Forsmark repository for spent nuclear fuel, central Sweden

The Swedish Nuclear Waste Management Company (SKB) plans to build a repository for storage of high-level radioactive spent nuclear fuel at the Forsmark site in central Sweden at a depth of about 470 m. The planned repository will cover an area of about 3.6 km2 at this depth. Prior to beginning excavation and tunneling, some detailed geophysical surveys are being performed at the planned site. One of these was a refraction seismic survey to determine depth to bedrock in the vicinity of the planned access ramp. Two lines, each about 300 m long and spaced about 35 m apart, were acquired in August 2011. Since the bedrock topography is known to be highly variable, a close receiver (2 m) and source (6 m) spacing was required to map it. This close spacing allowed the data also to be treated as reflection seismic data and some adjustments to the acquisition procedure were made in the field with this in mind to aid in the later processing. The main adjustment was that seismic data were recorded on all geophone stations simultaneously. That is, as shots were fired along one line, data were recorded along both that line and the other one. Likewise, when shots were fired along the other line, data were recorded along that line and the first line. This adjustment allowed semi-3D coverage between the lines. Results from first break traveltime tomography along the lines indicate a depth to bedrock that is greater than that found from geotechnical observations along the lines. This discrepancy is attributed to the uppermost bedrock being highly fractured and having a velocity significantly below that expected from the intact bedrock deeper down. Reflection seismic processing of the data shows a reflection at about 20 ms (about 60 m). The reflection is interpreted to have a gentle northwesterly dip component to it. Comparison with core data in the area suggests that the reflection is from a thin (a few metres thick) fracture zone, although none of the boreholes actually penetrate the reflector where it is mapped by the seismic data. This fracture zone may be part of a larger fracture zone mapped by core drilling further to the east. The newly mapped reflector may be crossed by the ramp when excavation begins. Further seismic surveying towards the west is required to verify if this will be the case.

Application of seismic waveform tomography to monitoring of CO2 injection: modeling and a real data example from the Ketzin site, Germany

Seismic full waveform inversion (waveform tomography) is a method to reconstruct the underground velocity field in high resolution using seismic data. The method was first introduced during the 1980’s and became computationally feasible during the late 1990’s when the method was implemented in the frequency domain. This work presents three case studies and one synthetic benchmark of full waveform inversion applications. Two of the case studies are focused on time-lapse cross-well and 2D reflection seismic data sets acquired at the Ketzin CO2 geological storage site. These studies are parts of the CO2SINK and CO2MAN projects. The results show that waveform tomography is more effective than traveltime tomography for the CO2 injection monitoring at the Ketzin site for the cross-well geometry. For the surface data sets we find it is difficult to recover the true value of the velocity anomaly due to the injection using the waveform inversion method, but it is possible to qualitatively locate the distribution of the injected CO2. The results agree well with expectations based upon conventional 2D CDP processing methods and more extensive 3D CDP processing methods in the area. A further investigation was done to study the feasibility and efficiency of seismic full waveform inversion for time-lapse monitoring of onshore CO2 geological storage sites using a reflection seismic geometry with synthetic data sets. The results show that waveform inversion may be a good complement to standard CDP processing when monitoring CO2 injection. The choice of method and strategy for waveform inversion is quite dependent on the goals of the time-lapse monitoring of the CO2 injection. The last case study is an application of the full waveform inversion method to two crooked profiles at the Forsmark site in eastern central Sweden. The main goal of this study was to help determine if the observed reflections are mainly due to fluid filled fracture zones or mafic sills. One main difficulty here is that the profiles have a crooked line geometry which corresponds to 3D seismic geometry, but a 2D based inversion method is being used. This is partly handled by a 3D to 2D coordinate projection method from traveltime inversion. The results show that these reflections are primarily due to zones of lower velocity, consistent with them being generated at water filled fracture zones.

Preliminary Seismic Time-lapse Results from the First Post-injection Survey at the Ketzin Pilot Site

Since the CO2SINK project start in April 2004, different methods involving seismics, geoelectrics andpressure-temperature monitoring have been applied at the Ketzin pilot site to map geological str ...

Feasibility Study of Seismic Interferometry for CO2 Storage Monitoring - An Example from the Ketzin Site

Seismic interferometry is a relatively low cost method compared with conventional seismic monitoring methods and can be performed together with microseismic and reservoir monitoring. These features make it a new potential tool for carbon dioxide storage monitoring. In this study, we acquired 6 nights of ambient noise data were recorded at the Ketzin experimental CO2 injection Site located in Ketzin, west of Berlin, Germany, in August 2013. An active survey was also acquired at the same time along of one line for a source test. Passive seismic interferometry was applied to the recorded noise data to reconstruct common shot gathers. The virtual shot gathers was processed to obtain a stacked section. The results show that the passive stacked sections are less coherent and low resolution compared with the active stacked section. However, we could also find some similar features to both the passive and active stacked sections in the shallow parts and even some agreement in the deeper parts.

Application of seismic full waveform inversion to monitor CO2injection: modelling and a real data example from the Ketzin site, Germany: Application of seismic full waveform inversion to monitor CO2injection

Seismic full waveform inversion (waveform tomography) is a method to reconstruct the underground velocity field in high resolution using seismic data. The method was first introduced during the 1980’s and became computationally feasible during the late 1990’s when the method was implemented in the frequency domain. This work presents three case studies and one synthetic benchmark of full waveform inversion applications. Two of the case studies are focused on time-lapse cross-well and 2D reflection seismic data sets acquired at the Ketzin CO2 geological storage site. These studies are parts of the CO2SINK and CO2MAN projects. The results show that waveform tomography is more effective than traveltime tomography for the CO2 injection monitoring at the Ketzin site for the cross-well geometry. For the surface data sets we find it is difficult to recover the true value of the velocity anomaly due to the injection using the waveform inversion method, but it is possible to qualitatively locate the distribution of the injected CO2. The results agree well with expectations based upon conventional 2D CDP processing methods and more extensive 3D CDP processing methods in the area. A further investigation was done to study the feasibility and efficiency of seismic full waveform inversion for time-lapse monitoring of onshore CO2 geological storage sites using a reflection seismic geometry with synthetic data sets. The results show that waveform inversion may be a good complement to standard CDP processing when monitoring CO2 injection. The choice of method and strategy for waveform inversion is quite dependent on the goals of the time-lapse monitoring of the CO2 injection. The last case study is an application of the full waveform inversion method to two crooked profiles at the Forsmark site in eastern central Sweden. The main goal of this study was to help determine if the observed reflections are mainly due to fluid filled fracture zones or mafic sills. One main difficulty here is that the profiles have a crooked line geometry which corresponds to 3D seismic geometry, but a 2D based inversion method is being used. This is partly handled by a 3D to 2D coordinate projection method from traveltime inversion. The results show that these reflections are primarily due to zones of lower velocity, consistent with them being generated at water filled fracture zones.