Application of 9-component S-wave 3D seismic data to study sedimentary facies and reservoirs in a biogasbearing area:A case study on the Pleistocene Qigequan Formation in Taidong area,Sanhu Depression,Qaidam Basin,NW China

To solve the problems in restoring sedimentary facies and predicting reservoirs in loose gas-bearing sediment,based on seismic sedimentologic analysis of the first 9-component S-wave 3D seismic dataset of China,a fourth-order isochronous stratigraphic framework was set up and then sedimentary facies and reservoirs in the Pleistocene Qigequan Formation in Taidong area of Qaidam Basin were studied by seismic geomorphology and seismic lithology.The study method and thought are as following.Firstly,techniques of phase rotation,frequency decomposition and fusion,and stratal slicing were applied to the 9-component S-wave seismic data to restore sedimentary facies of major marker beds based on sedimentary models reflected by satellite images.Then,techniques of seismic attribute extraction,principal component analysis,and random fitting were applied to calculate the reservoir thickness and physical parameters of a key sandbody,and the results are satisfactory and confirmed by blind testing wells.Study results reveal that the dominant sedimentary facies in the Qigequan Formation within the study area are delta front and shallow lake.The RGB fused slices indicate that there are two cycles with three sets of underwater distributary channel systems in one period.Among them,sandstones in the distributary channels of middle-low Qigequan Formation are thick and broad with superior physical properties,which are favorable reservoirs.The reservoir permeability is also affected by diagenesis.Distributary channel sandstone reservoirs extend further to the west of Sebei-1 gas field,which provides a basis to expand exploration to the western peripheral area.

Effect of the mixing of s-wave and chiral p-wave pairings on electrical shot noise properties of normal metal/superconductor tunnel junctions

We study theoretically the electrical shot noise properties of tunnel junctions between a normal metal and a superconductor with the mixture of singlet s-wave and chiral triplet p-wave pairing due to broken inversion symmetry. We investigate how the shot noise properties vary as the relative amplitude between the two parity components in the pairing potential is changed. It is demonstrated that some characteristics of the electrical shot noise properties of such tunnel junctions may depend sensitively on the relative amplitude between the two parity components in the pairing potential, and some significant changes may occur in the electrical shot noise properties when the relative amplitude between the two parity components is varied from the singlet s-wave pairing dominated regime to the chiral triplet p-wave pairing dominated regime. In the chiral triplet p-wave pairing dominated regime, the ratio of noise power to electric current is close to 2e both in the in-gap and in the out-gap region. In the singlet s-wave pairing dominated regime, the value of this ratio is close to 4e in the inner gap region but may reduce to about 2e in the outer gap region as the relative amplitude of the chiral triplet pairing component is increased. The variations of the differential shot noise with the bias voltage also exhibit some significantly different features in different regimes. Such different features can serve as useful diagnostic tools for the determination of the relative magnitude of the two parity components in the pairing potential.

Rock critical porosity inversion and S-wave velocity prediction

A critical porosity model is often used to calculate the dry frame elastic modulus by the rock critical porosity value which is affected by many factors. In practice it is hard for us to obtain an accurate critical porosity value and we can generally take only an empirical critical porosity value which often causes errors. In this paper, we propose a method to obtain the rock critical porosity value by inverting P-wave velocity and applying it to predict S-wave velocity. The applications of experiment and log data both show that the critical porosity inversion method can reduce the uncertainty resulting from using an empirical value in the past and provide the accurate critical porosity value for predicting S-wave velocity which significantly improves the prediction accuracy.

The bound weighted average method(BWAM)for predicting S-wave velocity

