scholarly journals Seismotectonics and 1-D velocity model of the Greater Geneva Basin, France–Switzerland

2020 ◽  
Vol 221 (3) ◽  
pp. 2026-2047 ◽  
Author(s):  
Verónica Antunes ◽  
Thomas Planès ◽  
Jiří Zahradník ◽  
Anne Obermann ◽  
Celso Alvizuri ◽  
...  
Keyword(s):  
Large Scale ◽  
Waveform Inversion ◽  
Broad Band ◽  
Velocity Model ◽  
Strike Slip ◽  
P Wave ◽  
Tectonic Structures ◽  
Stress Regime ◽  
Geothermal Exploration ◽  
Velocity Models

SUMMARY The Greater Geneva Basin (GGB), located in southwestern Switzerland and neighboring France, is enclosed by the rotating northwestern edge of the Alpine front and the Jura mountains chain. Recently, this basin has received increasing attention as a target for geothermal exploration. Historical and instrumental seismicity suggest that faults affecting the basin may still be active. Moderate-magnitude earthquakes have been located along the Vuache fault, a major strike-slip structure crossing the basin. Before geothermal exploration starts, it is key to evaluate the seismic rate in the region and identify possible seismogenic areas. In this context, we deployed a temporary seismic network of 20 broad-band stations (from September 2016 to January 2018) to investigate the ongoing seismic activity, its relationship with local tectonic structures, and the large-scale kinematics of the area. Our network lowered the magnitude of completeness of the permanent Swiss and French networks from 2.0 to a theoretical value of 0.5. Using a new coherence-based detector (LASSIE - particularly effective to detect microseismicity in noisy environments), we recorded scarce seismicity in the basin with local magnitudes ranging from 0.7 to 2.1 ML. No earthquakes were found in the Canton of Geneva where geothermal activities will take place. We constructed a local ’minimum 1-D P-wave velocity model’ adapted to the GGB using earthquakes from surrounding regions. We relocated the events of our catalogue obtaining deeper hypocentres compared to the locations obtained using the available regional velocity models. We also retrieved eight new focal mechanisms using a combination of polarities and waveform inversion techniques (CSPS). The stress inversion shows a pure strike-slip stress regime, which is in agreement with structural and geological data. Combining the background seismicity with our catalogue, we identified seismogenic areas offsetting the basin.

A seismotectonic study and minimum 1D velocity model for the Greater Geneva Basin, Western Switzerland

2020 ◽  
Author(s):  
Verónica Antunes ◽  
Thomas Planès ◽  
Jiří Zahradník ◽  
Anne Obermann ◽  
Celso Alvizuri ◽  
...  
Keyword(s):  
Seismic Activity ◽  
Large Scale ◽  
Detection Threshold ◽  
Velocity Model ◽  
Strike Slip ◽  
Tectonic Deformation ◽  
Geothermal Exploration ◽  
Area Of Interest ◽  
The Impact ◽  
1D Velocity Model

<p>In the framework of the Geothermie2020 project, the canton of Geneva and the Industrial Services of Geneva (SIG) are currently developing geothermal exploration in the Greater Geneva Basin (GGB), located in south-western Switzerland and neighbouring France. Before geothermal exploration begins, it is important to investigate the ongoing seismic activity, its relationship with local tectonic features, and the large-scale kinematics of the area. Background seismicity suggest that the local tectonic structures affecting the basin may still be active. Moderate-magnitude earthquakes have been identified along the Vuache fault, a major strike-slip structure crossing the basin. In this context we deployed a dense temporary network of 20 broadband stations around and within the GGB, during ~1.5 years, and reaching a detection threshold 0.5M<sub>L</sub>. </p><p>Using a new coherence-based detector (LASSIE), we detected and located 158 events in our area of interest. However, only 20 events were located in the GGB, with local magnitudes ranging from 0.7 to 2.2M<sub>L</sub>. We found no earthquakes in the Canton of Geneva where geothermal activities are taking place. We constructed a local minimum 1D velocity model with VELEST, using the recorded seismicity together with earthquakes from adjacent regions, in a total of 1263 P- and S-picks. The new velocity model allowed to relocate micro-seismic activity up to 11km depth along the main fault systems (i.e. Vuache, Cruseilles, Le Coin, and Arve) offsetting the GGB. We retrieved 8 new focal mechanisms for the area, using a combination of polarities and waveform inversion techniques (CSPS method). A stress inversion shows a tectonic deformation dominated by a quasi-pure strike-slip regime in the GGB, consistent with structural and geological data.</p><p>The study of microseismicity in a quiet sedimentary basin is challenging due to the scarce occurrence of seismic events combined with low signal-to-noise ratios and the often strong attenuation. However, the investigation of the sporadic (yet present) natural seismicity with dedicated dense networks could provide useful information about the GGB, even with a short-term experiment. We propose a newly-computed 1D velocity model that can be used in the GGB for seismic monitoring purposes throughout the geothermal project. This model can be easily improved later on, whenever more data is available. Monitoring the evolution and dispersion of the seismic-activity through the identified seismogenic areas during the geothermal project is essential. Quantifying the seismic rate in the basin before geothermal operations start will help to quantify the impact that geothermal energy extraction might have on the GGB. </p>

