Full-waveform analysis of core-mantle boundary structure using adjoint methods

Author(s):  
Anselme F.E. Borgeaud ◽  
Maria Koroni ◽  
Frédéric Deschamps

<div>This study presents a new approach for investigating the structure of the core-mantle boundary (CMB)</div><div>topography based on full-waveforms and adjoint methods. We compute intermediate period (10-20 seconds)</div><div>spectral-element seismograms using existing models of core-mantle boundary topography and we analyse the</div><div>sensitivity of relevant seismic phases. Our study adds new information about effects of CMB structure on</div><div>exact synthetics and observable traveltimes of seismic body waves by means of sensitivity kernels. It also</div><div>highlights the difficulty of imaging the boundary due to the strong trade-off between velocity and topography</div><div>variations, addressed by many previous investigators.</div><div> </div><div>Given the significance of CMB and its importance for many disciplines in geophysical research, there have</div><div>been many studies trying to understand and geographically map the variations of topography and velocity</div><div>above this boundary. However, the vast mantle area wherein seismic waves travel before and after they</div><div>interact with the CMB makes the identification of desired seismic phases somehow difficult. In addition, the</div><div>observable traveltimes can be hard to interpret as a result of the boundary’s topography only, due to the</div><div>approximate inverse methods and limited modelling methodologies. Despite considerable progress made the</div><div>past years, there is still a necessity for improving the understanding of effects of core-mantle boundary and</div><div>D″ structure on recorded waveforms.</div><div> </div><div>For our analyses, we perform comparisons between time shifts due to topography made on full-waveform</div><div>synthetics to ray theoretical predictions in order to assess methods usually deployed for imaging CMB.</div><div>Then, we calculate the corresponding sensitivity kernel for time windows around the theoretical arrival of</div><div>each phase. We focus on diffracted, core reflected and refracted <em>P</em> and <em>S</em> waves. The sensitivity kernels</div><div>depict the finite-frequency nature of these waves and possible contributions from other phases unpredictable</div><div>by ray theory. Results show that for most phases ray theory performs acceptably with some accuracy loss,</div><div>however comparisons of the effect of velocity variations to topography on traveltimes are discouraging due</div><div>to the low sensitivity to the latter.</div><div> </div><div>The conclusions drawn by our traveltime and sensitivity analyses are twofold. Firstly, using spectral-</div><div>element waveforms, the seismic phases which are frequently found in literature can be thoroughly investigated</div><div>and better understood, since their traveltime sensitivity through mantle and core is explicitly shown. The</div><div>full-waveform analysis allows us to assess the usability of phases which are informative for core-mantle</div><div>boundary structure and its topography. Secondly, we propose that using the analysed phases simultaneously</div><div>in a full-waveform inversion scheme will improve imaging of the CMB, while also allowing to jointly invert</div><div>for velocity variations along the D″ layer, which is generally poorly understood. From this study, we want</div><div>to promote advanced techniques of full-waveform inversion for improving CMB and lower mantle models.</div>

2021 ◽  
Author(s):  
Maria Koroni ◽  
Andreas Fichtner

<p>This study is a continuation of our efforts to connect adjoint methods and full-waveform inversion to common beamforming techniques, widely used and developed for signal enhancement. Our approach is focusing on seismic waves traveling in the Earth's mantle, which are phases commonly used to image internal boundaries, being however quite difficult to observe in real data. The main goal is to accentuate precursor waves arriving in well-known times before some major phase. These waves generate from interactions with global discontinuities in the mantle, thus being the most sensitive seismic phases and therefore most suitable for better understanding of discontinuity seismic structure. </p><p>Our work is based on spectral-element wave propagation which allows us to compute exact synthetic waveforms and adjoint methods for the calculation of sensitivity kernels. These tools are the core of full-waveform inversion and by our efforts we aim to incorporate more parts of the waveform in such inversion schemes. We have shown that targeted stacking of good quality waveforms arriving from various directions highlights the weak precursor waves. It additionally makes their traveltime finite frequency sensitivity prominent. This shows that we can benefit from using these techniques and exploit rather difficult parts of the seismogram.  It was also shown that wave interference is not easily avoided, but coherent phases arriving before the main phase also stack well and show on the sensitivity kernels. This does not hamper the evaluation of waveforms, as in a misfit measurement process one can exploit more phases on the body wave parts of seismograms.</p><p>In this study, we go a step forward and present recent developments of the approach relating to the effects of noise and a real data experiment. Realistic noise is added to synthetic waveforms in order to assess the methodology in a more pragmatic scenario. The addition of noise shows that stacking of coherent seismic phases is still possible and the sensitivity kernels of their traveltimes are not largely distorted, the precursor waves contribute sufficiently to their traveltime finite-frequency sensitivity kernels.<br>Using a well-located seismic array, we apply the method to real data and try to examine the possibilities of using non-ideal waveforms to perform imaging of the mantle discontinuity structure on the specific areas. In order to make the most out of the dense array configuration, we try subgroups of receivers for the targeted stacking and by moving along the array we aim at creating a cluster of stacks. The main idea is to use the subgroups as single receivers and create an evaluation of seismic discontinuity structure using information from each stack belonging to a subgroup. <br>Ideally, we aim at improving the tomographic images of discontinuities of selected regions by exploiting weaker seismic waves, which are nonetheless very informative.</p>


