Aimed at the poor performance of conventional geophones in exploration for deeper and complex targets, we present the principle and theoretical design of a new geophone based on the optical fiber Bragg grating （FBG）...Aimed at the poor performance of conventional geophones in exploration for deeper and complex targets, we present the principle and theoretical design of a new geophone based on the optical fiber Bragg grating （FBG） sensing technology. The important parameters such as response functions are calculated theoretically. Because of the advantages of FBG sensing technology, the new FBG geophone has a high dynamic range of 94dB at （10-200 Hz）. This new generation of geophones will have wide use in seismic prospecting due to its higher sensitivity, lighter weight, and lower cost.展开更多
The dynamic characteristics of ground soil using micro-tremor observation in Asia (Zushi and Ogasawara (Japan), Xi’an (China), Manila (Philippines), and Gujarat (India)) are studied. Ground micro-tremor signals were ...The dynamic characteristics of ground soil using micro-tremor observation in Asia (Zushi and Ogasawara (Japan), Xi’an (China), Manila (Philippines), and Gujarat (India)) are studied. Ground micro-tremor signals were observed and analyzed by fast Fourier transform method (FFT). The response of ground soil to frequency of ground micro-tremor is revealed, and func- tions with frequency-dependence and frequency-selection of micro-tremor for different foundation soil strata are also researched. The horizontal to vertical spectral ratio (H/V, Nakamura technique) of micro-tremor observed at the surface ground was used to evaluate the site’s predominant period. This paper also discusses the application of micro-tremor on site safety evaluation, and gives the observed calculation results obtained at multiple points. The experimental foundation and the deduction process of the method are described in detail. Some problems of the method are pointed out. Potential use of the technique’s good expandable nature makes it a useable means for preventing and reducing disaster’s harmful effects.展开更多
High-frequency seismic data components can be seriously attenuated during seismic wave propagation in unconsolidated （low-velocity） layers, resulting in reduced seismic resolution and signal-to-noise （S/N） ratio. ...High-frequency seismic data components can be seriously attenuated during seismic wave propagation in unconsolidated （low-velocity） layers, resulting in reduced seismic resolution and signal-to-noise （S/N） ratio. In this paper, first, based on Wiener filter theory, inverse filter calculations for near-surface absorption attenuation compensation were accomplished by analysis of the direct wave spectral components from different distances near the surface. The direct waves were generated by detonators in uphole shots and were acquired by receivers on the surface. The spatially varying inverse filters were designed to compensate for the frequency attenuation of 3D pre-stack CRG （common receiver-gather） data. After applying the filter to CRG data, the high frequency components were compensated with the low frequencies maintained. The seismic resolution and S/N ratio are enhanced and match better with synthetic seismograms and better meet the needs of geological interpretation.展开更多
A new seismic ray-tracing method is put forward based on parabolic travel-time interpolation(PTI) method, which is more accurate than the linear travel-time interpolation (LTI) method. Both PTI method and LTI method a...A new seismic ray-tracing method is put forward based on parabolic travel-time interpolation(PTI) method, which is more accurate than the linear travel-time interpolation (LTI) method. Both PTI method and LTI method are used to compute seismic travel-time and ray-path in a 2-D grid cell model. Firstly, some basic concepts are introduced. The calculations of travel-time and ray-path are carried out only at cell boundaries. So, the ray-path is always straight in the same cells with uniform velocity. Two steps are applied in PTI and LTI method, step 1 computes travel-time and step 2 traces ray-path. Then, the derivation of LTI formulas is described. Because of the presence of refraction wave in shot cell, the formula aiming at shot cell is also derived. Finally, PTI method is presented. The calculation of PTI method is more complex than that of LTI method, but the error is limited. The results of numerical model show that PTI method can trace ray-path more accurately and efficiently than LTI method does.展开更多
Based on Fermat’s principle, two-point ray tracing method was studied in three-dimensional structure. By means of first order Taylor’s incomplete series expansion (i.e. no expansion to the length of the ray), a symm...Based on Fermat’s principle, two-point ray tracing method was studied in three-dimensional structure. By means of first order Taylor’s incomplete series expansion (i.e. no expansion to the length of the ray), a symmetry block tridiagonal matrix equation set was deduced. Further, the positive definiteness of coefficient matrix was discussed, and the positive definiteness was accurately proved in a mathematical way. It assured that the algorithm was well-posed. Associated with iterative method, the solution to ray tracing can be got through step-by-step linearized iteration of the nonlinear problem. An algorithm of the whole path iterative ray tracing method in three-dimensional velocity structure was obtained. This method shows a clear and simple as well as explicit computation formula, which makes ray tracing computation easily applicable in practice. The correction vector is obtained through finding the solution to the positive definite block tridiagonal equation set, which ensures the method is robust convergence. This study offers a new kind of feasible and efficient ray tracing method for three dimensional seismic migration and tomography. Meanwhile, it also provides the prerequisite guarantee to design a fast algorithm.展开更多
Although earthquake lights have been known since ancient times,it has not been easy to study them.It was not until the 60s that the first photographs of them were taken.During the Peruvian earthquake in 2007,it was po...Although earthquake lights have been known since ancient times,it has not been easy to study them.It was not until the 60s that the first photographs of them were taken.During the Peruvian earthquake in 2007,it was possible to obtain the fi rst fi lm recording on earthquake lights.Likewise,during the earthquakes in Ecuador in 2016 and in Mexico in 2017,two fi lms of the earthquake lights were recorded.These fi lm recordings have helped in the study of earthquake lights,both for their objectivity and for their informational content.Several causal mechanisms have been proposed to explain earthquake lights:piezoelectricity,radon emanation,fluid diffusion,friction-vaporization,positive holes and dipole currents,among others.In this work a time difference correlation between earthquake lights and seismic ground accelerations was found and we use both seismic data and fi lm recordings of earthquake lights to explain its origin.In the discussion section it is suggested that fracturing of rocks manifest itself to some extent in the form of static electricity producing earthquake lights through induction The induction model proposed is new and it can explain the formation of EQL,even if the earth’s crust has layers of large electrical resistivity.The model also explains the formation of seismic lights without the need for special conditions on the earth’s surface or in the atmosphere.A better understanding of the earthquake lights generation process can improve our understanding of seismicity and help in the prediction of earthquakes.展开更多
By applying the empirical Green＇s function method, we first proved that the source spectra in the middle eastern area of North Tianshan follow the ω^-2 Brune model, and then, on this basis, obtained the source spect...By applying the empirical Green＇s function method, we first proved that the source spectra in the middle eastern area of North Tianshan follow the ω^-2 Brune model, and then, on this basis, obtained the source spectra for 105 earthquakes of ML 2.5 - 5.7 in the studied area upon removing the noise, instrument, propagation and site effects in the observational seismic S waveforms. Finally, we determined the source parameters such as the seismic moment, stress drop, source dimension, etc., based on the Brune model. The results show that there is a good linear correlation between scalar seismic moment, M0, and ML, which is Log10 M0 = 1.10 ML ＋ 17.20; The source radii range between 100- 1500m; the stress drop, 1-16MPa with a predominant range of 1 - 10MPa. The relationship of the seismic moment versus corner frequencies indicates that there may exist two source scaling, that is, when M0 〈 4×10^21 dyne. cm （equal to an ML4.0 event）, stress drop is weakly dependent on M0： whereas, when M0 〉 4×10^21 dyne·cm, stress drop is independent of M0.展开更多
Although the CTBT (Comprehensive Nuclear Test Ban Treaty) was passed in 1996, it is still necessary to develop new and highly efficient methods (Wu Zhongliang, Chen Yuntai, et al., 1993; Xu Shaoxie, et al.1994; Richar...Although the CTBT (Comprehensive Nuclear Test Ban Treaty) was passed in 1996, it is still necessary to develop new and highly efficient methods (Wu Zhongliang, Chen Yuntai, et al., 1993; Xu Shaoxie, et al.1994; Richard L. Garwin, 1994) to monitor possible events. Many discrimination criteria (Xu Shaoxie, et al.,1994; Institute of Geophysics, Chinese Academy of Sciences, 1976; Richard L. Garwin, 1994) have been put forward since the 1950s. The results show that each of the existing criteria has its own limitation, but the seismological method is an important and efficient method in the discrimination between nuclear explosion and earthquake. Especially in recent years, because of the little and little equivalent as well as the increasing hiding steps used in the test, a number of more efficient seismological methods have been worked out. In this paper, a new discrimination method, the Wavelet Packet Component Ratio (WPCR) method, is put forward. This method makes full use of the difference in variation with time between the spectra of nuclear explosions and earthquakes. Its discrimination efficiency is rather high.展开更多
基金sponsored by the National 863 Program(Grant No.2006AA06Z207&2006AA06Z213)the National Natural Science Foundation of China(Grant No.50674098)the National 973 Program(Grant No.2007CB209601)
文摘Aimed at the poor performance of conventional geophones in exploration for deeper and complex targets, we present the principle and theoretical design of a new geophone based on the optical fiber Bragg grating （FBG） sensing technology. The important parameters such as response functions are calculated theoretically. Because of the advantages of FBG sensing technology, the new FBG geophone has a high dynamic range of 94dB at （10-200 Hz）. This new generation of geophones will have wide use in seismic prospecting due to its higher sensitivity, lighter weight, and lower cost.
文摘The dynamic characteristics of ground soil using micro-tremor observation in Asia (Zushi and Ogasawara (Japan), Xi’an (China), Manila (Philippines), and Gujarat (India)) are studied. Ground micro-tremor signals were observed and analyzed by fast Fourier transform method (FFT). The response of ground soil to frequency of ground micro-tremor is revealed, and func- tions with frequency-dependence and frequency-selection of micro-tremor for different foundation soil strata are also researched. The horizontal to vertical spectral ratio (H/V, Nakamura technique) of micro-tremor observed at the surface ground was used to evaluate the site’s predominant period. This paper also discusses the application of micro-tremor on site safety evaluation, and gives the observed calculation results obtained at multiple points. The experimental foundation and the deduction process of the method are described in detail. Some problems of the method are pointed out. Potential use of the technique’s good expandable nature makes it a useable means for preventing and reducing disaster’s harmful effects.
基金supported by China Petroleum Technology Innovation Fund Project(Grant No.0610740122)
文摘High-frequency seismic data components can be seriously attenuated during seismic wave propagation in unconsolidated （low-velocity） layers, resulting in reduced seismic resolution and signal-to-noise （S/N） ratio. In this paper, first, based on Wiener filter theory, inverse filter calculations for near-surface absorption attenuation compensation were accomplished by analysis of the direct wave spectral components from different distances near the surface. The direct waves were generated by detonators in uphole shots and were acquired by receivers on the surface. The spatially varying inverse filters were designed to compensate for the frequency attenuation of 3D pre-stack CRG （common receiver-gather） data. After applying the filter to CRG data, the high frequency components were compensated with the low frequencies maintained. The seismic resolution and S/N ratio are enhanced and match better with synthetic seismograms and better meet the needs of geological interpretation.
文摘A new seismic ray-tracing method is put forward based on parabolic travel-time interpolation(PTI) method, which is more accurate than the linear travel-time interpolation (LTI) method. Both PTI method and LTI method are used to compute seismic travel-time and ray-path in a 2-D grid cell model. Firstly, some basic concepts are introduced. The calculations of travel-time and ray-path are carried out only at cell boundaries. So, the ray-path is always straight in the same cells with uniform velocity. Two steps are applied in PTI and LTI method, step 1 computes travel-time and step 2 traces ray-path. Then, the derivation of LTI formulas is described. Because of the presence of refraction wave in shot cell, the formula aiming at shot cell is also derived. Finally, PTI method is presented. The calculation of PTI method is more complex than that of LTI method, but the error is limited. The results of numerical model show that PTI method can trace ray-path more accurately and efficiently than LTI method does.
基金Project(40674071) supported by the National Natural Science Foundation of Chinaproject(KFAS2002-2003) supported by Korea Foundation for Advanced Studies
文摘Based on Fermat’s principle, two-point ray tracing method was studied in three-dimensional structure. By means of first order Taylor’s incomplete series expansion (i.e. no expansion to the length of the ray), a symmetry block tridiagonal matrix equation set was deduced. Further, the positive definiteness of coefficient matrix was discussed, and the positive definiteness was accurately proved in a mathematical way. It assured that the algorithm was well-posed. Associated with iterative method, the solution to ray tracing can be got through step-by-step linearized iteration of the nonlinear problem. An algorithm of the whole path iterative ray tracing method in three-dimensional velocity structure was obtained. This method shows a clear and simple as well as explicit computation formula, which makes ray tracing computation easily applicable in practice. The correction vector is obtained through finding the solution to the positive definite block tridiagonal equation set, which ensures the method is robust convergence. This study offers a new kind of feasible and efficient ray tracing method for three dimensional seismic migration and tomography. Meanwhile, it also provides the prerequisite guarantee to design a fast algorithm.
文摘Although earthquake lights have been known since ancient times,it has not been easy to study them.It was not until the 60s that the first photographs of them were taken.During the Peruvian earthquake in 2007,it was possible to obtain the fi rst fi lm recording on earthquake lights.Likewise,during the earthquakes in Ecuador in 2016 and in Mexico in 2017,two fi lms of the earthquake lights were recorded.These fi lm recordings have helped in the study of earthquake lights,both for their objectivity and for their informational content.Several causal mechanisms have been proposed to explain earthquake lights:piezoelectricity,radon emanation,fluid diffusion,friction-vaporization,positive holes and dipole currents,among others.In this work a time difference correlation between earthquake lights and seismic ground accelerations was found and we use both seismic data and fi lm recordings of earthquake lights to explain its origin.In the discussion section it is suggested that fracturing of rocks manifest itself to some extent in the form of static electricity producing earthquake lights through induction The induction model proposed is new and it can explain the formation of EQL,even if the earth’s crust has layers of large electrical resistivity.The model also explains the formation of seismic lights without the need for special conditions on the earth’s surface or in the atmosphere.A better understanding of the earthquake lights generation process can improve our understanding of seismicity and help in the prediction of earthquakes.
文摘By applying the empirical Green＇s function method, we first proved that the source spectra in the middle eastern area of North Tianshan follow the ω^-2 Brune model, and then, on this basis, obtained the source spectra for 105 earthquakes of ML 2.5 - 5.7 in the studied area upon removing the noise, instrument, propagation and site effects in the observational seismic S waveforms. Finally, we determined the source parameters such as the seismic moment, stress drop, source dimension, etc., based on the Brune model. The results show that there is a good linear correlation between scalar seismic moment, M0, and ML, which is Log10 M0 = 1.10 ML ＋ 17.20; The source radii range between 100- 1500m; the stress drop, 1-16MPa with a predominant range of 1 - 10MPa. The relationship of the seismic moment versus corner frequencies indicates that there may exist two source scaling, that is, when M0 〈 4×10^21 dyne. cm （equal to an ML4.0 event）, stress drop is weakly dependent on M0： whereas, when M0 〉 4×10^21 dyne·cm, stress drop is independent of M0.
文摘Although the CTBT (Comprehensive Nuclear Test Ban Treaty) was passed in 1996, it is still necessary to develop new and highly efficient methods (Wu Zhongliang, Chen Yuntai, et al., 1993; Xu Shaoxie, et al.1994; Richard L. Garwin, 1994) to monitor possible events. Many discrimination criteria (Xu Shaoxie, et al.,1994; Institute of Geophysics, Chinese Academy of Sciences, 1976; Richard L. Garwin, 1994) have been put forward since the 1950s. The results show that each of the existing criteria has its own limitation, but the seismological method is an important and efficient method in the discrimination between nuclear explosion and earthquake. Especially in recent years, because of the little and little equivalent as well as the increasing hiding steps used in the test, a number of more efficient seismological methods have been worked out. In this paper, a new discrimination method, the Wavelet Packet Component Ratio (WPCR) method, is put forward. This method makes full use of the difference in variation with time between the spectra of nuclear explosions and earthquakes. Its discrimination efficiency is rather high.
