scholarly journals Hybrid Kinematic-Dynamic Approach to Seismic Wave-Equation Modeling, Imaging, and Tomography

2015 ◽  
Vol 2015 ◽  
pp. 1-8 ◽  
Author(s):  
Alexandr S. Serdyukov ◽  
Anton A. Duchkov
Keyword(s):  
Mechanical Properties ◽  
Wave Equation ◽  
Efficient Implementation ◽  
Equation Modeling ◽  
Computationally Efficient ◽  
Time Step ◽  
Near Surface ◽  
Structure Response ◽  
Seismic Waveforms ◽  
First Arrival

Estimation of the structure response to seismic motion is an important part of structural analysis related to mitigation of seismic risk caused by earthquakes. Many methods of computing structure response require knowledge of mechanical properties of the ground which could be derived from near-surface seismic studies. In this paper we address computationally efficient implementation of the wave-equation tomography. This method allows inverting first-arrival seismic waveforms for updating seismic velocity model which can be further used for estimating mechanical properties. We present computationally efficient hybrid kinematic-dynamic method for finite-difference (FD) modeling of the first-arrival seismic waveforms. At every time step the FD computations are performed only in a moving narrowband following the first-arrival wavefront. In terms of computations we get two advantages from this approach: computation speedup and memory savings when storing computed first-arrival waveforms (it is not necessary to make calculations or store the complete numerical grid). Proposed approach appears to be specifically useful for constructing the so-called sensitivity kernels widely used for tomographic velocity update from seismic data. We then apply the proposed approach for efficient implementation of the wave-equation tomography of the first-arrival seismic waveforms.

scholarly journals Traveltime calculations from frequency‐domain downward‐continuation algorithms

Geophysics ◽  
2003 ◽  
Vol 68 (4) ◽  
pp. 1380-1388 ◽  
Author(s):  
Changsoo Shin ◽  
Seungwon Ko ◽  
Wonsik Kim ◽  
Dong‐Joo Min ◽  
Dongwoo Yang ◽  
...  
Keyword(s):  
Wave Equation ◽  
Frequency Domain ◽  
Downward Continuation ◽  
Calculation Algorithm ◽  
Near Surface ◽  
Absolute Value ◽  
Head Waves ◽  
First Arrival ◽  
The One ◽  
Derivatives Of

We present a new, fast 3D traveltime calculation algorithm that employs existing frequency‐domain wave‐equation downward‐continuation software. By modifying such software to solve for a few complex (rather than real) frequencies, we are able to calculate not only the first arrival and the approximately most energetic traveltimes at each depth point but also their corresponding amplitudes. We compute traveltimes by either taking the logarithm of displacements obtained by the one‐way wave equation at a frequency or calculating derivatives of displacements numerically. Amplitudes are estimated from absolute value of the displacement at a frequency. By using the one‐way downgoing wave equation, we also circumvent generating traveltimes corresponding to near‐surface upcoming head waves not often needed in migration. We compare the traveltimes computed by our algorithm with those obtained by picking the most energetic arrivals from finite‐difference solutions of the one‐way wave equation, and show that our traveltime calculation method yields traveltimes comparable to solutions of the one‐way wave equation. We illustrate the accuracy of our traveltime algorithm by migrating the 2D IFP Marmousi and the 3D SEG/EAGE salt models.

scholarly journals Computationally Efficient Solution of a 2D Diffusive Wave Equation Used for Flood Inundation Problems

Water ◽  
2019 ◽  
Vol 11 (10) ◽  
pp. 2195 ◽  
Author(s):  
Wojciech Artichowicz ◽  
Dariusz Gąsiorowski
Keyword(s):  
Wave Equation ◽  
Parallel Implementation ◽  
Time Integration ◽  
Computational Time ◽  
System Of Nonlinear Equations ◽  
Computationally Efficient ◽  
Time Step ◽  
Finite Element Scheme ◽  
Numerical Errors ◽  
Floodplain Inundation

This paper presents a study dealing with increasing the computational efficiency in modeling floodplain inundation using a two-dimensional diffusive wave equation. To this end, the domain decomposition technique was used. The resulting one-dimensional diffusion equations were approximated in space with the modified finite element scheme, whereas time integration was carried out using the implicit two-level scheme. The proposed algorithm of the solution minimizes the numerical errors and is unconditionally stable. Consequently, it is possible to perform computations with a significantly greater time step than in the case of the explicit scheme. An additional efficiency improvement was achieved using the symmetry of the tridiagonal matrix of the arising system of nonlinear equations, due to the application of the parallelization strategy. The computational experiments showed that the proposed parallel implementation of the implicit scheme is very effective, at about two orders of magnitude with regard to computational time, in comparison with the explicit one.

INFERRING THE COMPOSITION AND MECHANICAL PROPERTIES OF THE NEAR-SURFACE OF CERES FROM EMPLACEMENT MODELING OF LAYERED/PANCAKE-LIKE EJECTA DEPOSITS