2016 ◽  
Vol 4 (4) ◽  
pp. SU17-SU24 ◽  
Author(s):  
Vanessa Goh ◽  
Kjetil Halleland ◽  
René-Édouard Plessix ◽  
Alexandre Stopin

Reducing velocity inaccuracy in complex settings is of paramount importance for limiting structural uncertainties, therefore helping the geologic interpretation and reservoir characterization. Shallow velocity variations due, for instance, to gas accumulations or carbonate reefs, are a common issue offshore Malaysia. These velocity variations are difficult to image through standard reflection-based velocity model building. We have applied full-waveform inversion (FWI) to better characterize the upper part of the earth model for a shallow-water field, located in the Central Luconia Basin offshore Sarawak. We have inverted a narrow-azimuth data set with a maximum inline offset of 4.4 km. Thanks to dedicated broadband preprocessing of the data set, we could enhance the signal-to-noise ratio in the 2.5–10 Hz frequency band. We then applied a multiparameter FWI to estimate the background normal moveout velocity and the [Formula: see text]-parameter. Full-waveform inversion together with broadband data processing has helped to better define the faults and resolve the thin layers in the shallow clastic section. The improvements in the velocity model brought by FWI lead to an improved image of the structural closure and flanks. Moreover, the increased velocity resolution helps in distinguishing between two different geologic interpretations.


Geophysics ◽  
2016 ◽  
Vol 81 (6) ◽  
pp. R385-R397 ◽  
Author(s):  
Christian Boehm ◽  
Mauricio Hanzich ◽  
Josep de la Puente ◽  
Andreas Fichtner

Adjoint methods are a key ingredient of gradient-based full-waveform inversion schemes. While being conceptually elegant, they face the challenge of massive memory requirements caused by the opposite time directions of forward and adjoint simulations and the necessity to access both wavefields simultaneously for the computation of the sensitivity kernel. To overcome this bottleneck, we have developed lossy compression techniques that significantly reduce the memory requirements with only a small computational overhead. Our approach is tailored to adjoint methods and uses the fact that the computation of a sufficiently accurate sensitivity kernel does not require the fully resolved forward wavefield. The collection of methods comprises reinterpolation with a coarse temporal grid as well as adaptively chosen polynomial degree and floating-point precision to represent spatial snapshots of the forward wavefield on hierarchical grids. Furthermore, the first arrivals of adjoint waves are used to identify “shadow zones” that do not contribute to the sensitivity kernel. Numerical experiments show the high potential of this approach achieving an effective compression factor of three orders of magnitude with only a minor reduction in the rate of convergence. Moreover, it is computationally cheap and straightforward to integrate in finite-element wave propagation codes with possible extensions to finite-difference methods.


Geophysics ◽  
2020 ◽  
Vol 85 (5) ◽  
pp. B169-B179
Author(s):  
Majid Mirzanejad ◽  
Khiem T. Tran ◽  
Michael McVay ◽  
David Horhota ◽  
Scott J. Wasman

Sinkhole collapse may result in significant property damage and even loss of life. Early detection of sinkhole attributes (buried voids, raveling zones) is critical to limit the cost of remediation. One of the most promising ways to obtain subsurface imaging is 3D seismic full-waveform inversion. For demonstration, a recently developed 3D Gauss-Newton full-waveform inversion (3D GN-FWI) method is used to detect buried voids, raveling soils, and characterize variable subsurface soil/rock layering. It is based on a finite-difference solution of 3D elastic wave equations and Gauss-Newton optimization. The method is tested first on a data set constructed from the numerical simulation of a challenging synthetic model and subsequently on field data collected from two separate test sites in Florida. For the field tests, receivers and sources are placed in uniform 2D surface grids to acquire the seismic wavefields, which then are inverted to extract the 3D subsurface velocity structures. The inverted synthetic results suggest that the approach is viable for detecting voids and characterizing layering. The field seismic results reveal that the 3D waveform analysis identified a known manmade void (plastic culvert), unknown natural voids, raveling, as well as laterally variable soil/rock layering including rock pinnacles. The results are confirmed later by standard penetration tests, including depth to bedrock, two buried voids, and a raveling soil zone. Our study provides insight into the application of the 3D seismic FWI technique as a powerful tool in detecting shallow voids and other localized subsurface features.


2015 ◽  
Vol 6 (2) ◽  
pp. 5-16 ◽  
Author(s):  
Sergio Alberto Abreo Carrillo ◽  
Ana B. Ramirez ◽  
Oscar Reyes ◽  
David Leonardo Abreo-Carrillo ◽  
Herling González Alvarez

Sign in / Sign up

Export Citation Format

Share Document