The shear-wave velocity is a very important parameter in oil and gas seismic exploration, and vital in prestack elastic-parameters inversion and seismic attribute analysis. However, sheafing-velocity logging is seldom carried out because it is expensive. This paper presents a simple method for predicting S-wave velocity which covers the basic factors that influence seismic wave propagation velocity in rocks. The elastic modulus of a rock is expressed here as a weighted arithmetic average between Voigt and Reuss bounds, where the weighting factor, w, is a measurement of the geometric details of the pore space and mineral grains. The S-wave velocity can be estimated from w, which is derived from the P-wave modulus. The method is applied to process well-logging data for a carbonate reservoir in Sichuan Basin, and shows the predicted S-wave velocities agree well with the measured S-wave velocities.

Study of S-wave ray elastic impedance for identifying lithology and fluid

In this paper, we derive an approximation of the SS-wave reflection coefficient and the expression of S-wave ray elastic impedance （SREI） in terms of the ray parameter. The SREI can be expressed by the S-wave incidence angle or P-wave reflection angle, referred to as SREIS and SREIP, respectively. Our study using elastic models derived from real log measurements shows that SREIP has better capability for lithology and fluid discrimination than SREIS and conventional S-wave elastic impedance （SEI）. We evaluate the SREIP feasibility using 25 groups of samples from Castagna and Smith （1994）. Each sample group is constructed by using shale, brine-sand, and gas-sand. Theoretical evaluation also indicates that SRE1P at large incident angles is more sensitive to fluid than conventional fluid indicators. Real seismic data application also shows that SRE1P at large angles calculated using P-wave and S-wave impedance can efficiently characterize tight gas-sand.

Crustal S-wave velocity structure across the northeastern South China Sea continental margin: implications for lithology and mantle exhumation

The northeastern margin of the South China Sea (SCS), developed from continental rifting and breakup, is usually thought of as a non-volcanic margin. However, post-spreading volcanism is massive and lower crustal high-velocity anomalies are widespread, which complicate the nature of the margin here. To better understand crustal seismic velocities, lithology, and geophysical properties, we present an S-wave velocity (VS) model and a VP/VS model for the northeastern margin by using an existing P-wave velocity (VP) model as the starting model for 2-D kinematic S-wave forward ray tracing. The Mesozoic sedimentary sequence has lower VP/VS ratios than the Cenozoic sequence;in between is a main interface of P-S conversion. Two isolated high-velocity zones (HVZ) are found in the lower crust of the continental slope, showing S-wave velocities of 4.0–4.2 km/s and VP/VS ratios of 1.73–1.78. These values indicate a mafic composition, most likely of amphibolite facies. Also, a VP/VS versus VP plot indicates a magnesium-rich gabbro facies from post-spreading mantle melting at temperatures higher than normal. A third high-velocity zone (VP : 7.0–7.8 km/s;VP/VS: 1.85–1.96), 70-km wide and 4-km thick in the continent-ocean transition zone, is most likely to be a consequence of serpentinization of upwelled upper mantle. Seismic velocity structures and also gravity anomalies indicate that mantle upwelling/ serpentinization could be the most severe in the northeasternmost continent-ocean boundary of the SCS. Empirical relationships between seismic velocity and degree of serpentinization suggest that serpentinite content decreases with depth, from 43% in the lower crust to 37% into the mantle.

Pure S-waves in land P-wave source VSP data

Conventional land vertical seismic profiling （VSP） exploration usually uses P-wave sources and three-component geophones for receivers, emphasizing P- and converted S-waves. Previous studies show that both dynamite borehole shots and vertical vibrations from controllable seismic sources at the surface will produce relatively strong pure P-waves and weaker pure S-waves. Interfaces with a large Poisson＇s ratio difference have a positive influence on the formation of strong transmitted converted S-waves. By a comparative analysis of pure S-waves from sources and converted downgoing S-waves, we believe that the main frequency of pure S-waves is usually lower than pure P-waves while the main frequency of downgoing converted S-waves is close to that of P-waves. We have studied zero-offset and offset VSP data from land P-wave sources. Results show that pure S-waves commonly exist in these data with differences in wave intensity. S-wave velocity can be obtained from the P-wave source zero-offset VSP data. Finally, we discuss the bright future of joint application of VSP P-and S-waves and the full use of S-waves in P-wave source VSP data.

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.