2017 ◽  
Author(s):  
Kynan H.G. Hughson ◽  
◽  
Christopher T. Russell ◽  
Britney E. Schmidt ◽  
Heather Chilton ◽  
...  
Keyword(s):  
Mechanical Properties ◽  
Near Surface

scholarly journals A Computationally Efficient Shallow Water Model for Mixed Cohesive and Non-Cohesive Sediment Transport in the Yangtze Estuary

Water ◽  
2021 ◽  
Vol 13 (10) ◽  
pp. 1435
Author(s):  
Peng Hu ◽  
Junyu Tao ◽  
Aofei Ji ◽  
Wei Li ◽  
Zhiguo He
Keyword(s):  
Sediment Transport ◽  
Yangtze Estuary ◽  
Water Model ◽  
Cohesive Sediments ◽  
Shallow Water Model ◽  
Computationally Efficient ◽  
Time Step ◽  
The Yangtze Estuary ◽  
Erosion Coefficient ◽  
Cohesive Sediment Transport

In this paper, a computationally efficient shallow water model is developed for sediment transport in the Yangtze estuary by considering mixed cohesive and non-cohesive sediment transport. It is firstly shown that the model is capable of reproducing tidal-hydrodynamics in the estuarine region. Secondly, it is demonstrated that the observed temporal variation of suspended sediment concentration (SSC) for mixed cohesive and non-cohesive sediments can be well-captured by the model with calibrated parameters (i.e., critical shear stresses for erosion/deposition, erosion coefficient). Numerical comparative studies indicate that: (1) consideration of multiple sediment fraction (both cohesive and non-cohesive sediments) is important for accurate modeling of SSC in the Yangtze Estuary; (2) the critical shear stress and the erosion coefficient is shown to be site-dependent, for which intensive calibration may be required; and (3) the Deepwater Navigation Channel (DNC) project may lead to enhanced current velocity and thus reduced sediment deposition in the North Passage of the Yangtze Estuary. Finally, the implementation of the hybrid local time step/global maximum time step (LTS/GMaTS) (using LTS to update the hydro-sediment module but using GMaTS to update the morphodynamic module) can lead to a reduction of as high as 90% in the computational cost for the Yangtze Estuary. This advantage, along with its well-demonstrated quantitative accuracy, indicates that the present model should find wide applications in estuarine regions.

scholarly journals Evaluating the Observed Log-Likelihood Function in Two-Level Structural Equation Modeling with Missing Data: From Formulas to R Code

Psych ◽  
2021 ◽  
Vol 3 (2) ◽  
pp. 197-232
Author(s):  
Yves Rosseel
Keyword(s):  
Structural Equation ◽  
Structural Equation Models ◽  
Likelihood Function ◽  
Likelihood Estimation ◽  
Missing At Random ◽  
Equation Modeling ◽  
Computationally Efficient ◽  
Log Likelihood ◽  
Efficient Expression ◽  
Special Case

This paper discusses maximum likelihood estimation for two-level structural equation models when data are missing at random at both levels. Building on existing literature, a computationally efficient expression is derived to evaluate the observed log-likelihood. Unlike previous work, the expression is valid for the special case where the model implied variance–covariance matrix at the between level is singular. Next, the log-likelihood function is translated to R code. A sequence of R scripts is presented, starting from a naive implementation and ending at the final implementation as found in the lavaan package. Along the way, various computational tips and tricks are given.

Channel waves in cross‐borehole data

Geophysics ◽  
1992 ◽  
Vol 57 (2) ◽  
pp. 334-342 ◽  
Author(s):  
Larry R. Lines ◽  
Kenneth R. Kelly ◽  
John Queen
Keyword(s):  
Seismic Velocity ◽  
Test Facility ◽  
Equation Modeling ◽  
Borehole Data ◽  
Geological Information ◽  
Channel Wave ◽  
First Arrival ◽  
Borehole Seismic ◽  
Seismic Velocity Contrasts ◽  
Ray Bending

Layered geological formations with large seismic velocity contrasts can effectively create channel waves in cross‐borehole seismic data. The existence of channel waves for such waveguides can be confirmed by ray tracing, wave equation modeling, and modal analysis. Channel wave arrivals are identified in cross‐borehole data recorded at Conoco’s Newkirk test facility. For these data, where velocity contrasts are about 2 to 1, tomography based on first arrival traveltimes, is limited due to problems with extreme ray bending and seismic shadow zones. However, it may be possible to extract geological information using channel wave information. The seismometer differencing method appears to be a promising approach for detecting waveguide boundaries by use of cross‐borehole data.

GPU-Accelerated Efficient Implementation of FDTD Methods With Optimum Time-Step Selection

2014 ◽  
Vol 50 (2) ◽  
pp. 477-480 ◽  
Author(s):  
Theodoros T. Zygiridis ◽  
Nikolaos V. Kantartzis ◽  
Theodoros D. Tsiboukis
Keyword(s):  
Efficient Implementation ◽  
Time Step ◽  
Optimum Time ◽  
Step Selection ◽  
Fdtd Methods
Sign in / Sign up

Export Citation Format

Share Document