scholarly journals On accuracy of the finite-difference and finite-element schemes with respect to P-wave to S-wave speed ratio

Author(s):  
Peter Moczo ◽  
Jozef Kristek ◽  
Martin Galis ◽  
Peter Pazak
2009 ◽  
Vol 45 (3) ◽  
pp. 1614-1617 ◽  
Author(s):  
B. Van de Wiele ◽  
A. Manzin ◽  
L. Dupre ◽  
F. Olyslager ◽  
O. Bottauscio ◽  
...  

1974 ◽  
Vol 64 (6) ◽  
pp. 1789-1808 ◽  
Author(s):  
Robert Burridge ◽  
Cesar Levy

abstract It has recently been shown (Burridge, 1973) that in two dimensions plane-strain shear cracks lacking cohesion may run at speeds ranging from the Rayleigh-wave to the S-wave speed or possibly at the P-wave speed. On the other hand, it has long been known that in antiplane strain, cracks lacking cohesion must run at least at the S-wave speed. Since locally at the edge of a three-dimensional crack there is a combination of plane and antiplane strain, we have sought and found solutions for circular shear cracks expanding at the S-wave speed and at the P-wave speed. These have finite shear tractions ahead of the crack and so may correspond to frictional sliding in the absence of cohesion. The analysis combines the method of Kostrov (1964b) with the results of Burridge (1973). We carry out a complete evaluation for the displacements and tractions on the fault plane, and far-field radiation for the S-wave-speed crack. The corresponding evaluations for the P-wave speed are not elementary and are not attempted here. As far as the authors are aware, these are the first analytic solutions of three-dimensional crack problems which satisfy a physically plausible fracture criterion for failure under shear.


Geophysics ◽  
2006 ◽  
Vol 71 (5) ◽  
pp. C57-C67 ◽  
Author(s):  
Xian-Yun Wu ◽  
Ru-Shan Wu

We propose a dual-domain, one-way, elastic thin-slab method for fast and accurate amplitude variation with offset (AVO) modeling. In this method, the wavefield propagates in the wavenumber domain and interacts with heterogeneity in the space domain. The approach requires much less memory and is two to three orders of magnitude faster than a full-wave method using finite difference or finite element. The thin-bed AVO and AVOs with lateral parameter variations have been conducted using the thin-slab method and compared with reflectivity and finite-difference methods, respectively. It is shown that the thin-slab method can be used to accurately model reflections for most sedimentary rocks that have intermediate parameter perturbations ([Formula: see text] for P-wave velocity and [Formula: see text] for S-wave velocity). The combined effects of overburden structure and the scattering associated with heterogeneities on AVO have been investigated using the thin-slab method. Properties of the target zone and overburden structure control the AVO trends at overall offsets. Scattering associated with heterogeneities increases local variance in the reflected amplitudes and becomes significant for the sedimentary models with weak reflections. Interpretation of AVO observations based on homogeneous elastic models would therefore bias the estimated properties of the target. Furthermore, these effects can produce different apparent AVO trends in different offset ranges.


2021 ◽  
Author(s):  
Y.J. Zhou ◽  
Hongtao Wang

Abstract This study solidifies the aeolian sand by microbial induced carbonate precipitation (MICP) technique. The effects of cementation solution with different concentrations, particle size, and grouting batches are examined via the bender element, unconfined compressive test and scanning electron microscope (SEM). The bender element results show that the wave speed of loose aeolian sand is 200m/s; however, after solidify the aeolian sand, the speed of P-wave is about 450-600m/s and S-wave is about 350-500m/s. Additionally, the unconfined compressive strength (UCS) results indicate that when the concentration of cementation solution is 0.75mol/L, the UCS of bio-solidified sand sample is the highest. Then, compared with the aeolian sand with original grade, the particles ranging from 0.1-0.4mm have a better cementation effect. Moreover, the UCS of bio-solidified sand samples increases along with the grouting batch. From the SEM images, it can be seen that when the grouting batch reaches to five times, the particles are almost completely covered by CaCO 3 crystals compared with the three batches and four batches.


Geophysics ◽  
2009 ◽  
Vol 74 (5) ◽  
pp. T75-T95 ◽  
Author(s):  
Stéphane Operto ◽  
Jean Virieux ◽  
A. Ribodetti ◽  
J. E. Anderson

A 2D finite-difference, frequency-domain method was developed for modeling viscoacoustic seismic waves in transversely isotropic media with a tilted symmetry axis. The medium is parameterized by the P-wave velocity on the symmetry axis, the density, the attenuation factor, Thomsen’s anisotropic parameters [Formula: see text] and [Formula: see text], and the tilt angle. The finite-difference discretization relies on a parsimonious mixed-grid approach that designs accurate yet spatially compact stencils. The system of linear equations resulting from discretizing the time-harmonic wave equation is solved with a parallel direct solver that computes monochromatic wavefields efficiently for many sources. Dispersion analysis shows that four grid points per P-wavelength provide sufficiently accurate solutions in homogeneous media. The absorbing boundary conditions are perfectly matched layers (PMLs). The kinematic and dynamic accuracy of the method wasassessed with several synthetic examples which illustrate the propagation of S-waves excited at the source or at seismic discontinuities when [Formula: see text]. In frequency-domain modeling with absorbing boundary conditions, the unstable S-wave mode is not excited when [Formula: see text], allowing stable simulations of the P-wave mode for such anisotropic media. Some S-wave instabilities are seen in the PMLs when the symmetry axis is tilted and [Formula: see text]. These instabilities are consistent with previous theoretical analyses of PMLs in anisotropic media; they are removed if the grid interval is matched to the P-wavelength that leads to dispersive S-waves. Comparisons between seismograms computed with the frequency-domain acoustic TTI method and a finite-difference, time-domain method for the vertical transversely isotropic elastic equation show good agreement for weak to moderate anisotropy. This suggests the method can be used as a forward problem for viscoacoustic anisotropic full-waveform inversion.