Elastic reflection waveform inversion based on the decomposition of sensitivity kernels

Geophysics ◽  
2019 ◽  
Vol 84 (2) ◽  
pp. R235-R250 ◽  
Author(s):  
Zhiming Ren ◽  
Zhenchun Li ◽  
Bingluo Gu
Keyword(s):  
Wave Velocity ◽  
Waveform Inversion ◽  
Velocity Model ◽  
P Wave ◽  
Background Model ◽  
Reverse Time ◽  
S Wave ◽  
Velocity Models ◽  
Time Migration ◽  
Starting Model

Full-waveform inversion (FWI) has the potential to obtain an accurate velocity model. Nevertheless, it depends strongly on the low-frequency data and the initial model. When the starting model is far from the real model, FWI tends to converge to a local minimum. Based on a scale separation of the model (into the background model and reflectivity model), reflection waveform inversion (RWI) can separate out the tomography term in the conventional FWI kernel and invert for the long-wavelength components of the velocity model by smearing the reflected wave residuals along the transmission (or “rabbit-ear”) paths. We have developed a new elastic RWI method to build the P- and S-wave velocity macromodels. Our method exploits a traveltime-based misfit function to highlight the contribution of tomography terms in the sensitivity kernels and a sensitivity kernel decomposition scheme based on the P- and S-wave separation to suppress the high-wavenumber artifacts caused by the crosstalk of different wave modes. Numerical examples reveal that the gradients of the background models become sufficiently smooth owing to the decomposition of sensitivity kernels and the traveltime-based misfit function. We implement our elastic RWI in an alternating way. At each loop, the reflectivity model is generated by elastic least-squares reverse time migration, and then the background model is updated using the separated traveltime kernels. Our RWI method has been successfully applied in synthetic and real reflection seismic data. Inversion results demonstrate that the proposed method can retrieve preferable low-wavenumber components of the P- and S-wave velocity models, which are reliable to serve as a starting model for conventional elastic FWI. Also, our method with a two-stage inversion workflow, first updating the P-wave velocity using the PP kernels and then updating the S-wave velocity using the PS kernels, is feasible and robust even when P- and S-wave velocities have different structures.

scholarly journals Estimation of acoustic macro models using a genetic full-waveform inversion: Applications to the Marmousi model

Geophysics ◽  
2016 ◽  
Vol 81 (4) ◽  
pp. R173-R184 ◽  
Author(s):  
Angelo Sajeva ◽  
Mattia Aleardi ◽  
Eusebio Stucchi ◽  
Nicola Bienati ◽  
Alfredo Mazzotti
Keyword(s):  
Prior Information ◽  
Waveform Inversion ◽  
Velocity Model ◽  
Forward Modeling ◽  
Coarse Grid ◽  
P Wave ◽  
Stochastic Methods ◽  
Full Waveform Inversion ◽  
Velocity Models ◽  
Full Waveform

