A typical blasting vibration wave is a composite wave,and its attenuation law is affected by the type of dominant wave component.The purpose of the present study is to establish an attenuation equation of the peak par...A typical blasting vibration wave is a composite wave,and its attenuation law is affected by the type of dominant wave component.The purpose of the present study is to establish an attenuation equation of the peak particle velocity(PPV),taking into account the attenuation characteristics of P-,S-and R-waves in the blasting vibration wave.Field blasting tests were carried out as a case to specifically apply the proposed equation.In view of the fact that the discrete properties of rock mass will inevitably cause the uncertainty of blasting vibration,we also carried out a probability analysis of PPV uncertainty,and introduced the concept of reliability to evaluate blasting vibration.The results showed that the established attenuation equation had a higher prediction accuracy,and can be considered as a promising equation implemented on more complex sites.The adopted uncertainty analysis method can comprehensively take account of the attenuation law of blasting vibration measured on site and discrete properties of rock masses.The obtained distribution of the PPV uncertainty factor can quantitatively evaluate the reliability of blasting vibration,which is a powerful and necessary supplement to the PPV attenuation equation.展开更多
The recently built China Digital Seismic Network consists of the China National Digital Seismic Network (CNDSN),31 regional seismic networks and several small aperture arrays with more than 1 000 stations including 85...The recently built China Digital Seismic Network consists of the China National Digital Seismic Network (CNDSN),31 regional seismic networks and several small aperture arrays with more than 1 000 stations including 850+ broadband stations.It forms a gigantic seismic array that provides an unprecedented opportunity to study the Earth's deep interior besides its routine task of seismic monitoring.Many modern seismic studies rely on rotation of vertical and horizontal components in order to separate different types of seismic waves.Knowledge of the orientations of the two horizontal components thus is important to perform a correction rotation.We analyzed particle motions of teleseismic P waves recorded by the network and used them to estimate the northcomponent azimuth of each station.An SNR-weighted-multi-event method was introduced to obtain component azimuths that best explain the P-wave particle motions of all the events recorded at a station.The method provides robust estimates including a measurement error calculated from background noise levels.We found that about one third of the stations have some sort of problems,including misorientation of the two horizontal components,mislabeling and polarity reversal in one or more components.These problems need to be taken into account for any rotation based seismic studies.展开更多
Ultrasonic Lamb waves are considered as a sensitive and effective tool for nondestructive testing and evaluation of plate-like or pipe-like structures. The nature of multimode and dispersion causes the wave packets to...Ultrasonic Lamb waves are considered as a sensitive and effective tool for nondestructive testing and evaluation of plate-like or pipe-like structures. The nature of multimode and dispersion causes the wave packets to spread, and the modes overlap in both time and frequency domains as they propagate through the structures. By using a two-component laser interferometer technique, in combination with a priori knowledge of the dispersion characteristics and wave structure information of Lamb wave modes, a two-component signal processing technique is presented for implementing dispersion removal and mode separation simultaneously for two modes mixture signals of Lamb waves. The proposed algorithm is first processed and verified using synthetic Lamb wave signals. Then, the two-component displacements test experiment is conducted using different aluminum plate samples. Moreover, we confirm the effectiveness and robustness of this method.展开更多
The dispersion relation and electromagnetic polarization of the plasma waves are comprehensively studied in cold electron,proton,and heavy charged particle plasmas.Three modes are classified as the fast,intermediate,a...The dispersion relation and electromagnetic polarization of the plasma waves are comprehensively studied in cold electron,proton,and heavy charged particle plasmas.Three modes are classified as the fast,intermediate,and slow mode waves according to different phase velocities.When plasmas contain positively-charged particles,the fast and intermediate modes can interact at the small propagating angles,whereas the two modes are separate at the large propagating angles.The near-parallel intermediate and slow waves experience the linear polarization,circular polarization,and linear polarization again,with the increasing wave number.The wave number regime corresponding to the above circular polarization shrinks as the propagating angle increases.Moreover,the fast and intermediate modes cause the reverse change of the electromagnetic polarization at the special wave number.While the heavy particles carry the negative charges,the dispersion relations of the fast and intermediate modes are always separate,being independent of the propagating angles.Furthermore,this study gives new expressions of the three resonance frequencies corresponding to the highly-oblique propagation waves in the general three-component plasmas,and shows the dependence of the resonance frequencies on the propagating angle,the concentration of the heavy particle,and the mass ratio among different kinds of particles.展开更多
The nonlinear Schr?dinger(NLS) equation and Boussinesq equation are two very important integrable equations.They have widely physical applications. In this paper, we investigate a nonlinear system, which is the two-co...The nonlinear Schr?dinger(NLS) equation and Boussinesq equation are two very important integrable equations.They have widely physical applications. In this paper, we investigate a nonlinear system, which is the two-component NLS equation coupled to the Boussinesq equation. We obtain the bright–bright, bright–dark, and dark–dark soliton solutions to the nonlinear system. We discuss the collision between two solitons. We observe that the collision of bright–bright soliton is inelastic and two solitons oscillating periodically can happen in the two parallel-traveling bright–bright or bright–dark soliton solution. The general breather and rogue wave solutions are also given. Our results show again that there are more abundant dynamical properties for multi-component nonlinear systems.展开更多
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.展开更多
We study the generalized Darboux transformation to the three-component coupled nonlinear Schr ¨odinger equation.First-and second-order localized waves are obtained by this technique.In first-order localized wave,...We study the generalized Darboux transformation to the three-component coupled nonlinear Schr ¨odinger equation.First-and second-order localized waves are obtained by this technique.In first-order localized wave,we get the interactional solutions between first-order rogue wave and one-dark,one-bright soliton respectively.Meanwhile,the interactional solutions between one-breather and first-order rogue wave are also given.In second-order localized wave,one-dark-one-bright soliton together with second-order rogue wave is presented in the first component,and two-bright soliton together with second-order rogue wave are gained respectively in the other two components.Besides,we observe second-order rogue wave together with one-breather in three components.Moreover,by increasing the absolute values of two free parameters,the nonlinear waves merge with each other distinctly.These results further reveal the interesting dynamic structures of localized waves in the three-component coupled system.展开更多
Many different effective reflection information are often contaminated by exterior and random noise which concealed in the seismic data.Traditional single or fixed transform is not suit for exploiting their complicate...Many different effective reflection information are often contaminated by exterior and random noise which concealed in the seismic data.Traditional single or fixed transform is not suit for exploiting their complicated characteristics and attenuating the noise.Recent years,a novel method so-called morphological component analysis(MCA)is put forward to separate different geometrical components by amalgamating several irrelevance transforms.According to study the local singular and smooth linear components characteristics of seismic data,we propose a method of suppressing noise by integrating with the advantages of adaptive K-singular value decomposition(K-SVD)and wave atom dictionaries to depict the morphological features diversity of seismic signals.Numerical results indicate that our method can dramatically suppress the undesired noises,preserve the information of geologic body and geological structure and improve the signal-to-noise ratio of the data.We also demonstrate the superior performance of this approach by comparing with other novel dictionaries such as discrete cosine transform(DCT),undecimated discrete wavelet transform(UDWT),or curvelet transform,etc.This algorithm provides new ideas for data processing to advance quality and signal-to-noise ratio of seismic data.展开更多
The 10 920 stress indicators collected so far by the WSM (World Stress Map) project represent the observed ori-entations of the maximum horizontal principal stress (sHmax) in a certain region. Assuming that the long-w...The 10 920 stress indicators collected so far by the WSM (World Stress Map) project represent the observed ori-entations of the maximum horizontal principal stress (sHmax) in a certain region. Assuming that the long-wave component of sHmax is expressed by the absolute direction of plate motions, we can get the relative orientation and the magnitude of the short-wave component resulted from the local tectonic process or other factors with vector analytical technique. The global surface was divided into basic element bins by 2.52.5 dimensions and the WSM indicators were statistically analyzed for each element by weight coefficient method in order to determine the mean orientation of the stress. We calculated the long-wave component of the global stress field using HS2-NUVEL1 model. The relative magnitude or the direction limitation of short-wave component, which reflect the local contribution to the observed stresses, was determined by the angle between the mean sHmax and the orien-tation of the long-wave component. The results of this paper show that the contribution of either the long-wave component or the short-wave component is approximately equal to most of the global plates on the basis of the mean effect of the observed stresses. For some of continental regions, the local active tectonics plays an important role in the observed stresses and controls the generation and occurrence of earthquakes.展开更多
An accurate numerical simulation for wave equations is essential for understanding of wave propagation in the earth's interior as well as full waveform inversion and reverse time migration. However, due to computa...An accurate numerical simulation for wave equations is essential for understanding of wave propagation in the earth's interior as well as full waveform inversion and reverse time migration. However, due to computational cost and hardware capability limitations, numerical simulations are often performed within a finite domain. Thus, an adequate absorbing boundary condition (ABC) is indispensable for obtaining accurate numerical simulation results. In this study, we develop a hybrid ABC based on a transmitting boundary, which is referred to as THABC, to eliminate artificial boundary reflections in 3D second-order fractional viscoacoustic numerical simulations. Furthermore, we propose an adaptive weighted coefficient to reconcile the transmitting and viscoacoustic wavefields in THABC. Through several numerical examples, we determine that the proposed THABC approach is characterized by the following benefits. First, with the same number of absorbing layers, THABC exhibits a better ability in eliminating boundary reflection than traditional ABC schemes. Second, THABC is more effective in computation, since it only requires the wavefields at the current and last time steps to solve the transmitting formula within the absorbing layers. Benefiting from a simple but effective combination between the transmitting equation and the second-order wave equation, our scheme performs well in the 3D fractional Laplacian viscoacoustic numerical simulation.展开更多
基金financially supported by National Key R&D Program of China(Grant No.2020YFA0711802)National Nature Science Foundation of China(Grant Nos.51439008 and 51779248).
文摘A typical blasting vibration wave is a composite wave,and its attenuation law is affected by the type of dominant wave component.The purpose of the present study is to establish an attenuation equation of the peak particle velocity(PPV),taking into account the attenuation characteristics of P-,S-and R-waves in the blasting vibration wave.Field blasting tests were carried out as a case to specifically apply the proposed equation.In view of the fact that the discrete properties of rock mass will inevitably cause the uncertainty of blasting vibration,we also carried out a probability analysis of PPV uncertainty,and introduced the concept of reliability to evaluate blasting vibration.The results showed that the established attenuation equation had a higher prediction accuracy,and can be considered as a promising equation implemented on more complex sites.The adopted uncertainty analysis method can comprehensively take account of the attenuation law of blasting vibration measured on site and discrete properties of rock masses.The obtained distribution of the PPV uncertainty factor can quantitatively evaluate the reliability of blasting vibration,which is a powerful and necessary supplement to the PPV attenuation equation.
基金supported by NSF grant EAR-063566(F.N.)National Natural Science Foundation of China grant 40774042(J.L.)
文摘The recently built China Digital Seismic Network consists of the China National Digital Seismic Network (CNDSN),31 regional seismic networks and several small aperture arrays with more than 1 000 stations including 850+ broadband stations.It forms a gigantic seismic array that provides an unprecedented opportunity to study the Earth's deep interior besides its routine task of seismic monitoring.Many modern seismic studies rely on rotation of vertical and horizontal components in order to separate different types of seismic waves.Knowledge of the orientations of the two horizontal components thus is important to perform a correction rotation.We analyzed particle motions of teleseismic P waves recorded by the network and used them to estimate the northcomponent azimuth of each station.An SNR-weighted-multi-event method was introduced to obtain component azimuths that best explain the P-wave particle motions of all the events recorded at a station.The method provides robust estimates including a measurement error calculated from background noise levels.We found that about one third of the stations have some sort of problems,including misorientation of the two horizontal components,mislabeling and polarity reversal in one or more components.These problems need to be taken into account for any rotation based seismic studies.
基金Project supported by the National Natural Science Foundation of China(Grant No.11374230)
文摘Ultrasonic Lamb waves are considered as a sensitive and effective tool for nondestructive testing and evaluation of plate-like or pipe-like structures. The nature of multimode and dispersion causes the wave packets to spread, and the modes overlap in both time and frequency domains as they propagate through the structures. By using a two-component laser interferometer technique, in combination with a priori knowledge of the dispersion characteristics and wave structure information of Lamb wave modes, a two-component signal processing technique is presented for implementing dispersion removal and mode separation simultaneously for two modes mixture signals of Lamb waves. The proposed algorithm is first processed and verified using synthetic Lamb wave signals. Then, the two-component displacements test experiment is conducted using different aluminum plate samples. Moreover, we confirm the effectiveness and robustness of this method.
基金supported by National Natural Science Foundation of China(Nos.11303099,41531071 and 41574158)the Youth Innovation Promotion Association CAS
文摘The dispersion relation and electromagnetic polarization of the plasma waves are comprehensively studied in cold electron,proton,and heavy charged particle plasmas.Three modes are classified as the fast,intermediate,and slow mode waves according to different phase velocities.When plasmas contain positively-charged particles,the fast and intermediate modes can interact at the small propagating angles,whereas the two modes are separate at the large propagating angles.The near-parallel intermediate and slow waves experience the linear polarization,circular polarization,and linear polarization again,with the increasing wave number.The wave number regime corresponding to the above circular polarization shrinks as the propagating angle increases.Moreover,the fast and intermediate modes cause the reverse change of the electromagnetic polarization at the special wave number.While the heavy particles carry the negative charges,the dispersion relations of the fast and intermediate modes are always separate,being independent of the propagating angles.Furthermore,this study gives new expressions of the three resonance frequencies corresponding to the highly-oblique propagation waves in the general three-component plasmas,and shows the dependence of the resonance frequencies on the propagating angle,the concentration of the heavy particle,and the mass ratio among different kinds of particles.
基金supported by the National Natural Science Foundation of China(Grant Nos.11371248,11431008,11271254,11428102,and 11671255)the Fund from the Ministry of Economy and Competitiveness of Spain(Grant Nos.MTM2012-37070 and MTM2016-80276-P(AEI/FEDER,EU))
文摘The nonlinear Schr?dinger(NLS) equation and Boussinesq equation are two very important integrable equations.They have widely physical applications. In this paper, we investigate a nonlinear system, which is the two-component NLS equation coupled to the Boussinesq equation. We obtain the bright–bright, bright–dark, and dark–dark soliton solutions to the nonlinear system. We discuss the collision between two solitons. We observe that the collision of bright–bright soliton is inelastic and two solitons oscillating periodically can happen in the two parallel-traveling bright–bright or bright–dark soliton solution. The general breather and rogue wave solutions are also given. Our results show again that there are more abundant dynamical properties for multi-component nonlinear systems.
基金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.
基金Project supported by the Global Change Research Program of China(Grant No.2015CB953904)the National Natural Science Foundation of China(Grant Nos.11275072 and 11435005)+2 种基金the Doctoral Program of Higher Education of China(Grant No.20120076110024)the Network Information Physics Calculation of Basic Research Innovation Research Group of China(Grant No.61321064)Shanghai Collaborative Innovation Center of Trustworthy Software for Internet of Things,China(Grant No.ZF1213)
文摘We study the generalized Darboux transformation to the three-component coupled nonlinear Schr ¨odinger equation.First-and second-order localized waves are obtained by this technique.In first-order localized wave,we get the interactional solutions between first-order rogue wave and one-dark,one-bright soliton respectively.Meanwhile,the interactional solutions between one-breather and first-order rogue wave are also given.In second-order localized wave,one-dark-one-bright soliton together with second-order rogue wave is presented in the first component,and two-bright soliton together with second-order rogue wave are gained respectively in the other two components.Besides,we observe second-order rogue wave together with one-breather in three components.Moreover,by increasing the absolute values of two free parameters,the nonlinear waves merge with each other distinctly.These results further reveal the interesting dynamic structures of localized waves in the three-component coupled system.
基金sponsored by National Natural Science Foundation of China(No.41672325,41602334)National Key Research and Development Program of China(No.2017YFC0601505).
文摘Many different effective reflection information are often contaminated by exterior and random noise which concealed in the seismic data.Traditional single or fixed transform is not suit for exploiting their complicated characteristics and attenuating the noise.Recent years,a novel method so-called morphological component analysis(MCA)is put forward to separate different geometrical components by amalgamating several irrelevance transforms.According to study the local singular and smooth linear components characteristics of seismic data,we propose a method of suppressing noise by integrating with the advantages of adaptive K-singular value decomposition(K-SVD)and wave atom dictionaries to depict the morphological features diversity of seismic signals.Numerical results indicate that our method can dramatically suppress the undesired noises,preserve the information of geologic body and geological structure and improve the signal-to-noise ratio of the data.We also demonstrate the superior performance of this approach by comparing with other novel dictionaries such as discrete cosine transform(DCT),undecimated discrete wavelet transform(UDWT),or curvelet transform,etc.This algorithm provides new ideas for data processing to advance quality and signal-to-noise ratio of seismic data.
基金MOST contract of 2001BA601B02 and State Natural Science Foundation of China (49804006).
文摘The 10 920 stress indicators collected so far by the WSM (World Stress Map) project represent the observed ori-entations of the maximum horizontal principal stress (sHmax) in a certain region. Assuming that the long-wave component of sHmax is expressed by the absolute direction of plate motions, we can get the relative orientation and the magnitude of the short-wave component resulted from the local tectonic process or other factors with vector analytical technique. The global surface was divided into basic element bins by 2.52.5 dimensions and the WSM indicators were statistically analyzed for each element by weight coefficient method in order to determine the mean orientation of the stress. We calculated the long-wave component of the global stress field using HS2-NUVEL1 model. The relative magnitude or the direction limitation of short-wave component, which reflect the local contribution to the observed stresses, was determined by the angle between the mean sHmax and the orien-tation of the long-wave component. The results of this paper show that the contribution of either the long-wave component or the short-wave component is approximately equal to most of the global plates on the basis of the mean effect of the observed stresses. For some of continental regions, the local active tectonics plays an important role in the observed stresses and controls the generation and occurrence of earthquakes.
基金National Natural Science Foundation of China under Grant Nos.41930431 and 41974116Natural Science Foundation of Heilongjiang Province No.YQ2021D008CNPC Innovation Found No.2021DQ02-0302 for supporting this work.
文摘An accurate numerical simulation for wave equations is essential for understanding of wave propagation in the earth's interior as well as full waveform inversion and reverse time migration. However, due to computational cost and hardware capability limitations, numerical simulations are often performed within a finite domain. Thus, an adequate absorbing boundary condition (ABC) is indispensable for obtaining accurate numerical simulation results. In this study, we develop a hybrid ABC based on a transmitting boundary, which is referred to as THABC, to eliminate artificial boundary reflections in 3D second-order fractional viscoacoustic numerical simulations. Furthermore, we propose an adaptive weighted coefficient to reconcile the transmitting and viscoacoustic wavefields in THABC. Through several numerical examples, we determine that the proposed THABC approach is characterized by the following benefits. First, with the same number of absorbing layers, THABC exhibits a better ability in eliminating boundary reflection than traditional ABC schemes. Second, THABC is more effective in computation, since it only requires the wavefields at the current and last time steps to solve the transmitting formula within the absorbing layers. Benefiting from a simple but effective combination between the transmitting equation and the second-order wave equation, our scheme performs well in the 3D fractional Laplacian viscoacoustic numerical simulation.