用宽频带远震体波波形研究中国及邻近地区地壳构造及地质意义初探

The theory and approach of the broadband teleseismic body waveform inversion are expatiated in this paper, and the defining the crust structure's methods are developed. Based on the teleseismic P-wave data, the theoretic image of the P-wave radical component is calculated via the convolution of the teleseismic P-wave vertical component and the transform function, and thereby a P-wavefrom inversion method is built. The applied results show the approach effective, stable and its resolution high. The exact and reliable teleseismic P waveforms recorded by CDSN and IRIS and its geodynamics are utilized to obtain China and its vicinage lithospheric transfer functions, this region ithospheric structure is inverted through the inversion of reliable transfer functions, the new knowledge about the deep structure of China and its vicinage is obtained, and the reliable seismological evidence is provided to reveal the geodynamic evolution processes and set up the continental collisional theory. The major studies are as follows: Two important methods to study crustal and upper mantle structure -- body wave travel-time inversion and waveform modeling are reviewed systematically. Based on ray theory, travel-time inversion is characterized by simplicity, crustal and upper mantle velocity model can be obtained by using 1-D travel-time inversion preliminary, which introduces the reference model for studying focal location, focal mechanism, and fine structure of crustal and upper mantle. The large-scale lateral inhomogeneity of crustal and upper mantle can be obtained by three-dimensional t ravel-time seismic tomography. Based on elastic dynamics, through the fitting between theoretical seismogram and observed seismogram, waveform modeling can interpret the detail waveform and further uncover one-dimensional fine structure and lateral variation of crustal and upper mantle, especially the media characteristics of singular zones of ray. Whatever travel-time inversion and waveform modeling is supposed under certain approximate conditions, with respective advantages and disadvantages, and provide convincing structure information for elucidating physical and chemical features and geodynamic processes of crustal and upper mantle. Because the direct wave, surface wave, and refraction wave have lower resolution in investigating seismic velocity transitional zone, which is inadequate to study seismic discontinuities. On the contrary, both the converse and reflected wave, which sample the discontinuities directly, must be carefully picked up from seismogram to constrain the velocity transitional zones. Not only can the converse wave and reflected wave study the crustal structure, but also investigate the upper mantle discontinuities. There are a number of global and regional seismic discontinuities in the crustal and upper mantle, which plays a significant role in understanding physical and chemical properties and geodynamic processes of crustal and upper mantle. The broadband teleseismic P waveform inversion is studied particularly. The teleseismic P waveforms contain a lot of information related to source time function, near-source structure, propagation effect through the mantle, receiver structure, and instrument response, receiver function is isolated form teleseismic P waveform through the vector rotation of horizontal components into ray direction and the deconvolution of vertical component from the radial and tangential components of ground motion, the resulting time series is dominated by local receiver structure effect, and is hardly irrelevant to source and deep mantle effects. Receiver function is horizontal response, which eliminate multiple P wave reflection and retain direct wave and P-S converted waves, and is sensitive to the vertical variation of S wave velocity. Velocity structure beneath a seismic station has different response to radial and vertical component of an accident teleseismic P wave. To avoid the limits caused by a simplified assumption on the vertical response, the receiver function method is mended. In the frequency domain, the transfer function is showed by the ratio of radical response and vertical response of the media to P wave. In the time domain, the radial synthetic waveform can be obtained by the convolution of the transfer function with the vertical wave. In order to overcome the numerical instability, generalized reflection and transmission coefficient matrix method is applied to calculate the synthetic waveform so that all multi-reflection and phase conversion response can be included. A new inversion method, VFSA-LM method, is used in this study, which successfully combines very fast simulated annealing method (VFSA) with damped least square inversion method (LM). Synthetic waveform inversion test confirms its effectiveness and efficiency. Broadband teleseismic P waveform inversion is applied in lithospheric velocity study of China and its vicinage. According to the data of high quality CDSN and IRIS, we obtained an outline map showing the distribution of Asian continental crustal thickness. Based on these results gained, the features of distribution of the crustal thickness and outline of crustal structure under the Asian continent have been analyzed and studied. Finally, this paper advances the principal characteristics of the Asian continental crust. There exist four vast areas of relatively minor variations in the crustal thickness, namely, northern, eastern southern and central areas of Asian crust. As a byproduct, the earthquake location is discussed, Which is a basic issue in seismology. Because of the strong trade-off between the assumed initial time and focal depth and the nonlinear of the inversion problems, this issue is not settled at all. Aimed at the problem, a new earthquake location method named SAMS method is presented, In which, the objective function is the absolute value of the remnants of travel times together with the arrival times and use the Fast Simulated Annealing method is used to inverse. Applied in the Chi-Chi event relocation of Taiwan occurred on Sep 21, 2000, the results show that the SAMS method not only can reduce the effects of the trade-off between the initial time and focal depth, but can get better stability and resolving power. At the end of the paper, the inverse Q filtering method for compensating attenuation and frequency dispersion used in the seismic section of depth domain is discussed. According to the forward and inverse results of synthesized seismic records, our Q filtrating operator of the depth domain is consistent with the seismic laws in the absorbing media, which not only considers the effect of the media absorbing of the waves, but also fits the deformation laws, namely the frequency dispersion of the body wave. Two post stacked profiles about 60KM, a neritic area of China processed, the result shows that after the forward Q filtering of the depth domain, the wide of the wavelet of the middle and deep layers is compressed, the resolution and signal noise ratio are enhanced, and the primary sharp and energy distribution of the profile are retained.

