An efficient multiscale method for time-domain waveform tomography

Geophysics ◽  
2009 ◽  
Vol 74 (6) ◽  
pp. WCC59-WCC68 ◽  
Author(s):  
Chaiwoot Boonyasiriwat ◽  
Paul Valasek ◽  
Partha Routh ◽  
Weiping Cao ◽  
Gerard T. Schuster ◽  
...  

This efficient multiscale method for time-domain waveform tomography incorporates filters that are more efficient than Hamming-window filters. A strategy for choosing optimal frequency bands is proposed to achieve computational efficiency in the time domain. A staggered-grid, explicit finite-difference method with fourth-order accuracy in space and second-order accuracy in time is used for forward modeling and the adjoint calculation. The adjoint method is utilized in inverting for an efficient computation of the gradient directions. In the multiscale approach, multifrequency data and multiple grid sizes are used to overcome somewhat the severe local minima problem of waveform tomography. The method is applied successfully to 1D and 2D heterogeneous models; it can accurately recover low- and high-wavenumber components of the velocity models. The inversion result for the 2D model demonstrates that the multiscale method is computationally efficient and converges faster than a conventional, single-scale method.

Geophysics ◽  
2021 ◽  
pp. 1-50
Author(s):  
German Garabito ◽  
José Silas dos Santos Silva ◽  
Williams Lima

In land seismic data processing, the prestack time migration (PSTM) image remains the standard imaging output, but a reliable migrated image of the subsurface depends on the accuracy of the migration velocity model. We have adopted two new algorithms for time-domain migration velocity analysis based on wavefield attributes of the common-reflection-surface (CRS) stack method. These attributes, extracted from multicoverage data, were successfully applied to build the velocity model in the depth domain through tomographic inversion of the normal-incidence-point (NIP) wave. However, there is no practical and reliable method for determining an accurate and geologically consistent time-migration velocity model from these CRS attributes. We introduce an interactive method to determine the migration velocity model in the time domain based on the application of NIP wave attributes and the CRS stacking operator for diffractions, to generate synthetic diffractions on the reflection events of the zero-offset (ZO) CRS stacked section. In the ZO data with diffractions, the poststack time migration (post-STM) is applied with a set of constant velocities, and the migration velocities are then selected through a focusing analysis of the simulated diffractions. We also introduce an algorithm to automatically calculate the migration velocity model from the CRS attributes picked for the main reflection events in the ZO data. We determine the precision of our diffraction focusing velocity analysis and the automatic velocity calculation algorithms using two synthetic models. We also applied them to real 2D land data with low quality and low fold to estimate the time-domain migration velocity model. The velocity models obtained through our methods were validated by applying them in the Kirchhoff PSTM of real data, in which the velocity model from the diffraction focusing analysis provided significant improvements in the quality of the migrated image compared to the legacy image and to the migrated image obtained using the automatically calculated velocity model.


2020 ◽  
Vol 10 (11) ◽  
pp. 3750 ◽  
Author(s):  
Takumi Yoshida ◽  
Takeshi Okuzono ◽  
Kimihiro Sakagami

This paper presents a proposal of a time domain room acoustic solver using novel fourth-order accurate explicit time domain finite element method (TD-FEM), with demonstration of its applicability for practical room acoustic problems. Although time domain wave acoustic methods have been extremely attractive in recent years as room acoustic design tools, a computationally efficient solver is demanded to reduce their overly large computational costs for practical applications. Earlier, the authors proposed an efficient room acoustic solver using explicit TD-FEM having fourth-order accuracy in both space and time using low-order discretization techniques. Nevertheless, this conventional method only achieves fourth-order accuracy in time when using only square or cubic elements. That achievement markedly impairs the benefits of FEM with geometrical flexibility. As described herein, that difficulty is solved by construction of a specially designed time-integration method for time discretization. The proposed method can use irregularly shaped elements while maintaining fourth-order accuracy in time without additional computational complexity compared to the conventional method. The dispersion and dissipation characteristics of the proposed method are examined respectively both theoretically and numerically. Moreover, the practicality of the method for solving room acoustic problems at kilohertz frequencies is presented via two numerical examples of acoustic simulations in a rectangular sound field including complex sound diffusers and in a complexly shaped concert hall.


