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 ...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.展开更多
On September 16,2021,a MS6.0 earthquake struck Luxian County,one of the shale gas blocks in the Southeastern Sichuan Basin,China.To understand the seismogenic environment and its mechanism,we inverted a fine three-dim...On September 16,2021,a MS6.0 earthquake struck Luxian County,one of the shale gas blocks in the Southeastern Sichuan Basin,China.To understand the seismogenic environment and its mechanism,we inverted a fine three-dimensional S-wave velocity model from ambient noise tomography using data from a newly deployed dense seismic array around the epicenter,by extracting and jointly inverting the Rayleigh phase and group velocities in the period of 1.6–7.2 s.The results showed that the velocity model varied significantly beneath different geological units.The Yujiasi syncline is characterized by low velocity at depths of~3.0–4.0 km,corresponding to the stable sedimentary layer in the Sichuan Basin.The eastern and western branches of the Huayingshan fault belt generally exhibit high velocities in the NE-SW direction,with a few local low-velocity zones.The Luxian MS6.0 earthquake epicenter is located at the boundary between the high-and low-velocity zones,and the earthquake sequences expand eastward from the epicenter at depths of 3.0–5.0 km.Integrated with the velocity variations around the epicenter,distribution of aftershock sequences,and focal mechanism solution,it is speculated that the seismogenic mechanism of the main shock might be interpreted as the reactivation of pre-existing faults by hydraulic fracturing.展开更多
The Ningdu basin,located in southern Jiangxi province of southwest China,is one of the Mesozoic basin groups which has exploration prospects for geothermal energy.A study on the detailed velocity structure of the Ning...The Ningdu basin,located in southern Jiangxi province of southwest China,is one of the Mesozoic basin groups which has exploration prospects for geothermal energy.A study on the detailed velocity structure of the Ningdu basin can provide important information for geothermal resource exploration.In this study,we deployed a dense seismic array in the Ningdu basin to investigate the 3D velocity structure and discuss implications for geothermal exploration and geological evolution.Based on the dense seismic array including 35 short-period(5 s-100 Hz)seismometers with an average interstation distance of~5 km,Rayleigh surface wave dispersion curves were extracted from the continuous ambient noise data for surface wave tomographic inversion.Group velocity tomography was conducted and the 3D S-wave velocity structure was inverted by the neighborhood algorithm.The results revealed obvious low-velocity anomalies in the center of the basin,consistent with the low-velocity Cretaceous sedimentary rocks.The basement and basin-controlling fault can also be depicted by the S-wave velocity anomalies.The obvious seismic interface is about 2 km depth in the basin center and decreases to 700 m depth near the basin boundary,suggesting spatial thickness variations of the Cretaceous sediment.The fault features of the S-wave velocity profile coincide with the geological cognition of the western boundary basincontrolling fault,which may provide possible upwelling channels for geothermal fluid.This study suggests that seismic tomography with a dense array is an effective method and can play an important role in the detailed investigations of sedimentary basins.展开更多
Accurately picking P-and S-wave arrivals of microseismic(MS)signals in real-time directly influences the early warning of rock mass failure.A common contradiction between accuracy and computation exists in the current...Accurately picking P-and S-wave arrivals of microseismic(MS)signals in real-time directly influences the early warning of rock mass failure.A common contradiction between accuracy and computation exists in the current arrival picking methods.Thus,a real-time arrival picking method of MS signals is constructed based on a convolutional-recurrent neural network(CRNN).This method fully utilizes the advantages of convolutional layers and gated recurrent units(GRU)in extracting short-and long-term features,in order to create a precise and lightweight arrival picking structure.Then,the synthetic signals with field noises are used to evaluate the hyperparameters of the CRNN model and obtain an optimal CRNN model.The actual operation on various devices indicates that compared with the U-Net method,the CRNN method achieves faster arrival picking with less performance consumption.An application of large underground caverns in the Yebatan hydropower station(YBT)project shows that compared with the short-term average/long-term average(STA/LTA),Akaike information criterion(AIC)and U-Net methods,the CRNN method has the highest accuracy within four sampling points,which is 87.44%for P-wave and 91.29%for S-wave,respectively.The sum of mean absolute errors(MAESUM)of the CRNN method is 4.22 sampling points,which is lower than that of the other methods.Among the four methods,the MS sources location calculated based on the CRNN method shows the best consistency with the actual failure,which occurs at the junction of the shaft and the second gallery.Thus,the proposed method can pick up P-and S-arrival accurately and rapidly,providing a reference for rock failure analysis and evaluation in engineering applications.展开更多
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...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.展开更多
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 predi...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.展开更多
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 functio...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.展开更多
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...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.展开更多
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 veloci...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.展开更多
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...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展开更多
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 techni...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.展开更多
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 quadratica...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.展开更多
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 ...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.展开更多
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 necessa...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.展开更多
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 prot...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.展开更多
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...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.展开更多
Joint inversion for crustal velocity structure using only surface wave dispersion and receiver function data often suffers the non-unique solution problem due to lack of P wave velocity constraint in the data.Here we ...Joint inversion for crustal velocity structure using only surface wave dispersion and receiver function data often suffers the non-unique solution problem due to lack of P wave velocity constraint in the data.Here we developed a method for measuring PmP travel time using teleseismic S wave waveform.The major improvement by this method over the previous one-layer-crust search method is its use of a more realistic multi-layer-crust(MLC)velocity model in the study region based on available information.Numerical tests show that compared with the previous search method the MLC search method is faster and more reliable and accurate.One limit of the MLC search method is the requirement of a multi-layer background model,thus this method might not be applicable in regions where the velocity structure is unknown.Nevertheless,our numerical tests show that the MLC search method is robust in the sense that small deviations of the background model from the true velocity model do not influence the results severely.We applied the MLC search method to data from a temporary linear array across the Wabash Valley Seismic Zone in the central USA.We obtained 157 PmP travel-time measurements that show lateral crustal structure variation in the region.The measurements provide additional constraints in joint inversion for crustal velocity structure in this seismically active region.展开更多
Using the seismic records of 83 temporary and 17 permanent broadband seismic stations deployed in Tangshan earthquake region and its adjacent areas(39°N–41.5°N,115.5°E–119.5°E),we conducted a non...Using the seismic records of 83 temporary and 17 permanent broadband seismic stations deployed in Tangshan earthquake region and its adjacent areas(39°N–41.5°N,115.5°E–119.5°E),we conducted a nonlinear joint inversion of receiver functions and surface wave dispersion.We obtained some detailed information about the Tangshan earthquake region and its adjacent areas,including sedimentary thickness,Moho depth,and crustal and upper mantle S-wave velocity.Meanwhile,we also obtained the vP/vS structure along two sections across the Tangshan region.The results show that:(1)the Moho depth ranges from 30 km to 38 km,and it becomes shallower from Yanshan uplift area to North China basin;(2)the thickness of sedimentary layer ranges from 0 km to 3 km,and it thickens from Yanshan uplift region to North China basin;(3)the S-wave velocity structure shows that the velocity distribution of the upper crust has obvious correlation with the surface geological structure,while the velocity characteristics of the middle and lower crust are opposite to that of the upper crust.Compared with the upper crust,the heterogeneity of the middle and lower crust is more obvious;(4)the discontinuity of Moho on the two sides of Tangshan fault suggests that Tangshan fault cut the whole crust,and the low vS and high vP/vS beneath the Tangshan earthquake region may reflect the invasion of mantle thermal material through Tangshan fault.展开更多
Based on the waveform data of 5,076 local earthquakes recorded at 25 stations in Xinjiang during the period from 2009 to 2014 and the observation reports provided by the Xinjiang Digital Seismic Network,a data set of ...Based on the waveform data of 5,076 local earthquakes recorded at 25 stations in Xinjiang during the period from 2009 to 2014 and the observation reports provided by the Xinjiang Digital Seismic Network,a data set of 19,140 attenuation factors t*is obtained by fitting the high-frequency attenuation of S-wave spectra with a genetic algorithm. The spatial distribution of Q_S is determined by inverting the t*data with seismic tomography. The results show that the average Q0 in eastern Tianshan is 520,and there is a significant correlation between the Q_S value distribution or attenuation characteristics it disclosed and the surface structure of the study area. The Q_S value is lower in the intersection area of the mountain basin which is located on the north and south sides of the Tianshan Mountains,and the high Q_S distribution is more concentrated inside the Tianshan orogenic belt. The M≥6. 0 earthquakes have been basically located in the Low-Q_S region since 1900. 24 high heat flow points in eastern Tianshan are located at the north and south of Tianshan Mountains where low Q_S exists,indicating a negative correlation. In addition,there is a positive correlation between the velocity structure and the attenuation structure in the study area,which reflects the consistency of the 2-D attenuation structure with the velocity structure and the two-dimensional density structure.展开更多
The origins of superconductivity and pairing symmetry of order parameters are still controversial problems for FeSe thin films up to date.Under the Neumann boundary conditions,the electromagnetic properties of this sy...The origins of superconductivity and pairing symmetry of order parameters are still controversial problems for FeSe thin films up to date.Under the Neumann boundary conditions,the electromagnetic properties of this system are investigated using the two-band Ginzburg-Landau theory.We calculate the temperature dependence of upper critical field in arbitrary direction and critical supercurrent density through the FeSe film.It is revealed that the normalized upper critical field is independent of the film thickness and all of our theoretical results are in accordance with the experimental data.These thus strongly indicate the existence of two-gap s-wave superconductivity in this material.展开更多
文摘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.
基金This work was supported by the Special Fund of the Institute of Geophysics,China Earthquake Administration(Nos.DQJB22B19,DQJB22R29 and DQJB22B26)the National Natural Science Foundation of China(Nos.41974066,U1839209 and 42074053)。
文摘On September 16,2021,a MS6.0 earthquake struck Luxian County,one of the shale gas blocks in the Southeastern Sichuan Basin,China.To understand the seismogenic environment and its mechanism,we inverted a fine three-dimensional S-wave velocity model from ambient noise tomography using data from a newly deployed dense seismic array around the epicenter,by extracting and jointly inverting the Rayleigh phase and group velocities in the period of 1.6–7.2 s.The results showed that the velocity model varied significantly beneath different geological units.The Yujiasi syncline is characterized by low velocity at depths of~3.0–4.0 km,corresponding to the stable sedimentary layer in the Sichuan Basin.The eastern and western branches of the Huayingshan fault belt generally exhibit high velocities in the NE-SW direction,with a few local low-velocity zones.The Luxian MS6.0 earthquake epicenter is located at the boundary between the high-and low-velocity zones,and the earthquake sequences expand eastward from the epicenter at depths of 3.0–5.0 km.Integrated with the velocity variations around the epicenter,distribution of aftershock sequences,and focal mechanism solution,it is speculated that the seismogenic mechanism of the main shock might be interpreted as the reactivation of pre-existing faults by hydraulic fracturing.
基金supported by China Geological Survey (DD20190083, DD20221662)National Natural Science Foundation of China (41904044, 41974064, 42174076, 41874069)Youth Innovation Promotion Association CAS (2019330).
文摘The Ningdu basin,located in southern Jiangxi province of southwest China,is one of the Mesozoic basin groups which has exploration prospects for geothermal energy.A study on the detailed velocity structure of the Ningdu basin can provide important information for geothermal resource exploration.In this study,we deployed a dense seismic array in the Ningdu basin to investigate the 3D velocity structure and discuss implications for geothermal exploration and geological evolution.Based on the dense seismic array including 35 short-period(5 s-100 Hz)seismometers with an average interstation distance of~5 km,Rayleigh surface wave dispersion curves were extracted from the continuous ambient noise data for surface wave tomographic inversion.Group velocity tomography was conducted and the 3D S-wave velocity structure was inverted by the neighborhood algorithm.The results revealed obvious low-velocity anomalies in the center of the basin,consistent with the low-velocity Cretaceous sedimentary rocks.The basement and basin-controlling fault can also be depicted by the S-wave velocity anomalies.The obvious seismic interface is about 2 km depth in the basin center and decreases to 700 m depth near the basin boundary,suggesting spatial thickness variations of the Cretaceous sediment.The fault features of the S-wave velocity profile coincide with the geological cognition of the western boundary basincontrolling fault,which may provide possible upwelling channels for geothermal fluid.This study suggests that seismic tomography with a dense array is an effective method and can play an important role in the detailed investigations of sedimentary basins.
基金We acknowledge the funding support from National Natural Science Foundation of China(Grant No.42077263).
文摘Accurately picking P-and S-wave arrivals of microseismic(MS)signals in real-time directly influences the early warning of rock mass failure.A common contradiction between accuracy and computation exists in the current arrival picking methods.Thus,a real-time arrival picking method of MS signals is constructed based on a convolutional-recurrent neural network(CRNN).This method fully utilizes the advantages of convolutional layers and gated recurrent units(GRU)in extracting short-and long-term features,in order to create a precise and lightweight arrival picking structure.Then,the synthetic signals with field noises are used to evaluate the hyperparameters of the CRNN model and obtain an optimal CRNN model.The actual operation on various devices indicates that compared with the U-Net method,the CRNN method achieves faster arrival picking with less performance consumption.An application of large underground caverns in the Yebatan hydropower station(YBT)project shows that compared with the short-term average/long-term average(STA/LTA),Akaike information criterion(AIC)and U-Net methods,the CRNN method has the highest accuracy within four sampling points,which is 87.44%for P-wave and 91.29%for S-wave,respectively.The sum of mean absolute errors(MAESUM)of the CRNN method is 4.22 sampling points,which is lower than that of the other methods.Among the four methods,the MS sources location calculated based on the CRNN method shows the best consistency with the actual failure,which occurs at the junction of the shaft and the second gallery.Thus,the proposed method can pick up P-and S-arrival accurately and rapidly,providing a reference for rock failure analysis and evaluation in engineering applications.
基金South China Sea Institute of Oceanology (SCSIO) for providing R/V Shiyan-2 to carry out this experiment,sponsored by Oceanographic Research Vessel Sharing Plan (NORC2016-08) of National Natural Science Foundation of Chinafunded by National Natural Science Foundation of China (Grant Nos. 41776057, 41761134051, 91858213, 41730532 and 91428039)
文摘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.
基金supported by the National Key Basic Research Program of China(973 Program)(Grant No.2011CB013501)the National Science Fund for Distinguished Young Scholars of China(Grant No.51125037)
文摘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.
基金State Key Basic Research Development and Programming Project (G199804070201) State Natural Science Foundation (40074008).
文摘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.
基金supported by the Special Fund for Earthquake Research in the Public Interest(No.201108009)
文摘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.
基金Supported by National Natural Science Foundation of China(No.50378032and No.50538030)Associated Foundation of Earthquake Science(No.201009)Foundation of Heilongjiang Institute of Science and Technology(No.04-15).
文摘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.
文摘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
基金Project supported by Shahid Chamran University of Ahvaz
文摘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.
基金supported by the Research Fellow of Indian School of Mines in Dhanbad (No. 2010DR0016)
文摘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.
基金the Directorate General of Resources for Science Technologythe Higher Education of the Republic of Indonesia for granting a PMDSU scholarship to SR
文摘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.
基金supported by the National Natural Science Foundation of China,Grant No.41774142
文摘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.
基金supported by National Natural Science Foundation of China (No.10947106)
文摘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.
基金a part of author’s M.Sc Research under the project:‘‘Strengthening of Earthquake Engineering Center’’,funded by Higher Education Commission,Government of Pakistan
文摘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.
基金supported by US NSF grants EAR-1249701 and EAR-1661519.
文摘Joint inversion for crustal velocity structure using only surface wave dispersion and receiver function data often suffers the non-unique solution problem due to lack of P wave velocity constraint in the data.Here we developed a method for measuring PmP travel time using teleseismic S wave waveform.The major improvement by this method over the previous one-layer-crust search method is its use of a more realistic multi-layer-crust(MLC)velocity model in the study region based on available information.Numerical tests show that compared with the previous search method the MLC search method is faster and more reliable and accurate.One limit of the MLC search method is the requirement of a multi-layer background model,thus this method might not be applicable in regions where the velocity structure is unknown.Nevertheless,our numerical tests show that the MLC search method is robust in the sense that small deviations of the background model from the true velocity model do not influence the results severely.We applied the MLC search method to data from a temporary linear array across the Wabash Valley Seismic Zone in the central USA.We obtained 157 PmP travel-time measurements that show lateral crustal structure variation in the region.The measurements provide additional constraints in joint inversion for crustal velocity structure in this seismically active region.
基金This research is supported by Spark Program of Earthquake Science(No.XH18065Y)National Natural Science Foundation of China(Nos.41774066 and 41604049)。
文摘Using the seismic records of 83 temporary and 17 permanent broadband seismic stations deployed in Tangshan earthquake region and its adjacent areas(39°N–41.5°N,115.5°E–119.5°E),we conducted a nonlinear joint inversion of receiver functions and surface wave dispersion.We obtained some detailed information about the Tangshan earthquake region and its adjacent areas,including sedimentary thickness,Moho depth,and crustal and upper mantle S-wave velocity.Meanwhile,we also obtained the vP/vS structure along two sections across the Tangshan region.The results show that:(1)the Moho depth ranges from 30 km to 38 km,and it becomes shallower from Yanshan uplift area to North China basin;(2)the thickness of sedimentary layer ranges from 0 km to 3 km,and it thickens from Yanshan uplift region to North China basin;(3)the S-wave velocity structure shows that the velocity distribution of the upper crust has obvious correlation with the surface geological structure,while the velocity characteristics of the middle and lower crust are opposite to that of the upper crust.Compared with the upper crust,the heterogeneity of the middle and lower crust is more obvious;(4)the discontinuity of Moho on the two sides of Tangshan fault suggests that Tangshan fault cut the whole crust,and the low vS and high vP/vS beneath the Tangshan earthquake region may reflect the invasion of mantle thermal material through Tangshan fault.
基金jointly funded by the Contract Oriented Work Task for Seismic Situation in 2017(2017010104)Science for Earthquake Resilience(XH17041Y)Fund of Earthquake Agency of Xinjiang Uygur Autonomous Region(201401)
文摘Based on the waveform data of 5,076 local earthquakes recorded at 25 stations in Xinjiang during the period from 2009 to 2014 and the observation reports provided by the Xinjiang Digital Seismic Network,a data set of 19,140 attenuation factors t*is obtained by fitting the high-frequency attenuation of S-wave spectra with a genetic algorithm. The spatial distribution of Q_S is determined by inverting the t*data with seismic tomography. The results show that the average Q0 in eastern Tianshan is 520,and there is a significant correlation between the Q_S value distribution or attenuation characteristics it disclosed and the surface structure of the study area. The Q_S value is lower in the intersection area of the mountain basin which is located on the north and south sides of the Tianshan Mountains,and the high Q_S distribution is more concentrated inside the Tianshan orogenic belt. The M≥6. 0 earthquakes have been basically located in the Low-Q_S region since 1900. 24 high heat flow points in eastern Tianshan are located at the north and south of Tianshan Mountains where low Q_S exists,indicating a negative correlation. In addition,there is a positive correlation between the velocity structure and the attenuation structure in the study area,which reflects the consistency of the 2-D attenuation structure with the velocity structure and the two-dimensional density structure.
文摘The origins of superconductivity and pairing symmetry of order parameters are still controversial problems for FeSe thin films up to date.Under the Neumann boundary conditions,the electromagnetic properties of this system are investigated using the two-band Ginzburg-Landau theory.We calculate the temperature dependence of upper critical field in arbitrary direction and critical supercurrent density through the FeSe film.It is revealed that the normalized upper critical field is independent of the film thickness and all of our theoretical results are in accordance with the experimental data.These thus strongly indicate the existence of two-gap s-wave superconductivity in this material.