Elastic wave on seafloor caused by low frequency noise radiated from ship is called ship seismic wave which can be used to identify ship target. In order to analyze the wave components and the propagating properties o...Elastic wave on seafloor caused by low frequency noise radiated from ship is called ship seismic wave which can be used to identify ship target. In order to analyze the wave components and the propagating properties of ship seismic wave, the numerical calculation of synthetic seismograms on seafloor aroused by a low frequency point sound source is carried out using a wave number integration technique combined with inverse Fourier transform. According to the numerical example of hard seafloor, the time series of seismic wave on seafloor are mostly composed of interface waves and normal mode waves. Each normal mode wave has a well defined low cut-off frequency, while the interface wave doesn't have. The frequency dispersion of normal mode wave is obvious when frequency is lower than 100Hz, while the interface wave is dispersive only in the infra-sound frequency range. The time series of seismic wave is dominated by the interface wave when the source frequency is less than the minimal cut-off frequency of normal mode wave.展开更多
In the past two decades numerous studies were made to develop and improve the theory and practical computation techniques of synthesizing theoretical seismograms for the model of multi-layered half-space. Today, synth...In the past two decades numerous studies were made to develop and improve the theory and practical computation techniques of synthesizing theoretical seismograms for the model of multi-layered half-space. Today, synthesizing theoretical seismograms in multi-layered half-space is an important tool for understanding the structure of the Earth’s interior as well as the seismic source process from well-recorded seismograms data. As part of a review of the state-of-the-art, in this article I shall present a systematic and self-contained theory of elastic waves in multi-layered half-space media based on the developments in the past two decades.展开更多
A novel technique for automatic seismic data processing using both integral and local feature of seismograms was presented in this paper. Here, the term integral feature of seismograms refers to feature which may depi...A novel technique for automatic seismic data processing using both integral and local feature of seismograms was presented in this paper. Here, the term integral feature of seismograms refers to feature which may depict the shape of the whole seismograms. However, unlike some previous efforts which completely abandon the DIAL approach, i.e., signal detection, phase identifi- cation, association, and event localization, and seek to use envelope cross-correlation to detect seismic events directly, our technique keeps following the DIAL approach, but in addition to detect signals corresponding to individual seismic phases, it also detects continuous wave-trains and explores their feature for phase-type identification and signal association. More concrete ideas about how to define wave-trains and combine them with various detections, as well as how to measure and utilize their feature in the seismic data processing were expatiated in the paper. This approach has been applied to the routine data processing by us for years, and test results for a 16 days＇ period using data from the Xinjiang seismic station network were presented. The automatic processing results have fairly low false and missed event rate simultaneously, showing that the new technique has good application prospects for improvement of the automatic seismic data processing.展开更多
Synthetic seismograms generated by solving the seismic wave equation using numerical methods are being widely used in seismology. For fully three-dimensional seismic structure models, the generation of these synthetic...Synthetic seismograms generated by solving the seismic wave equation using numerical methods are being widely used in seismology. For fully three-dimensional seismic structure models, the generation of these synthetic seismograms may require large amount of computing resources. Conventional high-performance computer clusters may not provide a cost-effective solution to this type of applications. The newly emerging cloud-computing platform provides an alternative solution. In this paper, we describe our implementation of a synthetic seismogram generator based on the reciprocity principle using the Windows Azure cloud application framework. Our preliminary experiment shows that our cloud-based synthetic seismogram generator provides a costeffective and numerically efficient approach for computing synthetic seismograms based on the reciprocity principle.展开更多
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.展开更多
Based on the extensive-dilatancy anisotropy theory, the method of synthetic seismogram is used to estimate the anisotropic parameters. The advantages of the method lie in that it avoids the singularity resolution and ...Based on the extensive-dilatancy anisotropy theory, the method of synthetic seismogram is used to estimate the anisotropic parameters. The advantages of the method lie in that it avoids the singularity resolution and saves calculation time of computer by using the eigenvalue and eigenvector analytical expressions of Christoffel equation, at the same time, the result is tested by coherence function. The test result reveals that there exists a fine linear relation between original records and synthetic records, indicating the anisotropic parameters estimated by synthetic seismogram can reflect and describe the anisotropic characteristics of the given region medium.展开更多
A source-time function plays an important role in generating finite-difference synthetic seismograms.In this regard,we consider the amplitude spectra of five different time functions.This investigation shows that when...A source-time function plays an important role in generating finite-difference synthetic seismograms.In this regard,we consider the amplitude spectra of five different time functions.This investigation shows that when the high-frequency signals of the time function carry more power than the main signal,the tailing phenomena would be associated with the main signal which is not desirable in producing finite-difference synthetic seismograms.It is illustrated that the time function denoted by is the optimal source-time function among those examined in this paper for generating finite-difference synthetic seismograms.This matter is also verified by making finite-difference synthetic seismograms over an elastic earth model consisting of a layer over a half space.展开更多
In order to obtain a finite-difference synthetic seismogram, the finite earth model must be subdivided into cells.If the cell size becomes too large in comparison with wavelength of the source signal, waves disperse w...In order to obtain a finite-difference synthetic seismogram, the finite earth model must be subdivided into cells.If the cell size becomes too large in comparison with wavelength of the source signal, waves disperse with increasing traveltime. This phenomenon is known as grid dispersion. The grid dispersion produces a variation of velocity with different frequencies; that is the higher signal frequencies travel more slowly than the lower signal frequencies. Consequently, substantial tailing of the signal arises with increasing traveltime. This may happpen if: (1) the grid spacing is large; (2) the sampling rate is large; or (3) the source wavelength is too short compared with grid size. In other words, an important parameter in generating finite-difference synthetic seismograms is the number of grid points per wavelength of the source signal. In this paper, it is shown that the frequency of a source function has great effects on grid dispersion when P-Sv coupled waves propagate through elastic models by finite-differences. The two-dimensional elastic models considered in this paper consist of : (1) a normal fault, and (2) a layer over a half-space. This study verifies that, when generating finite-difference synthetic seismograms, the grid dispersion will be reduced to a satisfactory level if the grid points per wavelength at the half-power frequency of the source signal for the lowest velocity of the medium exceeds ten.展开更多
The extensive use of depth-imaged seismic data produced by pre-stack depth migration(PSDM)leads to the necessity to synthesize seismogram directly in depth domain.However,since seismic wavelet in depth domain is depen...The extensive use of depth-imaged seismic data produced by pre-stack depth migration(PSDM)leads to the necessity to synthesize seismogram directly in depth domain.However,since seismic wavelet in depth domain is dependent on media velocities.The time-domain convolution operation directly used in depth domain does not meet the linear time-invariant condition.In this paper,we present a new method for genuine depth-domain seismic synthesis.This method constructs the velocity-dependent seismic wavelets varying adaptively with the corresponding interval velocities in the depth direction and weights them by the reflectivities,then the synthetic seismic record is obtained by the superposition of these weighted depth-varying wavelets.We applied this method to synthesize the seismic record of both a multi-layered geological model and the field data.The results show that the method can accommodate the intrinsic velocity-dependent variation characteristics of seismic events in depth domain and avoid the redundant depth-to-time and time-to-depth transformations.展开更多
The S-wave velocity across the earth structure under Indonesia for Indonesia earthquakes has been investigated through seismogram analysis, simultaneously in the time domain and three Cartesian components. The data we...The S-wave velocity across the earth structure under Indonesia for Indonesia earthquakes has been investigated through seismogram analysis, simultaneously in the time domain and three Cartesian components. The data were recorded at DAV observational station, the Philippines. The main data set is the seismogram comparison between the measured and synthetic seismogram, instead of travel time data, as commonly used in other seismological research. The synthetic seismogram is calculated using the GEMINI method, which is equivalent to Mode Summation. The above seismogram comparison shows that the global earth mantle of PREMAN gives a deviating synthetic seismogram and has earlier arrival times than those of the measurement. The gradient of βh in the upper mantle layers is altered into a positive, rather than negative slope as stated in the PREMAN model, and negative corrections are imposed to the zero order of the polynomiars coefficients in all earth mantle layers. The excellent fitting, as well as travel time or waveform, is obtained from the surface waves of Love and Rayleigh, surface wave to the S and SS mantle waves as well as the core reflected waves. This result expresses that part of the earth mantle, due to a collision between India and Asia tectonic released zones, has a negative anomaly in S-wave velocity and vertical anisotropy in all of the earth mantle layers.展开更多
The S wave velocity structure between the hypocenter of C060394F earthquake,South Java and a series of observatory stations located in Australia and South-East Asia have been investigated through seismogram analysis i...The S wave velocity structure between the hypocenter of C060394F earthquake,South Java and a series of observatory stations located in Australia and South-East Asia have been investigated through seismogram analysis in the time domain and the three Cartesian components.The synthetic seismogram is constructed from the PREMAN global earth model.Seismogram comparison between the measured and synthetic seismograms shows large discrepancies.A correction to the S wave velocity structure is needed to solve these discrepancies.Seismogram analysis and fitting at the observatory stations located in Australia show a positive anomaly,whereas the stations in South-East Asia show a strong negative anomaly.Vertical anisotropy occurred in the upper mantle layer and earth mantle layers beneath it.The seismogram analysis in the time domain and the three components simultaneously provides a better picture of the earth model,compared to the analysis of the seismogram with arrival time,travel time difference and dispersion analysis.展开更多
The S wave velocity structure of the earth below Eastern Southeast Asia has been investigated by analyzing the seismogram from surface wave to multiple depth waves in the time domain and three Cartesian components sim...The S wave velocity structure of the earth below Eastern Southeast Asia has been investigated by analyzing the seismogram from surface wave to multiple depth waves in the time domain and three Cartesian components simultaneously.The wave passes across the front area of subduction zone between the Philippine plate and the Asian plate.The main data are waveform comparisons,instead of the arrival times.The synthetic seismogram is calculated using the GEMINI method.The synthetic seismogram constructed by PREMAN global earth model deviates greatly from the measured one.To solve this problem,corrections are needed for the β speed structure.Corrections cover the gradient change of β_h,which turns from negative to positive in upper mantle layers as in the PREMAN,change of earth crust depth and change of zero order coefficients of β velocity function in all earth mantle layers.So,the fitting is obtained,as well as the arrival time or the waveform of Love and Rayleigh surface waves,the S wave and the repetitive depth waves ScS_2 and ScS_3.This result reveals that the Southeast Asia,being stretched due to tectonic release,has a mantle in some parts with negative anomaly of S wave velocity and vertical anisotropy in all earth mantle layers.展开更多
Seismograms near source fault were synthesized using the hybrid empirical Green’s func-tion method where the discretely simulated seismic waveforms are used for Green’s functions instead of the observed waveforms of...Seismograms near source fault were synthesized using the hybrid empirical Green’s func-tion method where the discretely simulated seismic waveforms are used for Green’s functions instead of the observed waveforms of small earthquakes. The Green’s function seismic waveforms for small earth-quake were calculated by solving wave equation us-ing the pseudo-spectral method with the staggered grid real FFT strategy under a detailed 2-D velocity structure in Kobe region. Magnitude and seismic moment of simulated Green’s function waveforms were firstly determined by using the relationship be-tween fault length and corner frequency of source spectrum. The simulated Green’s function waveforms were employed to synthesize seismograms of strong ground motion near the earthquake fault. The syn-thetic seismograms of the target earthquake were performed based on the model with multiple source rupture processes. The results suggest that synthe-sized seismograms coincide well with observed seismic waveforms of the 1995 Hyogo-ken Nanbu earthquake. The simulated Green’s function method is very useful for prediction of the strong ground mo-tion in region without observed seismic waveforms. The present technique spreads application field of the empirical Green’s function method.展开更多
基金Sponsored by National Nature Science Foundation of China ( 51179195)National Defense Foundation of China ( 513030203-02)
文摘Elastic wave on seafloor caused by low frequency noise radiated from ship is called ship seismic wave which can be used to identify ship target. In order to analyze the wave components and the propagating properties of ship seismic wave, the numerical calculation of synthetic seismograms on seafloor aroused by a low frequency point sound source is carried out using a wave number integration technique combined with inverse Fourier transform. According to the numerical example of hard seafloor, the time series of seismic wave on seafloor are mostly composed of interface waves and normal mode waves. Each normal mode wave has a well defined low cut-off frequency, while the interface wave doesn't have. The frequency dispersion of normal mode wave is obvious when frequency is lower than 100Hz, while the interface wave is dispersive only in the infra-sound frequency range. The time series of seismic wave is dominated by the interface wave when the source frequency is less than the minimal cut-off frequency of normal mode wave.
文摘In the past two decades numerous studies were made to develop and improve the theory and practical computation techniques of synthesizing theoretical seismograms for the model of multi-layered half-space. Today, synthesizing theoretical seismograms in multi-layered half-space is an important tool for understanding the structure of the Earth’s interior as well as the seismic source process from well-recorded seismograms data. As part of a review of the state-of-the-art, in this article I shall present a systematic and self-contained theory of elastic waves in multi-layered half-space media based on the developments in the past two decades.
文摘A novel technique for automatic seismic data processing using both integral and local feature of seismograms was presented in this paper. Here, the term integral feature of seismograms refers to feature which may depict the shape of the whole seismograms. However, unlike some previous efforts which completely abandon the DIAL approach, i.e., signal detection, phase identifi- cation, association, and event localization, and seek to use envelope cross-correlation to detect seismic events directly, our technique keeps following the DIAL approach, but in addition to detect signals corresponding to individual seismic phases, it also detects continuous wave-trains and explores their feature for phase-type identification and signal association. More concrete ideas about how to define wave-trains and combine them with various detections, as well as how to measure and utilize their feature in the seismic data processing were expatiated in the paper. This approach has been applied to the routine data processing by us for years, and test results for a 16 days＇ period using data from the Xinjiang seismic station network were presented. The automatic processing results have fairly low false and missed event rate simultaneously, showing that the new technique has good application prospects for improvement of the automatic seismic data processing.
文摘Synthetic seismograms generated by solving the seismic wave equation using numerical methods are being widely used in seismology. For fully three-dimensional seismic structure models, the generation of these synthetic seismograms may require large amount of computing resources. Conventional high-performance computer clusters may not provide a cost-effective solution to this type of applications. The newly emerging cloud-computing platform provides an alternative solution. In this paper, we describe our implementation of a synthetic seismogram generator based on the reciprocity principle using the Windows Azure cloud application framework. Our preliminary experiment shows that our cloud-based synthetic seismogram generator provides a costeffective and numerically efficient approach for computing synthetic seismograms based on the reciprocity principle.
基金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.
基金National Natural Science Foundation of China (D4007410) and Joint Seismological Science Foundation of China (603006).
文摘Based on the extensive-dilatancy anisotropy theory, the method of synthetic seismogram is used to estimate the anisotropic parameters. The advantages of the method lie in that it avoids the singularity resolution and saves calculation time of computer by using the eigenvalue and eigenvector analytical expressions of Christoffel equation, at the same time, the result is tested by coherence function. The test result reveals that there exists a fine linear relation between original records and synthetic records, indicating the anisotropic parameters estimated by synthetic seismogram can reflect and describe the anisotropic characteristics of the given region medium.
文摘A source-time function plays an important role in generating finite-difference synthetic seismograms.In this regard,we consider the amplitude spectra of five different time functions.This investigation shows that when the high-frequency signals of the time function carry more power than the main signal,the tailing phenomena would be associated with the main signal which is not desirable in producing finite-difference synthetic seismograms.It is illustrated that the time function denoted by is the optimal source-time function among those examined in this paper for generating finite-difference synthetic seismograms.This matter is also verified by making finite-difference synthetic seismograms over an elastic earth model consisting of a layer over a half space.
文摘In order to obtain a finite-difference synthetic seismogram, the finite earth model must be subdivided into cells.If the cell size becomes too large in comparison with wavelength of the source signal, waves disperse with increasing traveltime. This phenomenon is known as grid dispersion. The grid dispersion produces a variation of velocity with different frequencies; that is the higher signal frequencies travel more slowly than the lower signal frequencies. Consequently, substantial tailing of the signal arises with increasing traveltime. This may happpen if: (1) the grid spacing is large; (2) the sampling rate is large; or (3) the source wavelength is too short compared with grid size. In other words, an important parameter in generating finite-difference synthetic seismograms is the number of grid points per wavelength of the source signal. In this paper, it is shown that the frequency of a source function has great effects on grid dispersion when P-Sv coupled waves propagate through elastic models by finite-differences. The two-dimensional elastic models considered in this paper consist of : (1) a normal fault, and (2) a layer over a half-space. This study verifies that, when generating finite-difference synthetic seismograms, the grid dispersion will be reduced to a satisfactory level if the grid points per wavelength at the half-power frequency of the source signal for the lowest velocity of the medium exceeds ten.
文摘The extensive use of depth-imaged seismic data produced by pre-stack depth migration(PSDM)leads to the necessity to synthesize seismogram directly in depth domain.However,since seismic wavelet in depth domain is dependent on media velocities.The time-domain convolution operation directly used in depth domain does not meet the linear time-invariant condition.In this paper,we present a new method for genuine depth-domain seismic synthesis.This method constructs the velocity-dependent seismic wavelets varying adaptively with the corresponding interval velocities in the depth direction and weights them by the reflectivities,then the synthetic seismic record is obtained by the superposition of these weighted depth-varying wavelets.We applied this method to synthesize the seismic record of both a multi-layered geological model and the field data.The results show that the method can accommodate the intrinsic velocity-dependent variation characteristics of seismic events in depth domain and avoid the redundant depth-to-time and time-to-depth transformations.
文摘The S-wave velocity across the earth structure under Indonesia for Indonesia earthquakes has been investigated through seismogram analysis, simultaneously in the time domain and three Cartesian components. The data were recorded at DAV observational station, the Philippines. The main data set is the seismogram comparison between the measured and synthetic seismogram, instead of travel time data, as commonly used in other seismological research. The synthetic seismogram is calculated using the GEMINI method, which is equivalent to Mode Summation. The above seismogram comparison shows that the global earth mantle of PREMAN gives a deviating synthetic seismogram and has earlier arrival times than those of the measurement. The gradient of βh in the upper mantle layers is altered into a positive, rather than negative slope as stated in the PREMAN model, and negative corrections are imposed to the zero order of the polynomiars coefficients in all earth mantle layers. The excellent fitting, as well as travel time or waveform, is obtained from the surface waves of Love and Rayleigh, surface wave to the S and SS mantle waves as well as the core reflected waves. This result expresses that part of the earth mantle, due to a collision between India and Asia tectonic released zones, has a negative anomaly in S-wave velocity and vertical anisotropy in all of the earth mantle layers.
文摘The S wave velocity structure between the hypocenter of C060394F earthquake,South Java and a series of observatory stations located in Australia and South-East Asia have been investigated through seismogram analysis in the time domain and the three Cartesian components.The synthetic seismogram is constructed from the PREMAN global earth model.Seismogram comparison between the measured and synthetic seismograms shows large discrepancies.A correction to the S wave velocity structure is needed to solve these discrepancies.Seismogram analysis and fitting at the observatory stations located in Australia show a positive anomaly,whereas the stations in South-East Asia show a strong negative anomaly.Vertical anisotropy occurred in the upper mantle layer and earth mantle layers beneath it.The seismogram analysis in the time domain and the three components simultaneously provides a better picture of the earth model,compared to the analysis of the seismogram with arrival time,travel time difference and dispersion analysis.
基金This research is funded by the DanaHibah Penelitian Dasar of DIKTI (Grant Fund of the Basic Research from Higher Education Agency) with grant number 019/SP3/PP/DP2M/Ⅱ/2006
文摘The S wave velocity structure of the earth below Eastern Southeast Asia has been investigated by analyzing the seismogram from surface wave to multiple depth waves in the time domain and three Cartesian components simultaneously.The wave passes across the front area of subduction zone between the Philippine plate and the Asian plate.The main data are waveform comparisons,instead of the arrival times.The synthetic seismogram is calculated using the GEMINI method.The synthetic seismogram constructed by PREMAN global earth model deviates greatly from the measured one.To solve this problem,corrections are needed for the β speed structure.Corrections cover the gradient change of β_h,which turns from negative to positive in upper mantle layers as in the PREMAN,change of earth crust depth and change of zero order coefficients of β velocity function in all earth mantle layers.So,the fitting is obtained,as well as the arrival time or the waveform of Love and Rayleigh surface waves,the S wave and the repetitive depth waves ScS_2 and ScS_3.This result reveals that the Southeast Asia,being stretched due to tectonic release,has a mantle in some parts with negative anomaly of S wave velocity and vertical anisotropy in all earth mantle layers.
基金Acknowledgements This work was supported by National Basic Research Program of China （Grant No. 2003CB716505） and the National Natural Science Foundation of China （Grant No. 40399141）.
文摘Seismograms near source fault were synthesized using the hybrid empirical Green’s func-tion method where the discretely simulated seismic waveforms are used for Green’s functions instead of the observed waveforms of small earthquakes. The Green’s function seismic waveforms for small earth-quake were calculated by solving wave equation us-ing the pseudo-spectral method with the staggered grid real FFT strategy under a detailed 2-D velocity structure in Kobe region. Magnitude and seismic moment of simulated Green’s function waveforms were firstly determined by using the relationship be-tween fault length and corner frequency of source spectrum. The simulated Green’s function waveforms were employed to synthesize seismograms of strong ground motion near the earthquake fault. The syn-thetic seismograms of the target earthquake were performed based on the model with multiple source rupture processes. The results suggest that synthe-sized seismograms coincide well with observed seismic waveforms of the 1995 Hyogo-ken Nanbu earthquake. The simulated Green’s function method is very useful for prediction of the strong ground mo-tion in region without observed seismic waveforms. The present technique spreads application field of the empirical Green’s function method.