三峡库区滑坡的时空分布、成因及其预报模式

Landslides are widely distributed along the main stream banks of the Three Gorges Reservoir area. Especially with the acceleration of the human economic activities in the recent 30 years, the occurrence of landslide hazards in the local area trends to be more serious. Because of the special geological, topographic and climatic conditions of the Three Gorges areas, many Paleo-landslides are found along the gentle slope terrain of the population relocation sites. Under the natural condition, the Paleo-landslides usually keep stable. The Paleo-landslides might revive while they are influenced under the strong rainfall, water storage and migration engineering disturbance. Therefore, the prediction and prevention of landslide hazards have become the important problem involving with the safety of migration engineering of the Three Gorges Reservoir area.The past research on the landslides of the Three Gorges area is mainly concentrated on the stability analysis of individual landslide, and importance was little attached to the knowledge on the geological environment background of the formation of regional landslides. So, the relationship between distribution and evolution of landslides and globe dynamic processes was very scarce in the past research. With further study, it becomes difficult to explain the reasons for the magnitude and frequency of major geological hazards in terms of single endogenic or exogenic processes. It is possible to resolve the causes of major landslides in the Three Gorges area through the systematic research of regional tectonics and river evolution history.In present paper, based on the view of coupling of earth's endogenic and exogenic processes, the author researches the temporal and spacial distribution and formation evolution of major landslides(Volume^lOOX 104m3) in the Three Gorges Reservoir area through integration of first-hand sources statistics, .geological evolution history, isotope dating and numerical simulation method etc. And considering the main formation factors of landslides (topography, geology and rainfall condition), the author discusses the occurrence probability and prediction model of rainfall induced landslides.The distribution and magnitude of Paleo-landslides in the Three Gorges area is mainly controlled by lithology, geological structure, bank slope shape and geostress field etc. The major Paleo-landslides are concentrated on the periods 2.7-15.0 X 104aB.R, which conrresponds to the warm and wettest Paleoclimate stages. In the same time, the Three Gorges area experiences with the quickest crust uplift phase since 15.0X 104aB.P. It is indicated that the dynamic factor of polyphase major Paleo-landslides is the coupling processes of neotectonic movement and Quaternary climate changes. Based on the numerical simulation results of the formation evolution of Baota landslide, the quick crust uplift makes the deep river incision and the geostress relief causes the rock body of banks flexible. Under the strong rainfall condition, the pore-water pressure resulted from rain penetration and high flood level can have the shear strength of weak structural plane decrease to a great degree. Therefore, the bank slope is easy to slide at the slope bottom where shear stress concentrates. Finally, it forms the composite draught-traction type landslide of dip stratified rocks.The susceptibility idea for the rainfall induced landslide is put forward in this paper and the degree of susceptibility is graded in terms of the topography and geological conditions of landslides. Base on the integration with geological environment factors and rainfall condition, the author gives a new probabilistic prediction model for rainfall induced landslides. As an example from Chongqing City of the Three Gorges area, selecting the 5 factors of topography, lithology combination, slope shape, rock structure and hydrogeology and 21 kinds of status as prediction variables, the susceptibility zonation is carried out by information methods. The prediction criterion of landslides is established by two factors: the maximum 24 hour rainfall and the antecedent effective precipitation of 15 days. The new prediction model is possible to actualize the real-time regional landslide prediction and improve accuracy of landslide forecast.