We have developed a stochastic full-waveform inversion that uses genetic algorithms (GA FWI) to estimate acoustic macro models of the P-wave velocity field. Stochastic methods such as GA severely suffer the curse of dimensionality, meaning that they require unaffordable computer resources for inverse problems with many unknowns and expensive forward modeling. To mitigate this issue, we have proposed a two-grid technique with a coarse grid to represent the subsurface for the GA inversion and a finer grid for the forward modeling. We have applied this procedure to invert synthetic acoustic data of the Marmousi model, and we have developed three different tests. The first two tests use a velocity model derived from standard stacking velocity analysis as prior information and differ only for the parameterization of the coarse grid. Their comparison indicates that a smart parameterization of the coarse grid may significantly improve the final result. The third test uses a linearly increasing 1D velocity model as prior information, a layer-stripping procedure, and a large number of model evaluations. All three tests return velocity models that fairly reproduce the long-wavelength structures of the Marmousi. First-break cycle skipping related to the seismograms of the final GA-FWI models is significantly reduced compared with that computed on the models used as prior information. Descent-based FWIs starting from final GA-FWI models yield velocity models with low and comparable model misfits and with an improved reconstruction of the structural details. The quality of the models obtained by GA FWI plus descent-based FWI is benchmarked against the models obtained by descent-based FWI started from a smoothed version of the Marmousi and started directly from the prior information models. Our results are promising and demonstrate the ability of the two-grid GA FWI to yield velocity models suitable as input to descent-based FWI.

scholarly journals Mitigating elastic effects in marine 3-D full-waveform inversion

2019 ◽  
Vol 220 (3) ◽  
pp. 2089-2104
Author(s):  
Òscar Calderón Agudo ◽  
Nuno Vieira da Silva ◽  
George Stronge ◽  
Michael Warner
Keyword(s):  
Wave Equation ◽  
Field Data ◽  
Waveform Inversion ◽  
Velocity Model ◽  
P Wave ◽  
Data Sets ◽  
Acoustic Wave Equation ◽  
Data Set ◽  
Velocity Models ◽  
Full Waveform

SUMMARY The potential of full-waveform inversion (FWI) to recover high-resolution velocity models of the subsurface has been demonstrated in the last decades with its application to field data. But in certain geological scenarios, conventional FWI using the acoustic wave equation fails in recovering accurate models due to the presence of strong elastic effects, as the acoustic wave equation only accounts for compressional waves. This becomes more critical when dealing with land data sets, in which elastic effects are generated at the source and recorded directly by the receivers. In marine settings, in which sources and receivers are typically within the water layer, elastic effects are weaker but can be observed most easily as double mode conversions and through their effect on P-wave amplitudes. Ignoring these elastic effects can have a detrimental impact on the accuracy of the recovered velocity models, even in marine data sets. Ideally, the elastic wave equation should be used to model wave propagation, and FWI should aim to recover anisotropic models of velocity for P waves (vp) and S waves (vs). However, routine three-dimensional elastic FWI is still commercially impractical due to the elevated computational cost of modelling elastic wave propagation in regions with low S-wave velocity near the seabed. Moreover, elastic FWI using local optimization methods suffers from cross-talk between different inverted parameters. This generally leads to incorrect estimation of subsurface models, requiring an estimate of vp/vs that is rarely known beforehand. Here we illustrate how neglecting elasticity during FWI for a marine field data set that contains especially strong elastic heterogeneities can lead to an incorrect estimation of the P-wave velocity model. We then demonstrate a practical approach to mitigate elastic effects in 3-D yielding improved estimates, consisting of using a global inversion algorithm to estimate a model of vp/vs, employing matching filters to remove elastic effects from the field data, and performing acoustic FWI of the resulting data set. The quality of the recovered models is assessed by exploring the continuity of the events in the migrated sections and the fit of the latter with the recovered velocity model.

scholarly journals Acoustic full-waveform inversion in an elastic world

Geophysics ◽  
2018 ◽  
Vol 83 (3) ◽  
pp. R257-R271 ◽  
Author(s):  
Òscar Calderón Agudo ◽  
Nuno Vieira da Silva ◽  
Michael Warner ◽  
Joanna Morgan
Keyword(s):  
Wave Velocity ◽  
Waveform Inversion ◽  
Velocity Model ◽  
P Wave ◽  
Data Sets ◽  
Full Waveform Inversion ◽  
Data Set ◽  
Velocity Models ◽  
Full Waveform ◽  
Elastic Data

