Based on the fifth-order Stokes regular wave theory, a simplified model for extreme-wave kinematics in deep sea was developed. In this model, from the wave records the average of two neighboring wave periods for the e...Based on the fifth-order Stokes regular wave theory, a simplified model for extreme-wave kinematics in deep sea was developed. In this model, from the wave records the average of two neighboring wave periods for the extreme crest or trough was defined as the period of the Stokes wave by the up and down zero-crossing methods. Then the input wave amplitude was deduced by substituting the wave period and extreme crest or trough into the expression for the fifth-order Stokes wave elevation. Thus the corresponding formula for the wave velocity can be used to describe kinematics beneath the extreme wave. By comparison with the published numerical models and experimental data, the proposed model is validated to be able to calculate the extreme wave velocity rather easily and accurately.展开更多
Engineering seismic exploration aims at shallow imaging which is confused by statics if the surface is uneven. Direct pre-stack depth migration （DPDM） is based on accurate elevations of sources and receivers, by whi...Engineering seismic exploration aims at shallow imaging which is confused by statics if the surface is uneven. Direct pre-stack depth migration （DPDM） is based on accurate elevations of sources and receivers, by which static correction is completely abandoned before migration and surely the imaging quality is remarkably improved. To obtain some artificial shot gathers, high-order staggered-grid finite-difference （FD） method is adapted to model acoustic wave propagation. Since the shot gathers are always disturbed by regular interferences, the statics still must be applied to supporting the interference elimination by apparent velocity filtering method. Then all the shot gathers should be removed back to their original positions by reverse statics. Finally, they are migrated by pre-stack reverse-time depth migration and imaged. The numerical experiments show that the DPDM can ideally avoid the mistakes caused by statics and increase imaging precision.展开更多
For microseisimic monitoring it is difficult to determine wave modes and their propagation velocity. In this paper, we propose a new method for automatically inverting in real time the source characteristics of micros...For microseisimic monitoring it is difficult to determine wave modes and their propagation velocity. In this paper, we propose a new method for automatically inverting in real time the source characteristics of microseismic events in mine engineering without wave mode identification and velocities. Based on the wave equation in a spherical coordinate system, we derive a tomographic imaging equation and formulate a scanning parameter selection criterion by which the microseisimic event maximum energy and corresponding parameters can be determined. By determining the maximum energy positions inside a given risk district, we can indentify microseismic events inside or outside the risk districts. The synthetic and field examples demonstrate that the proposed tomographic imaging method can automatically position microseismic events by only knowing the risk district dimensions and range of velocities without identifying the wavefield modes and accurate velocities. Therefore, the new method utilizes the full wavefields to automatically monitor microseismic events.展开更多
Deduced the propagation rule of longitudinal and transverse wave. On the basis of this, propagation rules in attenuated visco-elastic media and varied Lame coefficient were put forward as well. The subsequent numerica...Deduced the propagation rule of longitudinal and transverse wave. On the basis of this, propagation rules in attenuated visco-elastic media and varied Lame coefficient were put forward as well. The subsequent numerical analysis found that in a small scope longitudinal and transverse wave could be considered as homogeneously propagating when faultages and joints were not taken into account. The existence of lane hindered the wave＇s propagation, and it made the velocity gradient change in a local vicinity area. Therefore velocity varied in different direction.展开更多
An experimental study and theoretical analysis were carried out to explore the ground-borne vibration generated by elevated high-speed railway in rock strata.Taking a typical rail line constructed on rock area in Chin...An experimental study and theoretical analysis were carried out to explore the ground-borne vibration generated by elevated high-speed railway in rock strata.Taking a typical rail line constructed on rock area in China as the research object,a set of field tests was performed on Rizhao-Lankao High-Speed Railway,the bridge and ground vibrations were measured as trains passed at 330−340 km/h,then the transferring law and spatial distribution under individual frequencies were investigated.The experiment results indicate that the bridge frequency spectrum exhibited relatively high-frequency vibration peaks caused by short-wavelength irregularity;ground vibration farther than 30 m away can be amplified with a higher frequency and numerous components.Furthermore,the wave propagation equation of a stratified rock strata was established based on direct-stiffness method to explore the vibration attenuation rules via frequency-domain analysis.It is found that the rock area has a weaker correlation between vibration transmissibility and frequency,thicker and harder rock strata loss their vibration attenuation capacity.It can be concluded that the high-speed railways induced vibration on rock strata shows a wide frequency band and large amplitude,the design of reducing vibration aimed at specific frequency is important according to next more detailed numerical study.展开更多
The conventional fast converted-wave imaging method directly uses backward Pand converted S-wavefield to produce joint images. However, this image is accompanied by strong background noises, because the wavefi elds in...The conventional fast converted-wave imaging method directly uses backward Pand converted S-wavefield to produce joint images. However, this image is accompanied by strong background noises, because the wavefi elds in all propagation directions contribute to it. Given this issue, we improve the conventional imaging method in the two aspects. First, the amplitude-preserved P-and S-wavef ield are obtained by using an improved space-domain wavef ield separation scheme to decouple the original elastic wavef ield. Second, a convertedwave imaging condition is constructed based on the directional-wavefield separation and only the wavefields propagating in the same directions used for cross-correlation imaging, resulting in effectively eliminating the imaging artifacts of the wavefields with different directions;Complex-wavefi eld extrapolation is adopted to decompose the decoupled P-and S-wavefield into directional-wavefields during backward propagation, this improves the eff iciency of the directional-wavef ield separation. Experiments on synthetic data show that the improved method generates more accurate converted-wave images than the conventional one. Moreover, the improved method has application potential in micro-seismic and passive-source exploration due to its source-independent characteristic.展开更多
Based on the Tang-Othmer Ca^2＋ model, the drift behavior of intracellular Ca^2＋ spiral waves under the influence of weak electric field is investigated. Numerical results show that the dependence of drift velocity o...Based on the Tang-Othmer Ca^2＋ model, the drift behavior of intracellular Ca^2＋ spiral waves under the influence of weak electric field is investigated. Numerical results show that the dependence of drift velocity of the spiral tip on dc electric field is similar to experimental observations in BZ system. When an ac electric field is applied, interesting resonant-drift phenomenon is observed with ω=2ω0. All results can be explained analytically using a proximate method.展开更多
Theoretical analysis conducted of uniaxial compressive strength and tensile strength of artificially frozen soil and P waves and S waves and of the relationship between the two. Experiments are made on frozen sand and...Theoretical analysis conducted of uniaxial compressive strength and tensile strength of artificially frozen soil and P waves and S waves and of the relationship between the two. Experiments are made on frozen sand and frozen clay respectively at the temperature of - 7℃, - 12℃and - 17℃. Of the data obtained, regression analysis and gray-system correlation are conducted. As indicated by the results, the frozen soil tensile strength is closely correIated with the P wave velocity and the compressive with the S wave, hence the former is well described by the latter.展开更多
Recent studies using field case history data yielded new criteria for evaluating liquefaction potential in saturated granular deposits based on in situ, stress-corrected shear wave velocity. However, the conditions of...Recent studies using field case history data yielded new criteria for evaluating liquefaction potential in saturated granular deposits based on in situ, stress-corrected shear wave velocity. However, the conditions of relatively insufficient case histories and limited site conditions in this approach call for additional data to more reliably define liquefaction resistance as a function of shear wave velocity. In this study, a series of undrained cyclic triaxial tests were conducted on saturated sand with shear wave velocity Vs measured by bender element. By normalizing the data with respect to minimum void ratio, the test results, in-corporated with previously published laboratory data, statistically revealed good correlation of cyclic shear strength with small-strain shear modulus for sandy soils, which is almost irrespective of soil types and confining pressures. The consequently determined cyclic resistance ratio, CRR, was found to be approximately proportional to Vs4. Liquefaction resistance boundary curves were established by applying this relationship and compared to liquefaction criteria derived from seismic field measure-ments. Although in the range of Vs1>200 m/s the presented curves are moderately conservative, they are remarkably consistent with the published field performance criteria on the whole.展开更多
Under the excitation of elastic waves,local fluid flow in a complex porous medium is a major cause for wave dispersion and attenuation.When the local fluid flow process is simulated with wave propagation equations in ...Under the excitation of elastic waves,local fluid flow in a complex porous medium is a major cause for wave dispersion and attenuation.When the local fluid flow process is simulated with wave propagation equations in the double-porosity medium,two porous skeletons are usually assumed,namely,host and inclusions.Of them,the volume ratio of inclusion skeletons is low.All previous studies have ignored the consideration of local fluid flow velocity field in inclusions,and therefore they can not completely describe the physical process of local flow oscillation and should not be applied to the situation where the fluid kinetic energy in inclusions cannot be neglected.In this paper,we analyze the local fluid flow velocity fields inside and outside the inclusion,rewrite the kinetic energy function and dissipation function based on the double-porosity medium model containing spherical inclusions,and derive the reformulated Biot-Rayleigh(BR)equations of elastic wave propagation based on Hamilton’s principle.We present simulation examples with different rock and fluid types.Comparisons between BR equations and reformulated BR equations show that there are significant differences in wave response characteristics.Finally,we compare the reformulated BR equations with the previous theories and experimental data,and the results show that the theoretical results of this paper are correct and effective.展开更多
We present a scheme to simulate SH-wave propagation in a whole-Earth model with arbitrary lateral heterogeneities employing the Fourier pseudospectral method. Wave equations are defined in two-dimensional cylindrical ...We present a scheme to simulate SH-wave propagation in a whole-Earth model with arbitrary lateral heterogeneities employing the Fourier pseudospectral method. Wave equations are defined in two-dimensional cylindrical coordinates and the model is taken through a great circle of the Earth. Spatial derivatives in the wave equations are calculated in the wavenumber domain by multiplication, and the transformation between spatial and wavenumber domains is performed via fast Fourier transformation. Because of the high accuracy and high speed of the Fourier pseudospectral method, the scheme enables us to calculate a short-wavelength global SH-wavefield with accurate waveforms and arrival times for models with heterogeneities that can be approximated as azimuthally symmetric. Comparing with two-dimensional simulation methods based on an axisymmetric model, implementing the seismic source in the present scheme is more convenient. We calculated the global SH-wavefield for the preliminary reference Earth model to identify the generation, reflection and refraction of various seismic phases propagating in the Earth. Applications to a heterogeneous global model with low-velocity perturbation above the core-mantle boundary were conducted to analyze the effect of lateral heterogeneity on global SH-wave propagation.展开更多
By making use of the weak gravitational field approximation, we obtain a linearized solution of the gravitational vacuum field equation in an anisotropic spacetime. The plane-wave solution and dispersion relation of g...By making use of the weak gravitational field approximation, we obtain a linearized solution of the gravitational vacuum field equation in an anisotropic spacetime. The plane-wave solution and dispersion relation of gravitationaJ wave is presented explicitly. There is possibility that the speed of gravitational wave is larger than the speed of light and the easuality still holds. We show that the energy-momentum of gravitational wave in the ansiotropic spacetime is still well defined and conserved.展开更多
基金Supported by the NSFC (under Grant Nos.5070900 and 10772040)the National High Tech Research and Development Program of China (2006AA09A109-3)
文摘Based on the fifth-order Stokes regular wave theory, a simplified model for extreme-wave kinematics in deep sea was developed. In this model, from the wave records the average of two neighboring wave periods for the extreme crest or trough was defined as the period of the Stokes wave by the up and down zero-crossing methods. Then the input wave amplitude was deduced by substituting the wave period and extreme crest or trough into the expression for the fifth-order Stokes wave elevation. Thus the corresponding formula for the wave velocity can be used to describe kinematics beneath the extreme wave. By comparison with the published numerical models and experimental data, the proposed model is validated to be able to calculate the extreme wave velocity rather easily and accurately.
文摘Engineering seismic exploration aims at shallow imaging which is confused by statics if the surface is uneven. Direct pre-stack depth migration （DPDM） is based on accurate elevations of sources and receivers, by which static correction is completely abandoned before migration and surely the imaging quality is remarkably improved. To obtain some artificial shot gathers, high-order staggered-grid finite-difference （FD） method is adapted to model acoustic wave propagation. Since the shot gathers are always disturbed by regular interferences, the statics still must be applied to supporting the interference elimination by apparent velocity filtering method. Then all the shot gathers should be removed back to their original positions by reverse statics. Finally, they are migrated by pre-stack reverse-time depth migration and imaged. The numerical experiments show that the DPDM can ideally avoid the mistakes caused by statics and increase imaging precision.
基金support jointly by projects of the National Natural Science Fund Project (40674017 and 50774012)the National Key Basic Research and Development Plan 973 (2010CB226803)
文摘For microseisimic monitoring it is difficult to determine wave modes and their propagation velocity. In this paper, we propose a new method for automatically inverting in real time the source characteristics of microseismic events in mine engineering without wave mode identification and velocities. Based on the wave equation in a spherical coordinate system, we derive a tomographic imaging equation and formulate a scanning parameter selection criterion by which the microseisimic event maximum energy and corresponding parameters can be determined. By determining the maximum energy positions inside a given risk district, we can indentify microseismic events inside or outside the risk districts. The synthetic and field examples demonstrate that the proposed tomographic imaging method can automatically position microseismic events by only knowing the risk district dimensions and range of velocities without identifying the wavefield modes and accurate velocities. Therefore, the new method utilizes the full wavefields to automatically monitor microseismic events.
基金Supported by International Important Cooperation Project of National Natural Science Foundation of China(50320120001).
文摘Deduced the propagation rule of longitudinal and transverse wave. On the basis of this, propagation rules in attenuated visco-elastic media and varied Lame coefficient were put forward as well. The subsequent numerical analysis found that in a small scope longitudinal and transverse wave could be considered as homogeneously propagating when faultages and joints were not taken into account. The existence of lane hindered the wave＇s propagation, and it made the velocity gradient change in a local vicinity area. Therefore velocity varied in different direction.
基金Project(2016YFE0205200)supported by the National Key Research and Development Program of ChinaProjects(U1734207,51978585)supported by the National Natural Science Foundation of China。
文摘An experimental study and theoretical analysis were carried out to explore the ground-borne vibration generated by elevated high-speed railway in rock strata.Taking a typical rail line constructed on rock area in China as the research object,a set of field tests was performed on Rizhao-Lankao High-Speed Railway,the bridge and ground vibrations were measured as trains passed at 330−340 km/h,then the transferring law and spatial distribution under individual frequencies were investigated.The experiment results indicate that the bridge frequency spectrum exhibited relatively high-frequency vibration peaks caused by short-wavelength irregularity;ground vibration farther than 30 m away can be amplified with a higher frequency and numerous components.Furthermore,the wave propagation equation of a stratified rock strata was established based on direct-stiffness method to explore the vibration attenuation rules via frequency-domain analysis.It is found that the rock area has a weaker correlation between vibration transmissibility and frequency,thicker and harder rock strata loss their vibration attenuation capacity.It can be concluded that the high-speed railways induced vibration on rock strata shows a wide frequency band and large amplitude,the design of reducing vibration aimed at specific frequency is important according to next more detailed numerical study.
基金supported by the National Science and Technology Major Project of China(No.2017ZX05018-005)National Natural Science Foundation of China(No.41474110)
文摘The conventional fast converted-wave imaging method directly uses backward Pand converted S-wavefield to produce joint images. However, this image is accompanied by strong background noises, because the wavefi elds in all propagation directions contribute to it. Given this issue, we improve the conventional imaging method in the two aspects. First, the amplitude-preserved P-and S-wavef ield are obtained by using an improved space-domain wavef ield separation scheme to decouple the original elastic wavef ield. Second, a convertedwave imaging condition is constructed based on the directional-wavefield separation and only the wavefields propagating in the same directions used for cross-correlation imaging, resulting in effectively eliminating the imaging artifacts of the wavefields with different directions;Complex-wavefi eld extrapolation is adopted to decompose the decoupled P-and S-wavefield into directional-wavefields during backward propagation, this improves the eff iciency of the directional-wavef ield separation. Experiments on synthetic data show that the improved method generates more accurate converted-wave images than the conventional one. Moreover, the improved method has application potential in micro-seismic and passive-source exploration due to its source-independent characteristic.
基金Supported by National Natural Science Foundation of China under Grant Nos.10575041 and 10747005
文摘Based on the Tang-Othmer Ca^2＋ model, the drift behavior of intracellular Ca^2＋ spiral waves under the influence of weak electric field is investigated. Numerical results show that the dependence of drift velocity of the spiral tip on dc electric field is similar to experimental observations in BZ system. When an ac electric field is applied, interesting resonant-drift phenomenon is observed with ω=2ω0. All results can be explained analytically using a proximate method.
文摘Theoretical analysis conducted of uniaxial compressive strength and tensile strength of artificially frozen soil and P waves and S waves and of the relationship between the two. Experiments are made on frozen sand and frozen clay respectively at the temperature of - 7℃, - 12℃and - 17℃. Of the data obtained, regression analysis and gray-system correlation are conducted. As indicated by the results, the frozen soil tensile strength is closely correIated with the P wave velocity and the compressive with the S wave, hence the former is well described by the latter.
基金Project supported by the National Natural Science Foundation of China (No. 10372089), and Department of Education of Zhejiang Province (No. 20010572), China
文摘Recent studies using field case history data yielded new criteria for evaluating liquefaction potential in saturated granular deposits based on in situ, stress-corrected shear wave velocity. However, the conditions of relatively insufficient case histories and limited site conditions in this approach call for additional data to more reliably define liquefaction resistance as a function of shear wave velocity. In this study, a series of undrained cyclic triaxial tests were conducted on saturated sand with shear wave velocity Vs measured by bender element. By normalizing the data with respect to minimum void ratio, the test results, in-corporated with previously published laboratory data, statistically revealed good correlation of cyclic shear strength with small-strain shear modulus for sandy soils, which is almost irrespective of soil types and confining pressures. The consequently determined cyclic resistance ratio, CRR, was found to be approximately proportional to Vs4. Liquefaction resistance boundary curves were established by applying this relationship and compared to liquefaction criteria derived from seismic field measure-ments. Although in the range of Vs1>200 m/s the presented curves are moderately conservative, they are remarkably consistent with the published field performance criteria on the whole.
基金supported by the National Natural Science Foundation of China(Grant No.41104066)RIPED Youth Innovation Foundation(Grant No.2010-A-26-01)+1 种基金the National Basic Research Program of China(Grant No.2014CB239006)the Open fund of SINOPEC Key Laboratory of Geophysics(Grant No.WTYJY-WX2013-04-18)
文摘Under the excitation of elastic waves,local fluid flow in a complex porous medium is a major cause for wave dispersion and attenuation.When the local fluid flow process is simulated with wave propagation equations in the double-porosity medium,two porous skeletons are usually assumed,namely,host and inclusions.Of them,the volume ratio of inclusion skeletons is low.All previous studies have ignored the consideration of local fluid flow velocity field in inclusions,and therefore they can not completely describe the physical process of local flow oscillation and should not be applied to the situation where the fluid kinetic energy in inclusions cannot be neglected.In this paper,we analyze the local fluid flow velocity fields inside and outside the inclusion,rewrite the kinetic energy function and dissipation function based on the double-porosity medium model containing spherical inclusions,and derive the reformulated Biot-Rayleigh(BR)equations of elastic wave propagation based on Hamilton’s principle.We present simulation examples with different rock and fluid types.Comparisons between BR equations and reformulated BR equations show that there are significant differences in wave response characteristics.Finally,we compare the reformulated BR equations with the previous theories and experimental data,and the results show that the theoretical results of this paper are correct and effective.
基金supported by National Natural Science Foundation of China (Grant Nos. 40874020, 40474012 and 40821062)National R&D Special Fund for Public Welfare Industry (Grant No. 20070804)
文摘We present a scheme to simulate SH-wave propagation in a whole-Earth model with arbitrary lateral heterogeneities employing the Fourier pseudospectral method. Wave equations are defined in two-dimensional cylindrical coordinates and the model is taken through a great circle of the Earth. Spatial derivatives in the wave equations are calculated in the wavenumber domain by multiplication, and the transformation between spatial and wavenumber domains is performed via fast Fourier transformation. Because of the high accuracy and high speed of the Fourier pseudospectral method, the scheme enables us to calculate a short-wavelength global SH-wavefield with accurate waveforms and arrival times for models with heterogeneities that can be approximated as azimuthally symmetric. Comparing with two-dimensional simulation methods based on an axisymmetric model, implementing the seismic source in the present scheme is more convenient. We calculated the global SH-wavefield for the preliminary reference Earth model to identify the generation, reflection and refraction of various seismic phases propagating in the Earth. Applications to a heterogeneous global model with low-velocity perturbation above the core-mantle boundary were conducted to analyze the effect of lateral heterogeneity on global SH-wave propagation.
基金Supported by the National Natural Science Foundation of China under Grant Nos.1087512911075166and 11147176
文摘By making use of the weak gravitational field approximation, we obtain a linearized solution of the gravitational vacuum field equation in an anisotropic spacetime. The plane-wave solution and dispersion relation of gravitationaJ wave is presented explicitly. There is possibility that the speed of gravitational wave is larger than the speed of light and the easuality still holds. We show that the energy-momentum of gravitational wave in the ansiotropic spacetime is still well defined and conserved.
