Acoustic and elastic finite-difference modeling by optimal variable-length spatial operators

Geophysics ◽  
2020 ◽  
Vol 85 (2) ◽  
pp. T57-T70 ◽  
Author(s):  
Yang Liu
Keyword(s):  
Finite Difference ◽  
Finite Differences ◽  
Spatial Dispersion ◽  
Finite Difference Methods ◽  
High Order ◽  
Variable Length ◽  
Time Step ◽  
Heterogeneous Model ◽  
Fixed Length ◽  
Temporal Dispersion

Time-space domain finite-difference modeling has always had the problem of spatial and temporal dispersion. High-order finite-difference methods are commonly used to suppress spatial dispersion. Recently developed time-dispersion transforms can effectively eliminate temporal dispersion from seismograms produced by the conventional modeling of high-order spatial and second-order temporal finite differences. To improve the efficiency of the conventional modeling, I have developed optimal variable-length spatial finite differences to efficiently compute spatial derivatives involved in acoustic and elastic wave equations. First, considering that temporal dispersion can be removed, I prove that minimizing the relative error of the phase velocity can be approximately implemented by minimizing that of the spatial dispersion. Considering that the latter minimization depends on the wavelength that is dependent on the velocity, in this sense, this minimization is indirectly related to the velocity, and thus leads to variation of the spatial finite-difference operator with velocity for a heterogeneous model. Second, I use the Remez exchange algorithm to obtain finite-difference coefficients with the lowest spatial dispersion error over the largest possible wavenumber range. Then, dispersion analysis indicates the validity of the approximation and the algorithm. Finally, I use modeling examples to determine that the optimal variable-length spatial finite difference can greatly increase the modeling efficiency, compared to the conventional fixed-length one. Stability analysis and modeling experiments also indicate that the variable-length finite difference can adopt a larger time step to perform stable modeling than the fixed-length one for inhomogeneous models.

Time-space-domain implicit finite-difference methods for modeling acoustic wave equations

Geophysics ◽  
2018 ◽  
Vol 83 (4) ◽  
pp. T175-T193 ◽  
Author(s):  
Enjiang Wang ◽  
Jing Ba ◽  
Yang Liu
Keyword(s):  
Finite Difference ◽  
Acoustic Wave ◽  
Wave Equations ◽  
Finite Difference Methods ◽  
Second Order ◽  
Computational Time ◽  
Time Step ◽  
Space Domain ◽  
Time Space ◽  
Temporal Dispersion

It has been proved that the implicit spatial finite-difference (FD) method can obtain higher accuracy than explicit FD by using an even smaller operator length. However, when only second-order FD in time is used, the combined FD scheme is prone to temporal dispersion and easily becomes unstable when a relatively large time step is used. The time-space domain FD can suppress the temporal dispersion. However, because the spatial derivatives are solved explicitly, the method suffers from spatial dispersion and a large spatial operator length has to be adopted. We have developed two effective time-space-domain implicit FD methods for modeling 2D and 3D acoustic wave equations. First, the high-order FD is incorporated into the discretization for the second-order temporal derivative, and it is combined with the implicit spatial FD. The plane-wave analysis method is used to derive the time-space-domain dispersion relations, and two novel methods are proposed to determine the spatial and temporal FD coefficients in the joint time-space domain. First, we fix the implicit spatial FD coefficients and derive the quadratic convex objective function with respect to temporal FD coefficients. The optimal temporal FD coefficients are obtained by using the linear least-squares method. After obtaining the temporal FD coefficients, the SolvOpt nonlinear algorithm is applied to solve the nonquadratic optimization problem and obtain the optimized temporal and spatial FD coefficients simultaneously. The dispersion analysis, stability analysis, and modeling examples validate that the proposed schemes effectively increase the modeling accuracy and improve the stability conditions of the traditional implicit schemes. The computational efficiency is increased because the schemes can adopt larger time steps with little loss of spatial accuracy. To reduce the memory requirement and computational time for storing and calculating the FD coefficients, we have developed the representative velocity strategy, which only computes and stores the FD coefficients at several selected velocities. The modeling result of the 2D complicated model proves that the representative velocity strategy effectively reduces the memory requirements and computational time without decreasing the accuracy significantly when a proper velocity interval is used.

