For the purpose of investigating the influence of metastable olivine(MO) phase transformations on both deep seismicity and stagnation of slabs,we constructed a 2-dimensional finite element thermal model for a 120 Ma...For the purpose of investigating the influence of metastable olivine(MO) phase transformations on both deep seismicity and stagnation of slabs,we constructed a 2-dimensional finite element thermal model for a 120 Ma-old 50°dipping oceanic lithosphere descending at 10 cm/yr with velocity boundary layers,which would mitigate the interference of constant velocity field for the slab. The resulting temperatures show that most of intermediate and deep earthquakes occurring within the Tonga slab are occurring inside the 800℃and 1200℃isotherm,respectively.The elevation of olivine transformation near~410 km and respective persistence of metastable olivine and spinel within the transition zone and beneath 660 km would thus result in bimodal positive,zonal,negative density anomalies,respectively.These results together with the resulting pressure anomalies may reflect the stress pattern of the Tonga slab:(i) slab pull force exerts above a depth of~230 km;(ii) MO existence changes the buoyancy force within the transition zone and facilitates slab stagnation at a depth of 660 km;(iii) as the subducting materials accumulated over 660 km,deepest earthquakes occur due to MO transformation;(iv) a flattened‘slab’ may penetrate into the lower mantle due to the density increment of Sp transformation.展开更多
A rapid and accurate method for two-point ray tracing in horizontally layered velocity model is presented in this paper. Numerical experiments show that this method provides stable and rapid convergence with high accu...A rapid and accurate method for two-point ray tracing in horizontally layered velocity model is presented in this paper. Numerical experiments show that this method provides stable and rapid convergence with high accuracies, regardless of various 1-D velocity structures, takeoff angles and epicentral distances. This two-point ray tracing method is compared with the pseudobending technique and the method advanced by Kim and Baag (2002). It turns out that the method in this paper is much more efficient and accurate than the pseudobending technique, but is only applicable to 1-D velocity model. Kim's method is equivalent to ours for cases without large takeoff angles, but it fails to work when the takeoff angle is close to 90°. On the other hand, the method presented in this paper is applicable to cases with any takeoff angles with rapid and accurate convergence. Therefore, this method is a good choice for two-point ray tracing problems in horizontally layered velocity model and is efficient enough to be applied to a wide range of seismic problems.展开更多
The velocity ratio of a free shear layer has an important influence on the spatial development of the large scale coherent structures in the layer. In this study, numerical simulations are performed to get an insight ...The velocity ratio of a free shear layer has an important influence on the spatial development of the large scale coherent structures in the layer. In this study, numerical simulations are performed to get an insight into this problem. The obtained numerical results agree quite well with those of a linear inviscid stability theory and the available experimental data.展开更多
Based on S wave records of deep teleseisms on Digital Seismic Network of Shanxi Province, shear wave velocity structures beneath 6 stations were obtained by means of S wave waveform fitting. The result shows that the ...Based on S wave records of deep teleseisms on Digital Seismic Network of Shanxi Province, shear wave velocity structures beneath 6 stations were obtained by means of S wave waveform fitting. The result shows that the crust is thick in the studied region, reaching 40 km in thickness under 4 stations. The crust all alternatives high velocity layer with low velocity one. There appear varied velocity structures for different stations, and the stations around the same tectonic region exhibit similar structure characteristics. Combined with dominant depth distribution of many small-moderate earthquakes, the correlation between seismogenic layers and crustal structures of high and low velocity layers has been discussed.展开更多
The conventional acoustic logging interpretation method, which is based on vertical wells that penetrate isotropic formations, is not suitable for horizontal and deviated wells penetrating anisotropic formations. This...The conventional acoustic logging interpretation method, which is based on vertical wells that penetrate isotropic formations, is not suitable for horizontal and deviated wells penetrating anisotropic formations. This unsuitability is because during horizontal and deviated well drilling, cuttings will splash on the well wall or fall into the borehole bottom and form a thin bed of cuttings. In addition, the high velocity layers at different depths and intrinsic anisotropy may affect acoustic logging measurements. In this study, we examine how these factors affect the acoustic wave slowness measured in horizontal and deviated wells that are surrounded by an anisotropic medium using numerical simulation. We use the staggered-grid finite difference method in time domain (FDTD) combined with hybrid-PML. First, we acquire the acoustic slowness using a simulated array logging system, and then, we analyze how various factors affect acoustic slowness measurements and the differences between the effects of these factors. The factors considered are high-velocity layers, thin beds of cuttings, dipping angle, formation thickness, and anisotropy. The simulation results show that these factors affect acoustic wave slowness measurements differently. We observe that when the wavelength is much smaller than the distance between the borehole wall and high velocity layer, the true slowness of the formation could be acquired. When the wavelengths are of the same order (i.e., in the near-field scenarios), the geometrical acoustics theory is no longer applicable. Furthermore, when a thin bed of cuttings exists at the bottom of the borehole, Fermat's principle is still applicable, and true slowness can be acquired. In anisotropic formations, the measured slowness changes with increments in the dipping angle. Finally, for a measurement system with specific spacing, the slowness of a thin target layer can be acquired when the distance covered by the logging tool is sufficiently long. Based on systematical simulations with different dipping angles and anisotropy in homogenous TI media, slowness estimation charts are established to quantitatively determine the slowness at any dipping angle and for any value of the anisotropic ratio. Synthetic examples with different acoustic logging tools and different elastic parameters demonstrate that the acoustic slowness estimation method can be conveniently applied to horizontal and deviated wells in TI formations with high accuracy.展开更多
The Benzilan-Tangke deep seismic sounding profile in the western Sichuan region passes through the Song-pan-Garze orogenic belt with trend of NNE. Based on the travel times and the related amplitudes of phases in the ...The Benzilan-Tangke deep seismic sounding profile in the western Sichuan region passes through the Song-pan-Garze orogenic belt with trend of NNE. Based on the travel times and the related amplitudes of phases in the record sections, the 2-D P-wave crustal structure was ascertained in this paper. The velocity structure has quite strong lateral variation along the profile. The crust is divided into 5 layers, where the first, second and third layer belong to the upper crust, the forth and fifth layer belong to the lower crust. The low velocity anomaly zone gener-ally exists in the central part of the upper crust on the profile, and it integrates into the overlying low velocity basement in the area to the north of Ma'erkang. The crustal structure in the section can be divided into 4 parts: in the south of Garze-Litang fault, between Garze-Litang fault and Xianshuihe fault, between Xianshuihe fault and Longriba fault and in the north of Longriba fault, which are basically coincided with the regional tectonics division. The crustal thickness decreases from southwest to northeast along the profile, that is, from 62 km in the region of the Jinshajiang River to 52 km in the region of the Yellow River. The Moho discontinuity does not obviously change across the Xianshuihe fault based on the PmP phase analysis. The crustal average velocity along the profile is lower, about 6.30 km/s. The Benzilan-Tangke profile reveals that the crust in the study area is orogenic. The Xianshuihe fault belt is located in the central part of the profile, and the velocity is positive anomaly on the upper crust, and negative anomaly on the lower crust and upper mantle. It is considered as a deep tectonic setting in favor of strong earthquake's accumulation and occurrence.展开更多
The 1°×1° distribution map of crustmantle structural ratio R for the lithosphere along the Longitudinal Seismic Belt of China has been compiled using computer based on the results of geophysical prospec...The 1°×1° distribution map of crustmantle structural ratio R for the lithosphere along the Longitudinal Seismic Belt of China has been compiled using computer based on the results of geophysical prospecting by previous researchers, and the latest results by the present authors. Based on this map, an insight into the structural features of the crustmantle assemblage along the Longitudinal Seismic Belt has been gained, while their relation to seismic activity and the distributions of geothermal flux and intracrustal high conductivitylow velocity layers, as well as their tectonic effect to seismicity have been discussed.展开更多
In seismic exploration for coal, seismic waves are very difficult to transmit downward because of high velocity protective layers, making the reflection information very hard to receive above ground. Based on the Snel...In seismic exploration for coal, seismic waves are very difficult to transmit downward because of high velocity protective layers, making the reflection information very hard to receive above ground. Based on the Snell law and the Zoeppritz equation, we studied the relationship between the incidence angle and reflection seismic wave energy using a forward model of level media. The result shows that the seismic wave energy has a sudden increase at the critical angle. Based on the energy propagation rule, using big offset to receive the seismic wave energy under a protective layer can effectively reduce its protection effect.展开更多
Through deriving expressions relating the dip-angle ( m) of the lower boundary of a layer to the acoustic velocity (υm) of the layer and other pre-determinable parameters, υm and m can be taken as simultaneously ite...Through deriving expressions relating the dip-angle ( m) of the lower boundary of a layer to the acoustic velocity (υm) of the layer and other pre-determinable parameters, υm and m can be taken as simultaneously iterative variables while solving Shah 's equations . Consequently the previous method of computing υm and m presented by ZHANG S . is improved [1] , and the accuracy of solutions increased greatly.展开更多
It has been known that the error of measuring acoustic veloicities of thin sediment layers by the well-known T2-X2 approach is usually untolerable, and that this approach is unavailable in the case where sea-bed is ha...It has been known that the error of measuring acoustic veloicities of thin sediment layers by the well-known T2-X2 approach is usually untolerable, and that this approach is unavailable in the case where sea-bed is hard because no echo from any subsurface below sea-bottom can be received. Therefore applying the ray-parameter method to thin layers and the refraction method to hard layers need to be considered in an acoustic velocity measurement system composed of a sound source and a towed hydrophone streamer. Some problems of practical importance about the applications of the two methods, such as echo-data processing procedures and error estimations in measuring acoustic veloicities, are discussed, and the effectiveness of theoretical analyses has been verified through computer simulations.展开更多
基金supported by the National Natural Science Foundation of China(Nos.40574047 and 40628004)
文摘For the purpose of investigating the influence of metastable olivine(MO) phase transformations on both deep seismicity and stagnation of slabs,we constructed a 2-dimensional finite element thermal model for a 120 Ma-old 50°dipping oceanic lithosphere descending at 10 cm/yr with velocity boundary layers,which would mitigate the interference of constant velocity field for the slab. The resulting temperatures show that most of intermediate and deep earthquakes occurring within the Tonga slab are occurring inside the 800℃and 1200℃isotherm,respectively.The elevation of olivine transformation near~410 km and respective persistence of metastable olivine and spinel within the transition zone and beneath 660 km would thus result in bimodal positive,zonal,negative density anomalies,respectively.These results together with the resulting pressure anomalies may reflect the stress pattern of the Tonga slab:(i) slab pull force exerts above a depth of~230 km;(ii) MO existence changes the buoyancy force within the transition zone and facilitates slab stagnation at a depth of 660 km;(iii) as the subducting materials accumulated over 660 km,deepest earthquakes occur due to MO transformation;(iv) a flattened‘slab’ may penetrate into the lower mantle due to the density increment of Sp transformation.
基金National Natural Science Foundation of China (40074008 and 40134010).
文摘A rapid and accurate method for two-point ray tracing in horizontally layered velocity model is presented in this paper. Numerical experiments show that this method provides stable and rapid convergence with high accuracies, regardless of various 1-D velocity structures, takeoff angles and epicentral distances. This two-point ray tracing method is compared with the pseudobending technique and the method advanced by Kim and Baag (2002). It turns out that the method in this paper is much more efficient and accurate than the pseudobending technique, but is only applicable to 1-D velocity model. Kim's method is equivalent to ours for cases without large takeoff angles, but it fails to work when the takeoff angle is close to 90°. On the other hand, the method presented in this paper is applicable to cases with any takeoff angles with rapid and accurate convergence. Therefore, this method is a good choice for two-point ray tracing problems in horizontally layered velocity model and is efficient enough to be applied to a wide range of seismic problems.
基金Project supported by the National Natural Science Foundation of china
文摘The velocity ratio of a free shear layer has an important influence on the spatial development of the large scale coherent structures in the layer. In this study, numerical simulations are performed to get an insight into this problem. The obtained numerical results agree quite well with those of a linear inviscid stability theory and the available experimental data.
基金State Key Basic Development and Programming Project Mechanism and Prediction of Continental Strong Earthquakes (G1998040705).
文摘Based on S wave records of deep teleseisms on Digital Seismic Network of Shanxi Province, shear wave velocity structures beneath 6 stations were obtained by means of S wave waveform fitting. The result shows that the crust is thick in the studied region, reaching 40 km in thickness under 4 stations. The crust all alternatives high velocity layer with low velocity one. There appear varied velocity structures for different stations, and the stations around the same tectonic region exhibit similar structure characteristics. Combined with dominant depth distribution of many small-moderate earthquakes, the correlation between seismogenic layers and crustal structures of high and low velocity layers has been discussed.
基金supported by National Natural Science Foundation of China(No.41204094)Science Foundation of China University of Petroleum,Beijing(No.2462015YQ0506)
文摘The conventional acoustic logging interpretation method, which is based on vertical wells that penetrate isotropic formations, is not suitable for horizontal and deviated wells penetrating anisotropic formations. This unsuitability is because during horizontal and deviated well drilling, cuttings will splash on the well wall or fall into the borehole bottom and form a thin bed of cuttings. In addition, the high velocity layers at different depths and intrinsic anisotropy may affect acoustic logging measurements. In this study, we examine how these factors affect the acoustic wave slowness measured in horizontal and deviated wells that are surrounded by an anisotropic medium using numerical simulation. We use the staggered-grid finite difference method in time domain (FDTD) combined with hybrid-PML. First, we acquire the acoustic slowness using a simulated array logging system, and then, we analyze how various factors affect acoustic slowness measurements and the differences between the effects of these factors. The factors considered are high-velocity layers, thin beds of cuttings, dipping angle, formation thickness, and anisotropy. The simulation results show that these factors affect acoustic wave slowness measurements differently. We observe that when the wavelength is much smaller than the distance between the borehole wall and high velocity layer, the true slowness of the formation could be acquired. When the wavelengths are of the same order (i.e., in the near-field scenarios), the geometrical acoustics theory is no longer applicable. Furthermore, when a thin bed of cuttings exists at the bottom of the borehole, Fermat's principle is still applicable, and true slowness can be acquired. In anisotropic formations, the measured slowness changes with increments in the dipping angle. Finally, for a measurement system with specific spacing, the slowness of a thin target layer can be acquired when the distance covered by the logging tool is sufficiently long. Based on systematical simulations with different dipping angles and anisotropy in homogenous TI media, slowness estimation charts are established to quantitatively determine the slowness at any dipping angle and for any value of the anisotropic ratio. Synthetic examples with different acoustic logging tools and different elastic parameters demonstrate that the acoustic slowness estimation method can be conveniently applied to horizontal and deviated wells in TI formations with high accuracy.
基金National Key Research Development Project (No. G1998 040700/sub-project 95-13-02-03).
文摘The Benzilan-Tangke deep seismic sounding profile in the western Sichuan region passes through the Song-pan-Garze orogenic belt with trend of NNE. Based on the travel times and the related amplitudes of phases in the record sections, the 2-D P-wave crustal structure was ascertained in this paper. The velocity structure has quite strong lateral variation along the profile. The crust is divided into 5 layers, where the first, second and third layer belong to the upper crust, the forth and fifth layer belong to the lower crust. The low velocity anomaly zone gener-ally exists in the central part of the upper crust on the profile, and it integrates into the overlying low velocity basement in the area to the north of Ma'erkang. The crustal structure in the section can be divided into 4 parts: in the south of Garze-Litang fault, between Garze-Litang fault and Xianshuihe fault, between Xianshuihe fault and Longriba fault and in the north of Longriba fault, which are basically coincided with the regional tectonics division. The crustal thickness decreases from southwest to northeast along the profile, that is, from 62 km in the region of the Jinshajiang River to 52 km in the region of the Yellow River. The Moho discontinuity does not obviously change across the Xianshuihe fault based on the PmP phase analysis. The crustal average velocity along the profile is lower, about 6.30 km/s. The Benzilan-Tangke profile reveals that the crust in the study area is orogenic. The Xianshuihe fault belt is located in the central part of the profile, and the velocity is positive anomaly on the upper crust, and negative anomaly on the lower crust and upper mantle. It is considered as a deep tectonic setting in favor of strong earthquake's accumulation and occurrence.
文摘The 1°×1° distribution map of crustmantle structural ratio R for the lithosphere along the Longitudinal Seismic Belt of China has been compiled using computer based on the results of geophysical prospecting by previous researchers, and the latest results by the present authors. Based on this map, an insight into the structural features of the crustmantle assemblage along the Longitudinal Seismic Belt has been gained, while their relation to seismic activity and the distributions of geothermal flux and intracrustal high conductivitylow velocity layers, as well as their tectonic effect to seismicity have been discussed.
基金Projects 40574058 supported by the National Natural Science Foundation of China2005cb221500 by the National Basic Research Program of China 03(2007) by the Scientific and Technological Project about Geology and Mineral Resources of Henan Land Resources Department
文摘In seismic exploration for coal, seismic waves are very difficult to transmit downward because of high velocity protective layers, making the reflection information very hard to receive above ground. Based on the Snell law and the Zoeppritz equation, we studied the relationship between the incidence angle and reflection seismic wave energy using a forward model of level media. The result shows that the seismic wave energy has a sudden increase at the critical angle. Based on the energy propagation rule, using big offset to receive the seismic wave energy under a protective layer can effectively reduce its protection effect.
文摘Through deriving expressions relating the dip-angle ( m) of the lower boundary of a layer to the acoustic velocity (υm) of the layer and other pre-determinable parameters, υm and m can be taken as simultaneously iterative variables while solving Shah 's equations . Consequently the previous method of computing υm and m presented by ZHANG S . is improved [1] , and the accuracy of solutions increased greatly.
文摘It has been known that the error of measuring acoustic veloicities of thin sediment layers by the well-known T2-X2 approach is usually untolerable, and that this approach is unavailable in the case where sea-bed is hard because no echo from any subsurface below sea-bottom can be received. Therefore applying the ray-parameter method to thin layers and the refraction method to hard layers need to be considered in an acoustic velocity measurement system composed of a sound source and a towed hydrophone streamer. Some problems of practical importance about the applications of the two methods, such as echo-data processing procedures and error estimations in measuring acoustic veloicities, are discussed, and the effectiveness of theoretical analyses has been verified through computer simulations.