Full-waveform inversion (FWI) is a technique used to obtain high-quality velocity models of the subsurface. Despite the elastic nature of the earth, the anisotropic acoustic wave equation is typically used to model wave propagation in FWI. In part, this simplification is essential for being efficient when inverting large 3D data sets, but it has the adverse effect of reducing the accuracy and resolution of the recovered P-wave velocity models, as well as a loss in potential to constrain other physical properties, such as the S-wave velocity given that amplitude information in the observed data set is not fully used. Here, we first apply conventional acoustic FWI to acoustic and elastic data generated using the same velocity model to investigate the effect of neglecting the elastic component in field data and we find that it leads to a loss in resolution and accuracy in the recovered velocity model. Then, we develop a method to mitigate elastic effects in acoustic FWI using matching filters that transform elastic data into acoustic data and find that it is applicable to marine and land data sets. Tests show that our approach is successful: The imprint of elastic effects on the recovered P-wave models is mitigated, leading to better-resolved models than those obtained after conventional acoustic FWI. Our method requires a guess of [Formula: see text] and is marginally more computationally demanding than acoustic FWI, but much less so than elastic FWI.

scholarly journals Machine Learning as a Seismic Prior Velocity Model Building Method for Full-Waveform Inversion: A Case Study from Colombia

2021 ◽  
Vol 178 (2) ◽  
pp. 423-448
Author(s):  
Ursula Iturrarán-Viveros ◽  
Andrés M. Muñoz-García ◽  
Octavio Castillo-Reyes ◽  
Khemraj Shukla
Keyword(s):  
Machine Learning ◽  
Initial Velocity ◽  
Waveform Inversion ◽  
Synthetic Data ◽  
Velocity Model ◽  
Wave Impedance ◽  
P Wave ◽  
Well Log ◽  
Velocity Models ◽  
Full Waveform

AbstractWe use machine learning algorithms (artificial neural networks, ANNs) to estimate petrophysical models at seismic scale combining well-log information, seismic data and seismic attributes. The resulting petrophysical images are the prior inputs in the process of full-waveform inversion (FWI). We calculate seismic attributes from a stacked reflected 2-D seismic section and then train ANNs to approximate the following petrophysical parameters: P-wave velocity ($$V_\mathrm{{p}}$$ V p ), density ($$\rho $$ ρ ) and volume of clay ($$V_\mathrm{{clay}}$$ V clay ). We extend the use of the $$V_\mathrm{{clay}}$$ V clay by constraining it with the well lithology and we establish two classes: sands and shales. Consequently, machine learning allows us to build an initial estimate of the earth property model ($$V_\mathrm{{p}}$$ V p ), which is iteratively refined to produce a synthetic seismogram that matches the observed seismic data. We apply the 1-D Kennett method as a forward modeling tool to create synthetic data with the images of $$V_\mathrm{{p}}$$ V p , $$\rho $$ ρ and the thickness of layers (sands or shales) obtained with the ANNs. A nonlinear least-squares inversion algorithm minimizes the residual (or misfit) between observed and synthetic full-waveform data, which improves the $$V_\mathrm{{p}}$$ V p resolution. In order to show the advantage of using the ANN velocity model as the initial velocity model for the inversion, we compare the results obtained with the ANNs and two other initial velocity models. One of these alternative initial velocity models is computed via P-wave impedance, and the other is achieved by velocity semblance analysis: root-mean-square velocity (RMS). The results are in good agreement when we use $$\rho $$ ρ and $$V_\mathrm{{p}}$$ V p obtained by ANNs. However, the results are poor and the synthetic data do not match the real acquired data when using the semblance velocity model and the $$\rho $$ ρ from the well log (constant for the entire 2-D section). Nevertheless, the results improve when including $$\rho $$ ρ , the layered structure driven by the $$V_\mathrm{{clay}}$$ V clay (both obtained with ANNs) and the semblance velocity model. When doing inversion starting with the initial $$V_\mathrm{{p}}$$ V p model estimated using the P-wave impedance, there is some gain of the final $$V_\mathrm{{p}}$$ V p with respect to the RMS initial $$V_\mathrm{{p}}$$ V p . To assess the quality of the inversion of $$V_\mathrm{{p}}$$ V p , we use the information for two available wells and compare the final $$V_\mathrm{{p}}$$ V p obtained with ANNs and the final $$V_\mathrm{{p}}$$ V p computed with the P-wave impedance. This shows the benefit of employing ANNs estimations as prior models during the inversion process to obtain a final $$V_\mathrm{{p}}$$ V p that is in agreement with the geology and with the seismic and well-log data. To illustrate the computation of the final velocity model via FWI, we provide an algorithm with the detailed steps and its corresponding GitHub code.

