S-wave Component Research Articles

A new decoupling and elastic propagator for efficient elastic reverse time migration

Methods to decompose the elastic wavefield into compressional wave (P-wave) and shear wave (S-wave) components in heterogeneous media without wavefield distortions or energy leakage are the key issues in elastic imaging and inversion. We have introduced a decoupled P- and S-wave propagator to form an efficient elastic reverse time migration (RTM) framework, without assuming homogeneous Lamé parameters. Also, no wave-mode conversions occur using the proposed propagator in the presence of strong heterogeneities, which avoids the potential imaging artifacts caused by wave-mode conversions in the receiver-side backward extrapolation. In the proposed elastic RTM framework, the source-side forward wavefield is simulated with a P-wave propagator. The receiver-side wavefield is back extrapolated with the proposed propagator, using the recorded multicomponent seismic data as input. Compared to the conventional elastic RTM, the proposed framework reduces the computational complexity while preserving the imaging accuracy. We have determined its accuracy and efficiency using two synthetic examples.

Suppression of the s-Wave Order Parameter Near the Surface of the Infinite-Layer Electron-Doped Cuprate Superconductor Sr0.9La0.1CuO2

The temperature dependencies of the in-plane (λab) and out-of-plane (λc) components of the magnetic field penetration depth were investigated near the surface and in the bulk of the electron-doped superconductor Sr0.9La0.1CuO2 by means of magnetization measurements. The measured λab(T) and λc(T) were analyzed in terms of a two-gap model with mixed s+d-wave symmetry of the order parameter. λab(T) is well described by an almost pure anisotropic d-wave symmetry component (≃96%), mainly reflecting the surface properties of the sample. In contrast, λc(T) exhibits a mixed s+d-wave order parameter with a substantial s-wave component of more than 50%. The comparison of λab−2(T) measured near the surface with that determined in the bulk by means of the muon-spin rotation/relaxation technique demonstrates that the suppression of the s-wave component of the order parameter near the surface is associated with a reduction of the superfluid density by more than a factor of two.

Coda-wave decorrelation sensitivity kernels in 2-D elastic media: a numerical approach

SUMMARYProbabilistic sensitivity kernels based on the analytical solution of the diffusion and radiative transfer equations have been used to locate tiny changes detected in late arriving coda waves. These analytical kernels accurately describe the sensitivity of coda waves towards velocity changes located at a large distance from the sensors in the acoustic diffusive regime. They are also valid to describe the acoustic waveform distortions (decorrelations) induced by isotropically scattering perturbations. However, in elastic media, there is no analytical solution that describes the complex propagation of wave energy, including mode conversions, polarizations, etc. Here, we derive sensitivity kernels using numerical simulations of wave propagation in heterogeneous media in the acoustic and elastic regimes. We decompose the wavefield into P- and S-wave components at the perturbation location in order to construct separate P to P, S to S, P to S and S to P scattering sensitivity kernels. This allows us to describe the influence of P- and S-wave scattering perturbations separately. We test our approach using acoustic and elastic numerical simulations where localized scattering perturbations are introduced. We validate the numerical sensitivity kernels by comparing them with analytical kernel predictions and with measurements of coda decorrelations on the synthetic data.

Improvement of the objective function in the velocity structure inversion based on horizontal-to-vertical spectral ratio of earthquake ground motions

SUMMARYGiven the improvements that have been made in the forward calculations of seismic noise horizontal-to-vertical spectral ratios (NHVSRs) or earthquake ground motion HVSRs (EHVSRs), a number of HVSR inversion methods have been proposed to identify underground velocity structures. Compared with the studies on NHVSR inversion, the research on the EHVSR-based inversion methods is relatively rare. In this paper, to make full use of the widely available and constantly accumulating strong-motion observation data, we propose an S-wave HVSR inversion method based on diffuse-field approximation. Herein, the S-wave components of earthquake ground motion recordings are considered as data source. Improvements to the objective function has been achieved in this work. An objective function with the slope term is introduced. The new objective function can mitigate the multisolution phenomenon encountered when working with HVSR curves with multipeaks. Then, a synthetic case is used to show the verification of the proposed method and this method has been applied to invert underground velocity structures for six KiK-net stations based on earthquake observations. The results show that the proposed S-wave EHVSR inversion method is effective for identifying underground velocity structures.