A High-Order Flux Reconstruction Method for 2-D Vorticity Transport

2021 ◽  
Author(s):  
Adrin Gharakhani
Keyword(s):  
Finite Difference ◽  
Poisson Equation ◽  
Finite Difference Methods ◽  
Neumann Boundary ◽  
Normal Derivative ◽  
High Order ◽  
Benchmark Problems ◽  
Flux Reconstruction ◽  
Reconstruction Method ◽  
Vorticity Transport

Abstract A compact high-order finite difference method on unstructured meshes is developed for discretization of the unsteady vorticity transport equations (VTE) for 2-D incompressible flow. The algorithm is based on the Flux Reconstruction Method of Huynh [1, 2], extended to evaluate a Poisson equation for the streamfunction to enforce the kinematic relationship between the velocity and vorticity fields while satisfying the continuity equation. Unlike other finite difference methods for the VTE, where the wall vorticity is approximated by finite differencing the second wall-normal derivative of the streamfunction, the new method applies a Neumann boundary condition for the diffusion of vorticity such that it cancels the slip velocity resulting from the solution of the Poisson equation for the streamfunction. This yields a wall vorticity with order of accuracy consistent with that of the overall solution. In this paper, the high-order VTE solver is formulated and results presented to demonstrate the accuracy and convergence rate of the Poisson solution, as well as the VTE solver using benchmark problems of 2-D flow in lid-driven cavity and backward facing step channel at various Reynolds numbers.

A high order finite-difference scheme for time-domain modelling of time-varying seismoelectric waves

Geophysics ◽  
2021 ◽  
pp. 1-49
Author(s):  
Yanju Ji ◽  
Li Han ◽  
Xingguo Huang ◽  
Xuejiao Zhao ◽  
Kristian Jensen ◽  
...  
Keyword(s):  
Porous Media ◽  
Finite Difference ◽  
Time Domain ◽  
Calculation Method ◽  
High Order ◽  
Difference Approximation ◽  
Time Varying ◽  
Time Step ◽  
Static Approximation ◽  
High Order Finite Difference

Simulation of the seismoelectric effect serves as a useful tool to capture the observed seismoelectric conversion phenomenon in porous media, thus offering promising potential in underground exploration activities to detect pore fluids such as water, oil and gas. The static electromagnetic (EM) approximation is among the most widely used methods for numerical simulation of the seismoelectric responses. However, the static approximation ignores the accompanying electric field generated by the shear wave, resulting in considerable errors when compared to analytical results, particularly under high salinity conditions. To mitigate this problem, we propose a spatial high-order finite-difference time-domain (FDTD) method based on Maxwell's full equations of time-varying EM fields to simulate the seismoelectric response in 2D mode. To improve the computational efficiency influenced by the velocity differences between seismic and electromagnetic waves, different time steps are set according to the stability conditions, and the seismic feedback values of EM time nodes are obtained by linear approximation within the seismic unit time step. To improve the simulation accuracy of the seismoelectric response with the time-varying EM calculation method, finite-difference coefficients are obtained by solving the spatial high-order difference approximation based on Taylor expansion. The proposed method yields consistent simulation results compared to those obtained from the analytical method under different salinity conditions, thus indicating its validity for simulating seismoelectric responses in porous media. We further apply our method to both layered and anomalous body models and extend our algorithm to 3D. Results show that the time-varying EM calculation method could effectively capture the reflection and transmission phenomena of the seismic and EM wavefields at the interfaces of contrasting media. This may allow for the identification of abnormal locations, thus highlighting the capability of seismoelectric response simulation to detect subsurface properties.

Sign in / Sign up

Export Citation Format

Share Document