Reflected wave seismology has the following defects:the acquisition design is based on the assumption of layered media,the signal processing suppresses weak signals such as diffracted wave and scattered wave,and the s...Reflected wave seismology has the following defects:the acquisition design is based on the assumption of layered media,the signal processing suppresses weak signals such as diffracted wave and scattered wave,and the seismic wave band after the image processing is narrow.They limit the full utilization of broadband raw data.The concept of full wave seismic exploration is redefined based on the idea of balanced utilization of reflected wave,diffracted wave and scattered wave information,its characteristics and adaptive conditions are clarified.A set of key technologies suitable for full wave seismic exploration are put forward.During seismic acquisition period,it is necessary to adopt multi geometry,i.e.embed small bin,small offset and small channel interval data in conventional geometry.By discretizing of common midpoint(CMP)gathers,small offset with high coverage,the weak signals such as diffracted wave and scattered wave in the raw seismic data can be enhanced.During seismic processing,the signal and noise in the original seismic data need to be redefined at first.The effective signals of seismic data are enhanced through merging of multi-geometry data.By means of differential application of data with different bin sizes and different arrangement modes,different regimes of seismic waves can be effectively decomposed and imaged separately.During seismic interpretation stage,making the most of the full wave seismic data,and adopting well-seismic calibration on multi-scale and multi-dimension,the seismic attributes in multi-regimes and multi-domains are interpreted to reveal interior information of complex lithology bodies and improve the lateral resolution of non-layered reservoirs.展开更多
We apply the spectral-element method（SEM）,a high-order finite-element method（FEM） to simulate seismic wave propagation in complex media for exploration and geotechnical problems. The SEM accurately treats geometri...We apply the spectral-element method（SEM）,a high-order finite-element method（FEM） to simulate seismic wave propagation in complex media for exploration and geotechnical problems. The SEM accurately treats geometrical complexities through its flexible FEM mesh and accurately interpolates wavefields through high-order Lagrange polynomials. It has been a numerical solver used extensively in earthquake seismology. We demonstrate the applicability of SEM for selected 2D exploration and geotechnical velocity models with an open-source SEM software package SPECFEM2D. The first scenario involves a marine survey for a salt dome with the presence of major internal discontinuities,and the second example simulates seismic wave propagation for an open-pit mine with complex surface topography. Wavefield snapshots,synthetic seismograms,and peak particle velocity maps are presented to illustrate the promising use of SEM for industrial problems.展开更多
Pingtong Town is located on the fractured zone of the Wenchuan 8.0 earthquake, and is seriously damaged by the earthquake. Our observation line is centered at an earthquake exploration trench across the fractured zone...Pingtong Town is located on the fractured zone of the Wenchuan 8.0 earthquake, and is seriously damaged by the earthquake. Our observation line is centered at an earthquake exploration trench across the fractured zone in the NW-SE direction, and is about 400 m long. The results reveal trapped waves in the rup- tured fault zone of the earthquake, and indicate a great difference in physical property between the media inside and outside the fault zone. The predominant frequency of the fault-zone trapped waves is about 3 -4 Hz. The wave amplitudes are larger near the exploration trench. The width of the fault zone in the crust at this location is estimated to be 200 m. In some records, the waveforms and the arrival times of S waves are quite different between the two sides of the trench. The place of change coincides with the boundary of uplift at the surface.展开更多
Exploration and research of fault activities are the fundamentals of earthquake prediction and prevention and disaster reduction. In order to determine the location, characteristics and activities of the Zhengzhon-Lao...Exploration and research of fault activities are the fundamentals of earthquake prediction and prevention and disaster reduction. In order to determine the location, characteristics and activities of the Zhengzhon-Laoyachen fault, shallow seismic prospecting with different exploration depth across the Laoyachen fault was carried out in the northern suburbs of Zhengzhou city in 2006. The images of the subterranean structure and tectonics at depths of 30m- 6000m have been available by applying the combined methods of explosive seismic sources and vibrator seismic sources, as well as the combination of diverse observation systems with different parameters. The outcome indicates that the Laoyachen fault is a normal fault running NW and dipping NE, which offsets stratums ahead of Neogene （N）. However, no fault displacements are found in the interior stratums of Q ＋ N.展开更多
The main problems in seismic attribute technology are the redundancy of data and the uncertainty of attributes, and these problems become much more serious in multi-wave seismic exploration. Data redundancy will incre...The main problems in seismic attribute technology are the redundancy of data and the uncertainty of attributes, and these problems become much more serious in multi-wave seismic exploration. Data redundancy will increase the burden on interpreters, occupy large computer memory, take much more computing time, conceal the effective information, and especially cause the ＂curse of dimension＂. Uncertainty of attributes will reduce the accuracy of rebuilding the relationship between attributes and geological significance. In order to solve these problems, we study methods of principal component analysis （PCA）, independent component analysis （ICA） for attribute optimization and support vector machine （SVM） for reservoir prediction. We propose a flow chart of multi-wave seismic attribute process and further apply it to multi-wave seismic reservoir prediction. The processing results of real seismic data demonstrate that reservoir prediction based on combination of PP- and PS-wave attributes, compared with that based on traditional PP-wave attributes, can improve the prediction accuracy.展开更多
Seismological Bureau of Sichuan Province, Chengdu 610041, China2) Center for Analysis and Prediction, State Seismological Bureau, Beijing 100036, China3) Observation Center for Prediction of Earthquakes and Volcanic E...Seismological Bureau of Sichuan Province, Chengdu 610041, China2) Center for Analysis and Prediction, State Seismological Bureau, Beijing 100036, China3) Observation Center for Prediction of Earthquakes and Volcanic Eruptions, Faculty of Sciences, Tohoku University, Sendai 98077, Japan展开更多
Because ambient seismic noise provides estimated Green’s function （EGF） between two sites with high accuracy, Rayleigh wave propagation along the path connecting the two sites is well resolved. Therefore, earthquak...Because ambient seismic noise provides estimated Green’s function （EGF） between two sites with high accuracy, Rayleigh wave propagation along the path connecting the two sites is well resolved. Therefore, earthquakes which are close to one seismic station can be well located with calibration extracting from EGF. We test two algorithms in locating the 1998 Zhangbei earthquake, one algorithm is waveform-based, and the other is traveltime-based. We first compute EGF between station ZHB （a station about 40 km away from the epicenter） and five IC/IRIS stations. With the waveform-based approach, we calculate 1D synthetic single-force Green’s functions between ZHB and other four stations, and obtain traveltime corrections by correlating synthetic Green’s functions with EGFs in period band of 10–30 s. Then we locate the earthquake by minimizing the differential travel times between observed earthquake waveform and the 1D synthetic earthquake waveforms computed with focal mechanism provided by Global CMT after traveltime correction from EGFs. This waveform-based approach yields a location which error is about 13 km away from the location observed with InSAR. With the traveltime-based approach, we begin with measuring group velocity from EGFs as well as group arrival time on observed earthquake waveforms, and then locate the earthquake by minimizing the difference between observed group arrival time and arrival time measured on EGFs. This traveltime-based approach yields accuracy of 3 km, Therefore it is feasible to achieve GT5 （ground truth location with accuracy 5 km） with ambient seismic noises. The less accuracy of the waveform-based approach was mainly caused by uncertainty of focal mechanism.展开更多
Forecasting subtle traps by sequence stratigraphy and 3D seismic data is a sensitive topic in hydrocarbon exploration. Research on subtle traps by geophysical data is the most popular and difficult. Based on the suffi...Forecasting subtle traps by sequence stratigraphy and 3D seismic data is a sensitive topic in hydrocarbon exploration. Research on subtle traps by geophysical data is the most popular and difficult. Based on the sufficiently drilling data, log data, core data and 3D seismic data, sediment sequence of Qikou depression, Huanghua basin was partitioned by using sequence stratigraphy theory. Each sediment sequence system mode was built. Sediment faces of subtle traps were pointed out. Dominating factors forming subtle traps were analyzed. Sandstone seismic rock physics and its response were studied in Tertiary System. Sandstone geophysical response and elastic modulus vary laws with pressure, temperature, porosity, depth were built. Experimental result and practice shows that it is possible using seismic information forecasting subtle traps. Integrated using geology, log, drilling data, special seismic processing technique, interpretation technique, high precision horizon calibration technique, 3D seismic visualizing interpretation, seismic coherence analysis, attribute analysis, logging-constrained inversion, time frequency analysis, subtle trapsobject is identified and interpreted. Finally, advantage object of subtle trap in this area was determined. Bottomland sand stratigraphic and lithologic reservoirs in Qinan slope zone have been founded by means of high resolution 3D seismic data field technique, high resolution 3D seismic data processing technique and seismic wave impendence inversion technique.展开更多
The influence of the dispersion and uncertainty of the dynamic shear wave velocity and Poisson＇s ratio of soil in a hard rock site was investigated on the seismic response of reactor building structure. The analysis ...The influence of the dispersion and uncertainty of the dynamic shear wave velocity and Poisson＇s ratio of soil in a hard rock site was investigated on the seismic response of reactor building structure. The analysis is performed by considering the soil-structure interaction effects and based on the model of the reactor building in a typical pressurized water reactor nuclear power plant （NPP）. The numerical results show that for the typical floor selected, while the relative increment ratio of the dynamic shear wave velocity varies from -30% to 30% compared to the basis of 1 930 m/s, the relative variation of the horizontal response spectra peak value lies in the scope of ±10% for the internal structure, and the relative variation of the frequency corresponding to the spectra peak is 0.0% in most cases. The relative variation of the vertical response spectra peak value lies in the scope of - 10% to 22%, and the relative variation of the frequency corresponding to the Spectra peak lies in the scope of - 22% to 4%. The analysis indicates that the dynamic shear wave velocity and the Poisson＇s ratio of the rock would affect the seismic response of structure and the soil-structure interaction effects should be considered in seismic analysis and design of NPP even for a hard rock site.展开更多
Based on the inversion method of 2D velocity structure and interface, the crustal velocity structures of P-wave and S-wave along the profile L1 are determined simultaneously with deep seismic sounding data in Changbai...Based on the inversion method of 2D velocity structure and interface, the crustal velocity structures of P-wave and S-wave along the profile L1 are determined simultaneously with deep seismic sounding data in Changbaishan Tianchi volcanic region, and then its Poisson's ratio is obtained. Calculated results show that this technique overcomes some defects of traditional forward calculation method, and it is also very effective to determine Poisson's ratio distribution of deep seismic sounding profile, especially useful for study on volcanic magma and crustal fault zone. Study result indicates that there is an abnormally high Poisson's ratio body that is about 30 km wide and 12 km high in the low velocity region under Tianchi crater. Its value of Poisson's ratio is 8% higher than that of surrounding medium and it should be the magma chamber formed from melted rock with high temperature. There is a high Poisson's ratio zone ranging from magma chamber to the top of crust, which may be the uprise passage of hot substance. The lower part with high Poisson's ratio, which stretches downward to Moho, is possibly the extrusion way of hot substance from the uppermost mantle. The conclusions above are consistent with the study results of both tomographic determination of 3D crustal structure and magnetotelluric survey in this region.展开更多
An inversion method was applied to crustal earthquakes dataset to find S-wave attenuation characteristics beneath the Eastern Tohoku region of Japan. Accelerograms from 85 shallow crustal earthquakes up to 25 km depth...An inversion method was applied to crustal earthquakes dataset to find S-wave attenuation characteristics beneath the Eastern Tohoku region of Japan. Accelerograms from 85 shallow crustal earthquakes up to 25 km depth and magnitude range between 3.5 and 5.5 were analyzed to estimate the seismic quality factor Qs. A homogeneous attenuation model Qs for the wave propagation path was evaluated from spectral amplitudes, at 24 different frequencies between 0.5 and 20 Hz by using generalized inversion technique. To do this, non-parametric attenuation functions were calculated to observe spectral amplitude decay with hypocentral distance. Then, these functions were parameterized to estimate Qs. It was found that in Eastern Tohoku region, the Qs frequency dependence can be approximated with the function 33 f 1.22 within a frequency range between 0.5 and 20 Hz. However, the frequency dependence of Qs in the frequency range between 0.5 and 6 Hz is best approximated by Qs （f） = 36 f 0.94 showing relatively weaker frequency dependence as compared to the relation Qs （f） = 6 f^ 2.09 for the frequency range between 6 and 15 Hz. These results could be used to estimate source and site parameters for seismic hazard assessment in the region.展开更多
基金Supported by the Sinopec Ministry of Science and Technology Project(P21038-3)。
文摘Reflected wave seismology has the following defects:the acquisition design is based on the assumption of layered media,the signal processing suppresses weak signals such as diffracted wave and scattered wave,and the seismic wave band after the image processing is narrow.They limit the full utilization of broadband raw data.The concept of full wave seismic exploration is redefined based on the idea of balanced utilization of reflected wave,diffracted wave and scattered wave information,its characteristics and adaptive conditions are clarified.A set of key technologies suitable for full wave seismic exploration are put forward.During seismic acquisition period,it is necessary to adopt multi geometry,i.e.embed small bin,small offset and small channel interval data in conventional geometry.By discretizing of common midpoint(CMP)gathers,small offset with high coverage,the weak signals such as diffracted wave and scattered wave in the raw seismic data can be enhanced.During seismic processing,the signal and noise in the original seismic data need to be redefined at first.The effective signals of seismic data are enhanced through merging of multi-geometry data.By means of differential application of data with different bin sizes and different arrangement modes,different regimes of seismic waves can be effectively decomposed and imaged separately.During seismic interpretation stage,making the most of the full wave seismic data,and adopting well-seismic calibration on multi-scale and multi-dimension,the seismic attributes in multi-regimes and multi-domains are interpreted to reveal interior information of complex lithology bodies and improve the lateral resolution of non-layered reservoirs.
基金supported by the Natural Sciences and Engineering Research Council of Canada (NSERC)Center for Excellence in Mining Innovations (CEMI,through SUMIT project)+2 种基金Computations for this study were performed on hardwares purchased through the combined funding of Canada Foundation for Innovation (CFI)Ontario Research Fund (ORF)University of Toronto Startup Fund
文摘We apply the spectral-element method（SEM）,a high-order finite-element method（FEM） to simulate seismic wave propagation in complex media for exploration and geotechnical problems. The SEM accurately treats geometrical complexities through its flexible FEM mesh and accurately interpolates wavefields through high-order Lagrange polynomials. It has been a numerical solver used extensively in earthquake seismology. We demonstrate the applicability of SEM for selected 2D exploration and geotechnical velocity models with an open-source SEM software package SPECFEM2D. The first scenario involves a marine survey for a salt dome with the presence of major internal discontinuities,and the second example simulates seismic wave propagation for an open-pit mine with complex surface topography. Wavefield snapshots,synthetic seismograms,and peak particle velocity maps are presented to illustrate the promising use of SEM for industrial problems.
基金supported by the Natural Science Fundation of China(40774043,40674043,90814001)
文摘Pingtong Town is located on the fractured zone of the Wenchuan 8.0 earthquake, and is seriously damaged by the earthquake. Our observation line is centered at an earthquake exploration trench across the fractured zone in the NW-SE direction, and is about 400 m long. The results reveal trapped waves in the rup- tured fault zone of the earthquake, and indicate a great difference in physical property between the media inside and outside the fault zone. The predominant frequency of the fault-zone trapped waves is about 3 -4 Hz. The wave amplitudes are larger near the exploration trench. The width of the fault zone in the crust at this location is estimated to be 200 m. In some records, the waveforms and the arrival times of S waves are quite different between the two sides of the trench. The place of change coincides with the boundary of uplift at the surface.
基金sponsored by the State Development and Planning Commission(200197)the Key Projects of thetenth"Five-year Plan" of Henan People’s Government
文摘Exploration and research of fault activities are the fundamentals of earthquake prediction and prevention and disaster reduction. In order to determine the location, characteristics and activities of the Zhengzhon-Laoyachen fault, shallow seismic prospecting with different exploration depth across the Laoyachen fault was carried out in the northern suburbs of Zhengzhou city in 2006. The images of the subterranean structure and tectonics at depths of 30m- 6000m have been available by applying the combined methods of explosive seismic sources and vibrator seismic sources, as well as the combination of diverse observation systems with different parameters. The outcome indicates that the Laoyachen fault is a normal fault running NW and dipping NE, which offsets stratums ahead of Neogene （N）. However, no fault displacements are found in the interior stratums of Q ＋ N.
基金supported by China Important National Science & Technology Specific Projects (No.2011ZX05019-008)National Natural Science Foundation of China (No.40839901)
文摘The main problems in seismic attribute technology are the redundancy of data and the uncertainty of attributes, and these problems become much more serious in multi-wave seismic exploration. Data redundancy will increase the burden on interpreters, occupy large computer memory, take much more computing time, conceal the effective information, and especially cause the ＂curse of dimension＂. Uncertainty of attributes will reduce the accuracy of rebuilding the relationship between attributes and geological significance. In order to solve these problems, we study methods of principal component analysis （PCA）, independent component analysis （ICA） for attribute optimization and support vector machine （SVM） for reservoir prediction. We propose a flow chart of multi-wave seismic attribute process and further apply it to multi-wave seismic reservoir prediction. The processing results of real seismic data demonstrate that reservoir prediction based on combination of PP- and PS-wave attributes, compared with that based on traditional PP-wave attributes, can improve the prediction accuracy.
文摘Seismological Bureau of Sichuan Province, Chengdu 610041, China2) Center for Analysis and Prediction, State Seismological Bureau, Beijing 100036, China3) Observation Center for Prediction of Earthquakes and Volcanic Eruptions, Faculty of Sciences, Tohoku University, Sendai 98077, Japan
基金supported by Chinese Acadmy of Sciences Fund(No.KCZX-YW-116-1)Joint Seismological Science Fundation of China (Nos.20080878 and 200708035)
文摘Because ambient seismic noise provides estimated Green’s function （EGF） between two sites with high accuracy, Rayleigh wave propagation along the path connecting the two sites is well resolved. Therefore, earthquakes which are close to one seismic station can be well located with calibration extracting from EGF. We test two algorithms in locating the 1998 Zhangbei earthquake, one algorithm is waveform-based, and the other is traveltime-based. We first compute EGF between station ZHB （a station about 40 km away from the epicenter） and five IC/IRIS stations. With the waveform-based approach, we calculate 1D synthetic single-force Green’s functions between ZHB and other four stations, and obtain traveltime corrections by correlating synthetic Green’s functions with EGFs in period band of 10–30 s. Then we locate the earthquake by minimizing the differential travel times between observed earthquake waveform and the 1D synthetic earthquake waveforms computed with focal mechanism provided by Global CMT after traveltime correction from EGFs. This waveform-based approach yields a location which error is about 13 km away from the location observed with InSAR. With the traveltime-based approach, we begin with measuring group velocity from EGFs as well as group arrival time on observed earthquake waveforms, and then locate the earthquake by minimizing the difference between observed group arrival time and arrival time measured on EGFs. This traveltime-based approach yields accuracy of 3 km, Therefore it is feasible to achieve GT5 （ground truth location with accuracy 5 km） with ambient seismic noises. The less accuracy of the waveform-based approach was mainly caused by uncertainty of focal mechanism.
基金Project(2003034470) supported by the Postdoctoral Science Foundation of China project supported by the Postdoctoral Science Foundation of Central South University
文摘Forecasting subtle traps by sequence stratigraphy and 3D seismic data is a sensitive topic in hydrocarbon exploration. Research on subtle traps by geophysical data is the most popular and difficult. Based on the sufficiently drilling data, log data, core data and 3D seismic data, sediment sequence of Qikou depression, Huanghua basin was partitioned by using sequence stratigraphy theory. Each sediment sequence system mode was built. Sediment faces of subtle traps were pointed out. Dominating factors forming subtle traps were analyzed. Sandstone seismic rock physics and its response were studied in Tertiary System. Sandstone geophysical response and elastic modulus vary laws with pressure, temperature, porosity, depth were built. Experimental result and practice shows that it is possible using seismic information forecasting subtle traps. Integrated using geology, log, drilling data, special seismic processing technique, interpretation technique, high precision horizon calibration technique, 3D seismic visualizing interpretation, seismic coherence analysis, attribute analysis, logging-constrained inversion, time frequency analysis, subtle trapsobject is identified and interpreted. Finally, advantage object of subtle trap in this area was determined. Bottomland sand stratigraphic and lithologic reservoirs in Qinan slope zone have been founded by means of high resolution 3D seismic data field technique, high resolution 3D seismic data processing technique and seismic wave impendence inversion technique.
基金SUPPORTED BY NATIONAL NATURAL SCIENCE FOUNDATION FOR DISTINGUISHED YOUNG SCHOLARS OF CHINA (NO. 50425824).
文摘The influence of the dispersion and uncertainty of the dynamic shear wave velocity and Poisson＇s ratio of soil in a hard rock site was investigated on the seismic response of reactor building structure. The analysis is performed by considering the soil-structure interaction effects and based on the model of the reactor building in a typical pressurized water reactor nuclear power plant （NPP）. The numerical results show that for the typical floor selected, while the relative increment ratio of the dynamic shear wave velocity varies from -30% to 30% compared to the basis of 1 930 m/s, the relative variation of the horizontal response spectra peak value lies in the scope of ±10% for the internal structure, and the relative variation of the frequency corresponding to the spectra peak is 0.0% in most cases. The relative variation of the vertical response spectra peak value lies in the scope of - 10% to 22%, and the relative variation of the frequency corresponding to the Spectra peak lies in the scope of - 22% to 4%. The analysis indicates that the dynamic shear wave velocity and the Poisson＇s ratio of the rock would affect the seismic response of structure and the soil-structure interaction effects should be considered in seismic analysis and design of NPP even for a hard rock site.
基金Key Project from China Earthquake Administration and the Project (95-11-02-01) from Ministry of Science and Technology (2001DIA10003).
文摘Based on the inversion method of 2D velocity structure and interface, the crustal velocity structures of P-wave and S-wave along the profile L1 are determined simultaneously with deep seismic sounding data in Changbaishan Tianchi volcanic region, and then its Poisson's ratio is obtained. Calculated results show that this technique overcomes some defects of traditional forward calculation method, and it is also very effective to determine Poisson's ratio distribution of deep seismic sounding profile, especially useful for study on volcanic magma and crustal fault zone. Study result indicates that there is an abnormally high Poisson's ratio body that is about 30 km wide and 12 km high in the low velocity region under Tianchi crater. Its value of Poisson's ratio is 8% higher than that of surrounding medium and it should be the magma chamber formed from melted rock with high temperature. There is a high Poisson's ratio zone ranging from magma chamber to the top of crust, which may be the uprise passage of hot substance. The lower part with high Poisson's ratio, which stretches downward to Moho, is possibly the extrusion way of hot substance from the uppermost mantle. The conclusions above are consistent with the study results of both tomographic determination of 3D crustal structure and magnetotelluric survey in this region.
基金a part of author’s M.Sc Research under the project:‘‘Strengthening of Earthquake Engineering Center’’,funded by Higher Education Commission,Government of Pakistan
文摘An inversion method was applied to crustal earthquakes dataset to find S-wave attenuation characteristics beneath the Eastern Tohoku region of Japan. Accelerograms from 85 shallow crustal earthquakes up to 25 km depth and magnitude range between 3.5 and 5.5 were analyzed to estimate the seismic quality factor Qs. A homogeneous attenuation model Qs for the wave propagation path was evaluated from spectral amplitudes, at 24 different frequencies between 0.5 and 20 Hz by using generalized inversion technique. To do this, non-parametric attenuation functions were calculated to observe spectral amplitude decay with hypocentral distance. Then, these functions were parameterized to estimate Qs. It was found that in Eastern Tohoku region, the Qs frequency dependence can be approximated with the function 33 f 1.22 within a frequency range between 0.5 and 20 Hz. However, the frequency dependence of Qs in the frequency range between 0.5 and 6 Hz is best approximated by Qs （f） = 36 f 0.94 showing relatively weaker frequency dependence as compared to the relation Qs （f） = 6 f^ 2.09 for the frequency range between 6 and 15 Hz. These results could be used to estimate source and site parameters for seismic hazard assessment in the region.