scholarly journals Earthquake location in tectonic structures of the Alpine Chain: the case of the Constance Lake (Central Europe) seismic sequence

2021 ◽  
Author(s):  
Francesca D’Ajello Caracciolo ◽  
Rodolfo Console
Keyword(s):  
Central Europe ◽  
Velocity Model ◽  
Moho Depth ◽  
Seismic Sequence ◽  
Tectonic Structures ◽  
Study Region ◽  
Waveform Data ◽  
Seismic Stations ◽  
Velocity Models ◽  
Master Event

AbstractA set of four magnitude Ml ≥ 3.0 earthquakes including the magnitude Ml = 3.7 mainshock of the seismic sequence hitting the Lake Constance, Southern Germany, area in July–August 2019 was studied by means of bulletin and waveform data collected from 86 seismic stations of the Central Europe-Alpine region. The first single-event locations obtained using a uniform 1-D velocity model, and both fixed and free depths, showed residuals of the order of up ± 2.0 s, systematically affecting stations located in different areas of the study region. Namely, German stations to the northeast of the epicenters and French stations to the west exhibit negative residuals, while Italian stations located to the southeast are characterized by similarly large positive residuals. As a consequence, the epicentral coordinates were affected by a significant bias of the order of 4–5 km to the NNE. The locations were repeated applying a method that uses different velocity models for three groups of stations situated in different geological environments, obtaining more accurate locations. Moreover, the application of two methods of relative locations and joint hypocentral determination, without improving the absolute location of the master event, has shown that the sources of the four considered events are separated by distances of the order of one km both in horizontal coordinates and in depths. A particular attention has been paid to the geographical positions of the seismic stations used in the locations and their relationship with the known crustal features, such as the Moho depth and velocity anomalies in the studied region. Significant correlations between the observed travel time residuals and the crustal structure were obtained.

Interferometric imaging condition for wave-equation migration

Geophysics ◽  
2008 ◽  
Vol 73 (2) ◽  
pp. S47-S61 ◽  
Author(s):  
Paul Sava ◽  
Oleg Poliannikov
Keyword(s):  
Large Scale ◽  
Wigner Distribution ◽  
Random Noise ◽  
Velocity Model ◽  
Distribution Functions ◽  
Small Scale ◽  
Interferometric Imaging ◽  
Imaging Data ◽  
Imaging Condition ◽  
Velocity Models

The fidelity of depth seismic imaging depends on the accuracy of the velocity models used for wavefield reconstruction. Models can be decomposed in two components, corresponding to large-scale and small-scale variations. In practice, the large-scale velocity model component can be estimated with high accuracy using repeated migration/tomography cycles, but the small-scale component cannot. When the earth has significant small-scale velocity components, wavefield reconstruction does not completely describe the recorded data, and migrated images are perturbed by artifacts. There are two possible ways to address this problem: (1) improve wavefield reconstruction by estimating more accurate velocity models and image using conventional techniques (e.g., wavefield crosscorrelation) or (2) reconstruct wavefields with conventional methods using the known background velocity model but improve the imaging condition to alleviate the artifacts caused by the imprecise reconstruction. Wedescribe the unknown component of the velocity model as a random function with local spatial correlations. Imaging data perturbed by such random variations is characterized by statistical instability, i.e., various wavefield components image at wrong locations that depend on the actual realization of the random model. Statistical stability can be achieved by preprocessing the reconstructed wavefields prior to the imaging condition. We use Wigner distribution functions to attenuate the random noise present in the reconstructed wavefields, parameterized as a function of image coordinates. Wavefield filtering using Wigner distribution functions and conventional imaging can be lumped together into a new form of imaging condition that we call an interferometric imaging condition because of its similarity to concepts from recent work on interferometry. The interferometric imaging condition can be formulated both for zero-offset and for multioffset data, leading to robust, efficient imaging procedures that effectively attenuate imaging artifacts caused by unknown velocity models.