Mapping crustal S-wave velocity structure with SV-component receiver function method

In this article, we analyze the characters of SV-component receiver function of teleseismic body waves and its advantages in mapping the S-wave velocity structure of crust in detail. Similar to radial receiver function, SV-component receiver function can be obtained by directly deconvolving the P-component from the SV-component of teleseismic recordings. Our analyses indicate that the change of amplitude of SV-component receiver function against the change of epicentral distance is less than that of radial receiver function. Moreover, the waveform of SV-component receiver function is simpler than the radial receiver function and gives prominence to the PS converted phases that are the most sensitive to the shear wave velocity structure in the inversion. The synthetic tests show that the convergence of SV-component receiver function inversion is faster than that of the radial receiver function inversion. As an example, we investigate the S-wave velocity structure beneath HIA sta-tion by using the SV-component receiver function inversion method.

Estimation of S-wave Velocity for Gas Hydrate Reservoir in the Shenhu Area,North South China Sea

Estimation of S-wave velocity using logging data has mainly been performed for sandstone, mudstone and oil and gas strata, while its application to hydrate reservoirs has been largely overlooked. In this paper we present petxophysical methods to estimate the S-wave velocity of hydrate reservoirs with the P-wave velocity and the density as constraints. The three models used in this paper are an equivalent model （MBGL）, a three-phase model （TPBE）, and a thermo-elasticity model （TEM）. The MBGL model can effectively describe the internal relationship among the components of the rock, and the estimated P-wave velocities are in good agreement with the measured data （2.8% error）. However, in the TPBE model, the solid, liquid and gas phases axe considered to be independent of each other, and the estimation results are relatively low （46.6% error）. The TEM model is based on the sensitivity of the gas hydrate to temperature and pressure, and the accuracy of the estimation results is also high （3.6% error）. Before the estimation, the occurrence patterns of hydrates in the Shenhu area were examined, and occurrence state one （the hydrate is in solid form in the reservoir） was selected for analysis. By using the known P-wave velocity and density as constraints, a reasonable S-wave velocity value （ranging from 400 to 1100 m s 1 and for a hydrate layer of 1100 m s 1） can be obtained through multiple iterations. These methods and results provide new data and technical support for further research on hydrates and other geological features in the Shenhu area.

Estimation of Shallow S-Wave Velocity Structure of Two Practical Sites from Microtremors Array Observation in Tangshan Area

Microtremors array observation for estimating S-wave velocity structure from phase velocities of Rayleigh and Love wave on two practical sites in Tangshan area by a China-US joint group are researched.The phase velocities of Rayleigh wave are estimated from vertical component records and those of Love wave are estimated from three-component records of microtremors array using modified spatial auto-correlation method.Haskell matrix method is used in calculating Rayleigh and Love wave phase velocities,and the shallow S-wave velocity structure of two practical sites are estimated by means of a hybrid approach of Genetic Algorithm and Simplex.The results are compared with the PS logging data of the two sites,showing it is feasible to estimate the shallow S-wave velocity structure of practical site from the observation of microtremor array.

The strain seismograms of P- and S-waves of a local event recorded by four-gauge borehole strainmeter

At a sampling rate of 100 samples per second,the YRY-4 four-gauge borehole strainmeters（FGBS） are capable of recording transient strains caused by seismic waves such as P and S waves or strain seismograms. At such a high sampling rate, data from the YRY-4 strainmeters demonstrate fairly satisfactory self-consistency. The strain tensor seismograms demonstrate the senses of motion of P waves, that is, the type of seismic wave travels in the direction of the maximum normal strain change. The observed strain patterns of S waves significantly differ from those of P waves and should contain information about the source mechanism. Spectrum analysis shows that the strain seismograms are consistent with conventional broadband seismograms from the same site.

Interpretation of S-Wave Data from Fanshi-Taipusiqi DSS Profile and Analysis of Correlation between Deep Structural Characteristics and Seismicity