2021 ◽  
Author(s):  
Bojian Yin ◽  
Federico Corradi ◽  
Sander M. Bohté

ABSTRACTInspired by more detailed modeling of biological neurons, Spiking neural networks (SNNs) have been investigated both as more biologically plausible and potentially more powerful models of neural computation, and also with the aim of extracting biological neurons’ energy efficiency; the performance of such networks however has remained lacking compared to classical artificial neural networks (ANNs). Here, we demonstrate how a novel surrogate gradient combined with recurrent networks of tunable and adaptive spiking neurons yields state-of-the-art for SNNs on challenging benchmarks in the time-domain, like speech and gesture recognition. This also exceeds the performance of standard classical recurrent neural networks (RNNs) and approaches that of the best modern ANNs. As these SNNs exhibit sparse spiking, we show that they theoretically are one to three orders of magnitude more computationally efficient compared to RNNs with comparable performance. Together, this positions SNNs as an attractive solution for AI hardware implementations.


Geophysics ◽  
2014 ◽  
Vol 79 (3) ◽  
pp. R103-R119 ◽  
Author(s):  
Jianyong Bai ◽  
David Yingst ◽  
Robert Bloor ◽  
Jacques Leveille

Because of the conversion of elastic energy into heat, seismic waves are attenuated and dispersed as they propagate. The attenuation effects can reduce the resolution of velocity models obtained from waveform inversion or even cause the inversion to produce incorrect results. Using a viscoacoustic model consisting of a single standard linear solid, we discovered a theoretical framework of viscoacoustic waveform inversion in the time domain for velocity estimation. We derived and found the viscoacoustic wave equations for forward modeling and their adjoint to compensate for the attenuation effects in viscoacoustic waveform inversion. The wave equations were numerically solved by high-order finite-difference methods on centered grids to extrapolate seismic wavefields. The finite-difference methods were implemented satisfying stability conditions, which are also presented. Numerical examples proved that the forward viscoacoustic wave equation can simulate attenuative behaviors very well in amplitude attenuation and phase dispersion. We tested acoustic and viscoacoustic waveform inversions with a modified Marmousi model and a 3D field data set from the deep-water Gulf of Mexico for comparison. The tests with the modified Marmousi model illustrated that the seismic attenuation can have large effects on waveform inversion and that choosing the most suitable inversion method was important to obtain the best inversion results for a specific seismic data volume. The tests with the field data set indicated that the inverted velocity models determined from the acoustic and viscoacoustic inversions were helpful to improve images and offset gathers obtained from migration. Compared to the acoustic inversion, viscoacoustic inversion is a realistic approach for real earth materials because the attenuation effects are compensated.


Author(s):  
Anton Turk ◽  
Jasna Prpić-Oršić ◽  
Carlos Guedes Soares

A hybrid nonlinear time domain seakeeping analysis is applied to the study of a container ship advancing at different headings and encounter frequencies. A time-domain nonlinear strip theory in six degrees-of-freedom has been extended to predict ship motions by solving the unsteady hydrodynamic problem in the frequency domain and the equations of motion in the time domain which allows introducing nonlinearities in the linear model. The code is used to make parametric roll predictions for various speeds and headings and the results are summarized in a very intuitive 2D and 3D polar plots showing the full range of the parametric rolling realizations. The method developed is fairly accurate, robust, very computationally efficient, and can predict nonlinear ship motions. It is well suited to be used as a tool in ship design or as part of a path optimization model.


Sign in / Sign up

Export Citation Format

Share Document