Full-waveform inversion for full-wavefield imaging: Decades in the making

2021 ◽  
Vol 40 (5) ◽  
pp. 324-334
Author(s):  
Rongxin Huang ◽  
Zhigang Zhang ◽  
Zedong Wu ◽  
Zhiyuan Wei ◽  
Jiawei Mei ◽  
...  
Keyword(s):  
Model Building ◽  
Seismic Imaging ◽  
Structural Information ◽  
Waveform Inversion ◽  
Velocity Model ◽  
Data Sets ◽  
Full Waveform Inversion ◽  
Data Types ◽  
Velocity Models ◽  
Full Waveform

Seismic imaging using full-wavefield data that includes primary reflections, transmitted waves, and their multiples has been the holy grail for generations of geophysicists. To be able to use the full-wavefield data effectively requires a forward-modeling process to generate full-wavefield data, an inversion scheme to minimize the difference between modeled and recorded data, and, more importantly, an accurate velocity model to correctly propagate and collapse energy of different wave modes. All of these elements have been embedded in the framework of full-waveform inversion (FWI) since it was proposed three decades ago. However, for a long time, the application of FWI did not find its way into the domain of full-wavefield imaging, mostly owing to the lack of data sets with good constraints to ensure the convergence of inversion, the required compute power to handle large data sets and extend the inversion frequency to the bandwidth needed for imaging, and, most significantly, stable FWI algorithms that could work with different data types in different geologic settings. Recently, with the advancement of high-performance computing and progress in FWI algorithms at tackling issues such as cycle skipping and amplitude mismatch, FWI has found success using different data types in a variety of geologic settings, providing some of the most accurate velocity models for generating significantly improved migration images. Here, we take a step further to modify the FWI workflow to output the subsurface image or reflectivity directly, potentially eliminating the need to go through the time-consuming conventional seismic imaging process that involves preprocessing, velocity model building, and migration. Compared with a conventional migration image, the reflectivity image directly output from FWI often provides additional structural information with better illumination and higher signal-to-noise ratio naturally as a result of many iterations of least-squares fitting of the full-wavefield data.

Multi-parameter model building for the Qiuyue structure using four-component ocean-bottom seismometer data

Geophysics ◽  
2021 ◽  
pp. 1-52
Author(s):  
Yuzhu Liu ◽  
Xinquan Huang ◽  
Jizhong Yang ◽  
Xueyi Liu ◽  
Bin Li ◽  
...  
Keyword(s):  
Wave Velocity ◽  
Waveform Inversion ◽  
Velocity Model ◽  
P Wave ◽  
Velocity Analysis ◽  
Full Waveform Inversion ◽  
Ocean Bottom ◽  
Ocean Bottom Seismometer ◽  
Full Waveform ◽  
P Wave Velocity

Thin sand-mud-coal interbedded layers and multiples caused by shallow water pose great challenges to conventional 3D multi-channel seismic techniques used to detect the deeply buried reservoirs in the Qiuyue field. In 2017, a dense ocean-bottom seismometer (OBS) acquisition program acquired a four-component dataset in East China Sea. To delineate the deep reservoir structures in the Qiuyue field, we applied a full-waveform inversion (FWI) workflow to this dense four-component OBS dataset. After preprocessing, including receiver geometry correction, moveout correction, component rotation, and energy transformation from 3D to 2D, a preconditioned first-arrival traveltime tomography based on an improved scattering integral algorithm is applied to construct an initial P-wave velocity model. To eliminate the influence of the wavelet estimation process, a convolutional-wavefield-based objective function for the preprocessed hydrophone component is used during acoustic FWI. By inverting the waveforms associated with early arrivals, a relatively high-resolution underground P-wave velocity model is obtained, with updates at 2.0 km and 4.7 km depth. Initial S-wave velocity and density models are then constructed based on their prior relationships to the P-wave velocity, accompanied by a reciprocal source-independent elastic full-waveform inversion to refine both velocity models. Compared to a traditional workflow, guided by stacking velocity analysis or migration velocity analysis, and using only the pressure component or other single-component, the workflow presented in this study represents a good approach for inverting the four-component OBS dataset to characterize sub-seafloor velocity structures.

Sign in / Sign up

Export Citation Format

Share Document