By processing S-wave data from the Fanshi-Huai’an-Taipusiqi DSS profile,which is a three-component,wide-angle reflection/refraction profile,and in the light of the results from P-wave interpretation,two-dimensional(2-D)structures of the crust and upper mantle are presented,including S-wave velocity Vs and the physical parameter of medium-Poisson’s ratio a.Taking other geological and geophysical information into account,and with reference to the results from petrophysical experiments at home and abroad,we carried out interpretation and inference with respect to deep crustal structure,tectonics,and lithologic characters.It has been concluded that in the upper and middle crust,a values are mostly not greater than 0.25,and rocks,which generally assume brittle,are mainly composed of granite; the rocks in the lower layer of the upper crust between Yangyuan-Huai’an containing inorganic CO2 itself releases carbon; for the rocks in the lower crust and crust-mantle transitional zone,which are comparatively

Validation of the Wiedemann–Franz law in a granular s-wave superconductor in the nanometer scale

The present study tries to evaluate the validity of the Wiedemann-Franz law in a granular s-wave superconductor in the presence of concentrated impurities. By using Green＇s function method and the Kubo formula technique, three distinct contributions of the Aslamazov-Larkin, the Maid-Thompson and, the density of states are calculated for both the electrical conductivity and the thermal conductivity in a granular s-wave superconductor. It is demonstrated that these different con- tributions to the fluctuation conductivity depend differently on the tunneling because of their different natures. This study examines the transport in a granular superconductor system in three dimensions in the limit of large tunneling conductance, which makes it possible to ignore all localization effects and the Coulomb interaction. We find that the tunneling is efficient near the critical temperature and that there is a crossover to the characteristic behavior of a homogeneous system. When it is far from the critical temperature, the tunneling is not effective and the system behaves as an ensemble of real zero-dimensional grains. The results show that the Wiedemann-Franz law is violated in both temperature regions.

Propagation of P-and S-waves in initially stressed gravitating half space

The present paper contributes in studying the phase velocities of P- and S-waves in a half space subjected to a compressive initial stress and gravity field. The density and acceleration due to gravity vary quadratically along the depth. The dispersion equation is derived in a closed form. It is shown that the phase velocities depend not only on the initial stress, gravity, and direction of propagation but also on the inhomogeneity parameter associated with the density and acceleration due to gravity. Various particular cases are obtained, and the results match with the classical results. Numerical investigations on the phase velocities of P- and S-waves against the wave number are made for various sets of values of the material parameters, and the results are illustrated graphically. The graphical user interface model is developed to generalize the effect.

Double-difference tomography of P- and S-wave velocity structure beneath the western part of Java, Indonesia

West Java in the western part of the Sunda Arc has a relatively high seismicity due to subduction activity and faults.In this study,double-difference tomography was used to obtain the 3D velocity tomograms of P and S waves beneath the western part of Java.To infer the geometry of the structure beneath the study area,precise earthquake hypo・center determination was first performed before tomographic imaging.For this,earthquake waveform data were extracted from the regional Meteorological,Climatological,Geophysical Agency(BMKG)network of Indonesia from South Sumatra to Central Java.The P and S arrival times for about 1,000 events in the period April 2009 to July 2016 were selected,the key features being events of magnitude>3,azimuthal gap<210°and number of phases>8.A nonlinear method using the oct-tree sampling algorithm from the NonLinLoc program was employed to determine the earthquake hypocenters.The hypocenter locations were then relocated using double-difference tomography(tomoDD).A significant reduction of travel-time(root mean square basis)and a better clustering of earthquakes were achieved which correlated well with the geological structure in West Java.Double-difference tomography was found to give a clear velocity structure,especially beneath the volcanic arc area,i.e.,under Mt Anak Krakatau,Mt Salak and the mountains complex in the southern part of West Java.Low velocity anomalies for the P and S waves as well as the vp/vs ratio below the volcanoes indicated possible partial melting of the upper mantle which ascended from the subducted slab beneath the volcanic arc.