Geophysics ◽  
2006 ◽  
Vol 71 (5) ◽  
pp. S199-S207 ◽  
Author(s):  
Robert Sun ◽  
George A. McMechan ◽  
Chen-Shao Lee ◽  
Jinder Chow ◽  
Chen-Hong Chen

Using two independent, 3D scalar reverse-time depth migrations, we migrate the reflected P- and S-waves in a prestack 3D, three-component (3-C), elastic seismic data volume generated with a P-wave source in a 3D model and recorded at the top of the model. Reflected P- and S-waves are extracted by divergence (a scalar) and curl (a 3-C vector) calculations, respectively, during shallow downward extrapolation of the elastic seismic data. The imaging time for the migrations of both the reflected P- and P-S converted waves at each point is the one-way P-wave traveltime from the source to that point.The divergence (the extracted P-waves) is reverse-time extrapolated using a finite-difference solution of the 3D scalar wave equation in a 3D P-velocity modeland is imaged to obtain the migrated P-image. The curl (the extracted S-waves) is first converted into a scalar S-wavefield by taking the curl’s absolute value as the absolute value of the scalar S-wavefield and assigning a positive sign if the curl is counterclockwise relative to the source or a negative sign otherwise. This scalar S-wavefield is then reverse-time extrapolated using a finite-difference solution of the 3D scalar wave equation in a 3D S-velocity model, and it is imaged with the same one-way P-wave traveltime imaging condition as that used for the P-wave. This achieves S-wave polarity uniformity and ensures constructive S-wave interference between data from adjacent sources. The algorithm gives satisfactory results on synthetic examples for 3D laterally inhomogeneous models.


2018 ◽  
Vol 6 (3) ◽  
pp. T569-T584 ◽  
Author(s):  
Bob A. Hardage ◽  
Donald Wagner

We use finite-difference modeling to illustrate direct-P and direct-SV radiation produced by land-based, vertical-displacement sources such as vertical vibrators and vertical impacts, sources that reflection seismologists have used for decades to image subsurface targets. We are particularly interested in understanding the nature of SV radiation that these sources create. Such insight will be important for those in the seismic interpretation community who are interested in S-wave reflection seismology. One distinction between our work and that of earlier investigators is that we go beyond traditional assumptions that a propagation medium has a flat surface and is isotropic and homogeneous throughout the model space. Instead, we use models composed of small finite-difference cells that allow earth anomalies with dimensions of only a few feet to be near a source station. These small irregularities are expressed as variations in propagation velocities in a small population of model cells, or as slight changes in surface elevation. These small, realistic, near-source irregularities cause robust direct-SV illumination to travel in near-vertical takeoff angles from a surface-based source station and to often travel in true-vertical takeoff angles. We also depart from the conventional practice of calculating direct-P and direct-SV radiation produced by only a single vertical vibrator. Instead, we model direct-P and direct-SV illumination produced by arrays of 2, 3, and 4 inline vertical vibrators. We find that vibrator arrays create SV-mode beam forming that causes a rich amount of direct-SV radiation to occur inside the critical takeoff-angle cone below a source station that extends from [Formula: see text] to [Formula: see text] from vertical. We consider our source-radiation models to be more realistic representations of conditions encountered in actual seismic field practice than previous models that use oversimplified earth conditions. Our purpose is to provide interpreters information that will let them decide if common, surface-based, land sources such as vertical vibrators can be used to illuminate geologic targets with direct-P and direct-SV illuminating wavefields rather than limiting the use of these sources to P-wave imaging only, as has been standard practice for several decades. Our results indicate that vertical vibrators and vertical impacts produce illuminating wavefields that are appropriate for imaging deep geologic targets. Our models apply to the total family of surface-based, vertical-displacement sources, which includes vertical vibrators and accelerated weight drops.


2010 ◽  
Vol 20 (02) ◽  
pp. 265-295 ◽  
Author(s):  
JÉRÔME DRONIOU ◽  
ROBERT EYMARD ◽  
THIERRY GALLOUËT ◽  
RAPHAÈLE HERBIN

We investigate the connections between several recent methods for the discretization of anisotropic heterogeneous diffusion operators on general grids. We prove that the Mimetic Finite Difference scheme, the Hybrid Finite Volume scheme and the Mixed Finite Volume scheme are in fact identical up to some slight generalizations. As a consequence, some of the mathematical results obtained for each of the methods (such as convergence properties or error estimates) may be extended to the unified common framework. We then focus on the relationships between this unified method and nonconforming Finite Element schemes or Mixed Finite Element schemes. We also show that for isotropic operators, on particular meshes such as triangular meshes with acute angles, the unified method boils down to the well-known efficient two-point flux Finite Volume scheme.


Sign in / Sign up

Export Citation Format

Share Document