A staggered-grid lowrank finite-difference method for elastic wave extrapolation

Elastic wave extrapolation in the time domain is significant for an elastic wave equation-based processing. To improve the simulation reliability and accuracy of decoupled elastic P- and S- waves, we propose the staggered-grid lowrank finite-difference method based on the elastic wave decomposition. For elastic wave propagation, a lowrank finite-difference method based on the staggered grid is derived to improve the accuracy. Regarding the application of the decoupled elastic wave equation, we derive the finite-difference scheme coefficients which are dependent on velocity. Based on the elastic wave decomposition and plane wave theories, we formulate the elastic wave-extrapolation operators, which contain trigonometric adjustment factors. Accordingly, by applying the lowrank method to approximating the operators, the finite- difference scheme is designed to discretize the decoupled wave equation. The derivation processing implies the combination of elastic wave-mode decomposition and extrapolation. The proposed method enables elastic P- and S-waves to extrapolate in the time-space domain separately and produces accurate P-and S-wave components simultaneously. Dispersion analysis suggests that our proposed method is reliable and accurate in a wide range of wavenumber. Numerical simulation tests on a simple model and the Marmousi2 model validate the accuracy and effectiveness of the method, showing its ability in handling complex structures. Although the operators are accurate only when the medium is homogeneous, they are of high accuracy when the velocity gradient is quite small and are applicable when the velocity gradient is large. The subsequent results of reverse time migration for the Marmousi2 model also suggest that the proposed method is enough to serve as an extrapolator in elastic reverse time migration.

Nuclear spin-lattice relaxation rate in noncentrosymmetric superconductors Y2C3

For a noncentrosymmetric superconductor such as Y2C3, we consider a parity-mixing model composed of spin-singlet s-wave and spin-triplet f-wave pairing components. The d-vector in f-wave state is chosen to be parallel to the Dresselhaus asymmetric spin-orbit coupling vector. It is found that, the quasiparticle excitation spectrum exhibits distinct nodal structure as a consequence of parity-mixing. We demonstrate particularly that such a model can qualitatively account for the existing experimental results of the temperature dependence of the nuclear spin-lattice relaxation rate T1−1 in Y2C3.

D-wave effects in diffractive electroproduction of heavy quarkonia from the photon-like Vrightarrow Q{bar{Q}} transition