Non-artifact vector P- and S-wave separation for elastic reverse time migration

Elastic reverse time migration(RTM)uses the elastic wave equation to extrapolate multicomponent seismic data to the subsurface and separate the elastic wavefield into P-and S-waves.P-and S-wave separation is a necessary step in elastic RTM to avoid crosstalk between coupled wavefields.However,the current curl-divergence operator-based separation method has a polarity reversal problem in PS imaging,and vector separation methods often have separation artifacts at the interface,which affects the quality of the imaging stack.We propose a non-artifact P-and S-wave separation method based on the first-order velocity-strain equation.This equation is used for wavefield extrapolation and separation in the first-order staggered-grid finite-difference scheme,and the storage and calculation amounts are consistent with the classical first-order velocity-stress equation.The separation equation does not calculate the partial derivatives of the elastic parameters,and thus,there is no artifact in the separated Pand S-waves.During wavefield extrapolation,the dynamic characteristics of the reflected wave undergo some changes,but the transmitted wavefield is accurate;therefore,it does not affect the dynamic characteristics of the final migration imaging.Through numerical examples of 2 D simple models,part SEAM model,BP model,and 3 D 4-layer model,different wavefield separation methods and corresponding elastic RTM imaging results are analyzed.We found that the velocity-strain based elastic RTM can image subsurface structures well,without spike artifacts caused by separation artifacts,and without polarity reversal phenomenon of the PS imaging.

Localization of s-Wave and Quantum Effective Potential of a Quasi-free Particle with Position-Dependent Mass

The properties of the 8-wave for a quasl-free partide with position-dependent mass （PDM） have been discussed in details. Differed from the system with constant mass in which the localization of the s-wave for the free quantum particle around the origin only occurs in two dimensions, the quasi-free particle with PDM can experience attractive forces in D dimensions except D = 1 when its mass function satisfies some conditions. The effective mass of a particle varying with its position can induce effective interaction, which may be attractive in some cases. The analytical expressions of the eigenfunctions and the corresponding probability densities for the 8-waves of the two- and three-dimensional systems with a special PDM are given, and the existences of localization around the origin for these systems are shown.

Crustal attenuation characteristics of S-waves beneath the Eastern Tohoku region,Japan

An inversion method was applied to crustal earthquakes dataset to find S-wave attenuation characteristics beneath the Eastern Tohoku region of Japan. Accelerograms from 85 shallow crustal earthquakes up to 25 km depth and magnitude range between 3.5 and 5.5 were analyzed to estimate the seismic quality factor Qs. A homogeneous attenuation model Qs for the wave propagation path was evaluated from spectral amplitudes, at 24 different frequencies between 0.5 and 20 Hz by using generalized inversion technique. To do this, non-parametric attenuation functions were calculated to observe spectral amplitude decay with hypocentral distance. Then, these functions were parameterized to estimate Qs. It was found that in Eastern Tohoku region, the Qs frequency dependence can be approximated with the function 33 f 1.22 within a frequency range between 0.5 and 20 Hz. However, the frequency dependence of Qs in the frequency range between 0.5 and 6 Hz is best approximated by Qs （f） = 36 f 0.94 showing relatively weaker frequency dependence as compared to the relation Qs （f） = 6 f^ 2.09 for the frequency range between 6 and 15 Hz. These results could be used to estimate source and site parameters for seismic hazard assessment in the region.

Scaling Law of s-Wave Valence Proton Distributions

Inthe framework of the single-particle potential model, the root-meamsquare radii of s-wave valence proton distributions for P and S isotopes have been systematically calculated. The scaling law of s-wave valence proton distributions is obtained by analyzing the relation between radii and separation energies. Necessary condition for the occurrence of 2sl/2-state proton halos is deduced from the scaling law. Quantitative condition derived from this work can to some extent serve as reference for proton halo search.