We analyse the validity of a commonly used identification between structures of the virtual photon gamma ^*rightarrow Qbar{Q} and vector meson Vrightarrow Q{bar{Q}} transitions. In the existing studies of S-wave vector-meson photoproduction in the literature, such an identification is typically performed in the light-front (LF) frame while the radial component of the meson wave function is rather postulated than computed from the first principles. The massive photon-like Vrightarrow Q{bar{Q}} vertex, besides the S-wave component, also contains an extra D-wave admixture in the Q{bar{Q}} rest frame. However, the relative weight of these contributions cannot be justified by any reasonable nonrelativistic Q{bar{Q}} potential model. In this work, we investigate the relative role of the D-wave contribution starting from the photon-like quarkonium Vrightarrow Q{bar{Q}} transition in both frames: in the Q{bar{Q}} rest frame (with subsequent Melosh spin transform to the LF frame) and in the LF frame (without Melosh transform). In both cases, we employ the same radial wave functions found in the potential approach. We show that the photon-like transition imposed in the Q{bar{Q}} rest frame, with subsequent boosting to the LF frame, leads to significant discrepancies with the experimental data. In the second case, i.e. imposing the photon-like transition straight in the LF frame, we find that the corresponding total J/ psi (1S) photoproduction cross sections are very close to those obtained with the “S-wave only” Vrightarrow Q{bar{Q}} transition, both leading to a good description of the data. However, we find that the “S-wave only” transition leads to a better description of photoproduction data for excited (e.g. psi '(2S)) heavy quarkonium states, which represent a more effective tool for study of D-wave effects. Consequently, the predictions for production of excited states based on the photon-like structure of Vrightarrow Q{bar{Q}} transition should be treated with a great care due to a much stronger sensitivity of the D-wave contribution to the nodal structure of quarkonium wave functions.

Studies of the resonance components in the B_s decays into charmonia plus kaon pair

In this work, the decays of B_s meson to a charmonium state and a K^+K^- pair are carefully investigated in the perturbative QCD approach. Following the latest fit from the LHCb experiment, we restrict ourselves to the case where the produced K^+K^- pair interact in isospin zero S, P, and D wave resonances in the kinematically allowed mass window. Besides the dominant contributions of the phi (1020) resonance in the P-wave and f_2'(1525) in the D-wave, other resonant structures in the high mass region as well as the S-wave components are also included. The invariant mass spectra for most of the resonances in the B_srightarrow J/psi K^+K^- decay are well reproduced. The obtained three-body decay branching ratios can reach the order of 10^{-4}, which seem to be accessible in the near future experiments. The associated polarization fractions of those vector-vector and vector-tensor modes are also predicted, which are compared with the existing data from the LHCb Collaboration.

Dalitz plot analysis of the D+ \u2192 K\u2212K+K+ decay

The resonant structure of the doubly Cabibbo-suppressed decay D+→K−K+K+ is studied for the first time. The measurement is based on a sample of pp-collision data, collected at a centre-of-mass energy of 8 TeV with the LHCb detector and corresponding to an integrated luminosity of 2 fb−1. The amplitude analysis of this decay is performed with the isobar model and a phenomenological model based on an effective chiral Lagrangian. In both models the S-wave component in the K−K+ system is dominant, with a small contribution of the ϕ(1020) meson and a negligible contribution from tensor resonances. The K+K− scattering amplitudes for the considered combinations of spin (0,1) and isospin (0,1) of the two-body system are obtained from the Dalitz plot fit with the phenomenological decay amplitude.

Study and application of PS-wave pre-stack migration in HTI media and an anisotropic correction method

Anisotropy correction is necessary during the processing of converted PSwave seismic data to achieve accurate structural imaging, reservoir prediction, and fracture detection. To effectively eliminate the adverse effects of S-wave splitting and to improve PSwave imaging quality, we tested methods for pre-stack migration imaging and anisotropic correction of PS-wave data. We based this on the propagation rules of seismic waves in a horizontal transverse isotropy medium, which is a fractured medium model that reflects likely subsurface conditions in the field. We used the radial (R) and transverse (T) components of PS-wave data to separate the fast and slow S-wave components, after which their propagation moveout was effectively extracted. Meanwhile, corrections for the energies and propagation moveouts of the R and T components were implemented using mathematical rotation. The PS-wave imaging quality was distinctly improved, and we demonstrated the reliability of our methods through numerical simulations. Applying our methods to three-dimensional and three-component seismic field data from the Xinchang-Hexingchang region of the Western Sichuan Depression in China, we obtained high-quality seismic imaging with continuous reflection wave groups, distinct structural features, and specific stratigraphic contact relationships. This study provides an effective and reliable approach for data processing that will improve the exploration of complex, hidden lithologic gas reservoirs.

Sensitivity analysis of P-waves and S-waves to gas hydrate in the Shenhu area using OBS

Compared to towed streamers, ocean-bottom seismometers (OBS) obtain both S-wave data and richer wavefield information. In this paper, the induced polarization method is used to conduct wavefield separation on OBS data obtained from the Shenhu area in the South China Sea. A comparison of the changes in P- and S-waves, and a comprehensive analysis of geological factors within the area, enable analysis and description of the occurrence of natural gas hydrate in the study area. Results show an increase in P-wave velocity when natural gas hydrate exists in the formation, whereas the S-wave velocity remains almost constant, as S-waves can only propagate through the rock skeleton. Therefore, the bottom-simulating reflection (BSR) response of the P-wave is better than that of the S-wave in the frequency analysis profile. In a wide-angle section, the refractive wave of the hydrate layer is evident when using P-wave components but identification is difficult with S-wave components. This velocity model illustrates the sensitivity of P- and S-wave components to gas hydrate. The use of this polarization method and results of analysis provide technical and theoretical support for research on hydrate deposits and other geological features in the Shenhu area.

Nematic superconducting state in iron pnictide superconductors

Nematic order often breaks the tetragonal symmetry of iron-based superconductors. It arises from regular structural transition or electronic instability in the normal phase. Here, we report the observation of a nematic superconducting state, by measuring the angular dependence of the in-plane and out-of-plane magnetoresistivity of Ba0.5K0.5Fe2As2 single crystals. We find large twofold oscillations in the vicinity of the superconducting transition, when the direction of applied magnetic field is rotated within the basal plane. To avoid the influences from sample geometry or current flow direction, the sample was designed as Corbino-shape for in-plane and mesa-shape for out-of-plane measurements. Theoretical analysis shows that the nematic superconductivity arises from the weak mixture of the quasi-degenerate s-wave and d-wave components of the superconducting condensate, most probably induced by a weak anisotropy of stresses inherent to single crystals.

Compressional- and shear-wave studies of distributed acoustic sensing acquired vertical seismic profile data

Understanding the strengths and limitations of rapidly advancing distributed acoustic sensing (DAS) technology used for recording vertical seismic profile (VSP) data is achieved by comparing DAS and geophone data sets using both compressional-wave (P-wave) and shear-wave (S-wave) VSP data and their corresponding geophysical answer products. We validate the kinematics (time) and dynamics (amplitude) of DAS VSP data by examining the extracted slowness values, response-to-incident angles, corridor stacks, and common-depth-point (CDP) transforms. For kinematics validation, the slowness values computed from P- and S-wave components of DAS VSP data agree with the geophone slowness values. For dynamics validation, we confirm the cos2 θ response of the fiber to the incident angle of the seismic wavefield for P-waves and sin 2θ for S-waves. The amplitudes of the P-wave corridor stacks are comparable; the S-wave corridor stacks are similar for shallow events and differ for later events due to the limited response of the fiber to S-waves at near-vertical angles of incidence. High-quality CDP transform images are obtained for P- and S-waves. These analyses indicate that properly acquired DAS VSP data sets are reliable for the kinetics and dynamics of both P- and S-waves. Once the fiber-optic cable is installed in the well, VSP acquisition costs are greatly reduced because DAS data may be acquired with no additional well intervention. The extensive spatial coverage obtained using fiber-optic cables, the ease of acquiring time-lapse (4D) VSP data, and the reliability of the resulting DAS VSP data sets are making DAS technology an extremely important VSP acquisition tool.

Vector-based excitation amplitude imaging condition for elastic RTM

In recent years, many studies have focused on elastic reverse time migration (RTM). In response to the problems associated with elastic RTM, we propose a new procedure for 2D elastic multicomponent RTM. In this new method, decomposed P- and S-wave components are obtained from the decoupled propagation of the source and receiver wavefields, which allows the expedient calculation of the Poynting vectors and the incident and reflection angles of the P- and S-waves. In addition, we deduce the vector-based excitation amplitude imaging condition. This process automatically accounts for the particle vibration directions when determining the angle-dependent signed reflection coefficients, and does not require the sign to be determined apart from the value of the reflection coefficients. This concept was further extended to the source-normalized crosscorrelation imaging condition. The reflection coefficient of the layered model test was in agreement with the Zoeppritz theory, the PP and PS wave images of the Marmousi II model were clear, and the PS wave images had higher resolution and richer details. In addition, since the calculated reflection coefficients are angle-dependent, they can be easily used for the extraction of angle-domain common-image gathers. Moreover, the imaging condition avoids the polarization reversal in PS wave images and does not require all of the source wavefield data. Consequently, the computation and storage requirements are significantly reduced, which will facilitate the use of the elastic RTM in practice.

An evaluation of numerical approaches for S-wave component simulation in rock blasting

The shear wave (S-wave) component of the total blast vibration always plays an important role in damage to rock or adjacent structures. Numerical approach has been considered as an economical and effective tool in predicting blast vibration. However, S-wave has not yet attracted enough attention in previous numerical simulations. In this paper, three typical numerical models, i.e. the continuum-based elastic model, the continuum-based damage model, and the coupled smooth particle hydrodynamics (SPH)-finite element method (FEM) model, were first introduced and developed to simulate the blasting of a single cylindrical charge. Then, the numerical results from different models were evaluated based on a review on the generation mechanisms of S-wave during blasting. Finally, some suggestions on the selection of numerical approaches for simulating generation of the blast-induced S-wave were put forward. Results indicate that different numerical models produce different results of S-wave. The coupled numerical model was the best, for its outstanding capacity in producing S-wave component. It is suggested that the model that can describe the cracking, sliding or heaving of rock mass, and the movement of fragments near the borehole should be selected preferentially, and priority should be given to the material constitutive law that could record the nonlinear mechanical behavior of rock mass near the borehole.

A novel method for determining the small-strain shear modulus of soil using the bender elements technique

Bender elements technique has become a popular tool for determining shear wave velocity, Vs, hence the small-strain shear modulus of soils, Gmax, thanks to its simplicity and nondestructive character among other advantages. Several methods were proposed to determine the first arrival of Vs. However, none of them can be widely adopted as a standard and there is still an uncertainty on the detection of the first arrival. In this study, bender elements tests were performed on lime-treated soil and both shear wave and compression wave velocities at various frequencies were measured. In-depth analysis showed that the S-wave received signal presents an identical travel time and opposite polarity compared with that of the S-wave components in P-wave received signal, especially at high frequency. From this observation, a novel interpretation method based on the comparison between the S-wave and P-wave received signals at high frequency is proposed. This method enables the determination of the arrival time of the S-wave objectively, avoiding a less reliable first arrival pick-up point. Furthermore, the “π-point” method and cross-correlation method were also employed and the obtained results agree well with those from the proposed method, indicating the accuracy and reliability of the latter. The effects of frequency on the shear wave velocity are also discussed.

Effect of surface charge and s-wave component on (001) surface and bulk states of [formula omitted

The local densities of the (001) surface and bulk states (LDOS) of thin YBa2Cu3O7−δ films were studied based on the self–consistent tight binding model. The CuO-chain and CuO2-plane layers were considered to be a coupling between bands, and then the input surface charge density and the s-wave sub-dominant component were included in the material. In the calculation process, the surface and bulk states were determined, for comparison. It was found that the size of the superconducting gap of a plane layer at the surface increased with an increasing input of the surface charge density, but the LDOS in the bulk were not affected much by this change. When the s-wave component was included in the material, it affected several peaks in the LDOS. This caused some peaks to split while some other peaks became sharper.

Polarized proton-deuteron scattering as a test of time-reversal invariance

Scattering of protons with transversal polarization $p_y^p$ on deuterons with tensor polarization $P_{xz}$ provides a null-test signal for time-reversal (T) invariance violating but parity (P) conserving effects. We calculate the corresponding null-test observable at beam energies 100-1000 MeV within the spin-dependent Glauber theory considering T-violating P-conserving nucleon-nucleon interactions. The S-wave component of the deuteron wave function as well as the D-wave are taken into account and the latter is found to play an important role for the magnitude and the energy dependence of the observable in question. Specifically, with inclusion of the D wave the maximum of the obtained signal is shifted to higher beam energies, i.e. to 700-800 MeV.

Multigrid-based ‘shifted-Laplacian’ preconditioning for the time-harmonic elastic wave equation

We investigate the numerical performance of an iterative solver for a frequency-domain finite-difference discretization of the isotropic elastic wave equation. The solver is based on the 'shifted-Laplacian' preconditioner, originally designed for the acoustic wave equation. This preconditioner represents a discretization of a heavily damped wave equation and can be efficiently inverted by a multigrid iteration. However, the application of multigrid to the elastic case is not straightforward because standard methods, such as point-Jacobi, fail to smooth the S-wave wavenumber components of the error when high P-to-S velocity ratios are present. We consider line smoothers as an alternative and apply local-mode analysis to evaluate the performance of the various components of the multigrid preconditioner. Numerical examples in 2-D demonstrate the efficacy of our method.

Skyrme interaction to second order in nuclear matter

Based on the phenomenological Skyrme interaction various density-dependent nuclear matter quantities are calculated up to second order in many-body perturbation theory. The spin–orbit term as well as two tensor terms contribute at second order to the energy per particle. The simultaneous calculation of the isotropic Fermi-liquid parameters provides a rigorous check through the validity of the Landau relations. It is found that published results for these second order contributions are incorrect in most cases. In particular, interference terms between s-wave and p-wave components of the interaction can contribute only to (isospin or spin) asymmetry energies. Even with nine adjustable parameters, one does not obtain a good description of the empirical nuclear matter saturation curve in the low density region . The reason for this feature is the too strong density-dependence of several second-order contributions. The inclusion of the density-dependent term is therefore indispensable for a realistic description of nuclear matter in the Skyrme framework.