特大地震激发地球自由振荡的数值模拟研究

Large earthquakes, such as the Chile earthquake in 1960 and the Sumatra-Andaman earthquake on Dec 26, 2004 in Indonesia, have generated the Earth’s free oscillations. The eigenfrequencies of the Earth’s free oscillations are closely related to the Earth’s internal structures. The conventional methods, which mainly focus on calculating the eigenfrequecies by analytical ways, and the analysis on observations can not easily study the whole processes from earthquake occurrence to the Earth’s free oscillation inspired. Therefore, we try to use numerical method incorporated with large-scale parallel computing to study on the Earth’s free oscillations excited by giant earthquakes. We first give a review of researches and developments of the Earth’s free oscillation, and basical theories under spherical coordinate system. We then give a review of the numerical simulation of seismic wave propagation and basical theories of spectral element method to simulate global seismic wave propagation. As a first step to study the Earth’s free oscillations, we use a finite element method to simulate the propagation of elastic waves and the generation of oscillations of the chime bell of Marquis Yi of Zeng, by striking different parts of the bell, which possesses the oval crosssection. The bronze chime bells of Marquis Yi of Zeng are precious cultural relics of China. The bells have a two-tone acoustic characteristic, i.e., striking different parts of the bell generates different tones. By analysis of the vibration in the bell and the spectrum analysis, we further help the understanding of the mechanism of two-tone acoustic characteristics of the chime bell of Marquis Yi of Zeng. The preliminary calculations have clearly shown that two different modes of oscillation can be generated by striking different parts of the bell, and indicate that finite element numerical simulation of the processes of wave propagation and two-tone generation of the chime bell of Marquis Yi of Zeng is feasible. These analyses provide a new quantitative and visual way to explain the mystery of the two-tone acoustic characteristics. The method suggested by this study can be applied to simulate free oscillations excited by great earthquakes with complex Earth structure. Taking into account of such large-scale structure of the Earth, small-scale low-precision numerical simulation can not simply meet the requirement. The increasing capacity in high-performance parallel computing and progress on fully numerical solutions for seismic wave fields in realistic three-dimensional spherical models, Spectral element method and high-performance parallel computing were incorporated to simulate the seismic wave propagation processes in the Earth’s interior, without the effects of the Earth’s gravitational potential. The numerical simulation shows that, the results of the toroidal modes of our calculation agree well with the theoretical values, although the accuracy of our results is much limited, the calculated peaks are little distorted due to three-dimensional effects. There exist much great differences between our calculated values of spheroidal modes and theoretical values, because we don’t consider the effect the Earth’ gravitation in numerical model, which leads our values are smaller than the theoretical values. When , is much smaller, the effect of the Earth’s gravitation make the periods of spheroidal modes become shorter. However, we now can not consider effects of the Earth’s gravitational potential into the numerical model to simulate the spheroidal oscillations, but those results still demonstrate that, the numerical simulation of the Earth’s free oscillation is very feasible. We make the numerical simulation on processes of the Earth’s free oscillations under spherically symmetric Earth model using different special source mechanisms. The results quantitatively show that Earth’s free oscillations excited by different earthquakes are different, and oscillations at different locations are different for free oscillation excited by the same earthquake. We also explore how the Earth’s medium attenuation will take effects on the Earth’s free oscillations, and take comparisons with the observations. The medium attenuation can make influences on the Earth’s free oscillations, though the effects on lower-frequency fundamental oscillations are weak. At last, taking 2008 Wenchuan earthquake for example, we employ spectral element method incorporated with large-scale parallel computing technology to investigate the characteristics of seismic wave propagation excited by Wenchuan earthquake. We calculate synthetic seismograms with one-point source model and three-point source model respectively. Full 3-D visualization of the numerical results displays the profile of the seismic wave propagation with respect to time. The three-point source, which was proposed by the latest investigations through field observation and reverse estimation, can better demonstrate the spatial and temporal characteristics of the source rupture processes than one-point source. Primary results show that those synthetic signals calculated from three-point source agree well with the observations. This can further reveal that the source rupturing process of Wenchuan earthquake is a multi-rupture process, which is composed by at least three or more stages of rupture processes. In conclusion, the numerical simulation can not only solve some problems concluding the Earth’s ellipticity and anisotropy, which can be easily solved by conventional methods, but also finally solve the problems concluding topography model and lateral heterogeneity. We will try to find a way to fully implement self-gravitation in spectral element method in future, and do our best to continue researching the Earth’s free oscillations using the numerical simulations to see how the Earth’ lateral heterogeneous will affect the Earth’s free oscillations. These will make it possible to bring modal spectral data increasingly to bear on furthering our understanding of the Earth’s three-dimensional structure.

南海张裂大陆边缘数值模拟研究

The South China Sea (SCS) is one of the largest marginal seas in the western Pacific, which is located at the junction of Eurasian plate, Pacific plate and Indian-Australian plate. It was formed by continent breakup and sea-floor spreading in Cenozoic. The complicated interaction among the three major plates made tectonic movement complex and geological phenomena very rich in this area. The SCS is an ideal place to study the formation and evolution of rifted continental margin and sea-floor spreading since it is old enough to have experienced the major stages of the basin evolution but still young enough to have preserved its original nature. As the demand for energy grows day by day in our country, the deep water region of the northern continental margin in the SCS has become a focus of oil and gas exploration because of its huge hydrocarbon potential. Therefore, to study the rifted continental margin of the SCS not only can improve our understanding of the formation and evolution processes of rifted continental margin, but also can provide theoretical support for hydrocarbon exploration in rifted continental margin. This dissertation mainly includes five topics as follows: (1) Various classic lithosphere stretching models are reviewed, and the continuous non-uniform stretching model is modified to make it suitable for the case where the extension of lithopheric mantle exceeds that of the crust. Then simple/pure shear flexural cantilever model is applied to model the basement geometries of SO49-18 profile in the northern continental margin of the SCS. By fitting the basements obtained by using 2DMove software with modeling results, it is found that the reasonable effective elastic thickness is less than 5km in this region. According to this result, it is assumed that there is weak lower crust in the northern continental margin in the SCS. (2) We research on the methods for stretching factor estimation based on various lithosphere stretching models, and apply the method based on multiple finite rifting model to estimate the stretching factors of several wells and profiles in the northern continental margin of the SCS. (3) We improve one-dimension strain rate inversion method with conjugate gradient method, and apply it to invert the strain rate of several wells in the northern continental margin of the SCS. Two-dimension strain rate forward modeling is carried out, and the modeling results show that effective elastic thickness is a key parameter to control basin’s geometry. (4) We simulate divergent upwelling mantle flow model using finite difference method, and apply this newly developed model to examine the formation mechanism of the northwest and central sub-basin in the SCS. (5) We inverse plate thickness and basal temperature of oceanic lithosphere using sea-floor ages and bathymetries of the North Pacific and the North Atlantic based on varied-parameters plate model, in which the heat conductivity, heat capacity and coefficient of thermal expansion depend on temperature or depth. A new empirical formula is put forward based the inversed parameters, which depicts the relation among sea-floor age, bathymetry and heat flow. Then various similar empirical formulae, including the newly developed one, are applied to examine the sea-floor spread issue in the SCS based on the heat flow and bathymetry data of the abyssal sub-basin.

活动构造区区域地壳稳定性研究——以金沙江水电梯级工程规划区为例

The study of regional crustal stability of active tectonic region basically includes analysis of recent activity of Earth's crust, single factor assessment, study of complexity, and comprehensive assessment of crust stability. In this thesis, some work are made as follows: · Based on abundant data from gravity field, aeromagnetic survey, magnetism, magnetotelluric deep sounding, remote sensing and geotectonic as well as earthquakes observed in recent years around this region and adjacent zones, we can get a through understanding about the structural features and activity of the earth's crust in Chuan-Dian region. The results from explosion earthquake and telluric electromagnetic sounding are consistent with the structural features of the crust manifested by the geophysical field. The data of deep geologic structures are important for us to work out a vivid three-dimensional model of the earth's crustal structure of the Jinsha River region. According to a synthesis, the author of this thesis proposes some indicators for dividing the faulted blocks. It can also be inferred that the movement of the Chuan-Dian faulted block, which is the relatively active part of southwestern China, is controlled by the boundary faults, and the intensive activity and deformation are concentrated along the boundaries of the block. · Mainly discussing respectively the mechanism and laws of active faults, earthquakes, and geological hazards activity, and their influences on the stability and security of engineering, also trying to probe into the way to assess the risk of single factor in this section. Especially with the method of fractal geometry, the thesis has discussed how to study the complexity of each factor. These geologic hazards which are distributed at the uppermost part of the crust in this region form a typical mountainous set of the active tectonic areas. The results of survey show that some slopes are liable- to -sliding with a weak layer of low shear strength. Occurrences of landslides are to a great extent related to local geological structures, in particular active faults. This is why numerous landslides have occurred simultaneously around the epicenter of a strong earthquake or the center of a strong rainfall, which are related to active faults. · The analysis of the crustal stability is based on a regional grid division, and a fuzzy comprehensive analysis method is used to determine the grade of the quality in each grid. The evaluation factors and their weights are taken from the results of the hierarchical analysis. The evaluation indexes consist of qualitative and quantitative ones. The qualitative ones can be quantified through the experts weighing system, while the quantitative ones can be obtained from statistical analysis. For quality grades, four levels are used: stable, essentially stable, sub-stable, and unstable. The results of the evaluation on Jinshajiang region demonstrate that the crustal stability become distinctly worse in the areas controlled by active deep faults. Therefore, detailed investigations on the active faulting and geologic hazards, include earthquake activity are especially necessary for those areas adjacent to the deep fault belts. On the bases of the data available and the survey results, we have made a preliminary assessment for the construction conditions and adaptability of every planned site in the middle or lower reaches of Jinsha River. Finally, the thesis prospected the vista of the study of crustal stability.

华南地壳岩石弹性和变形性质的高温高压实验研究

The South continent of China lies to southeast of Eurasia block. It is an active area from the view of crust growth and continent spread and is a transition zone between continental crust and oceanic crust. The compressional wave velocities and anisotropies of typical crustal metamorphic rocks were determined at high temperature (up to 1000 ℃) and high pressure(up to 800MPa). The experimental results show that the velocities generally increase with pressure, and is unaffected by temperature up to around 550 ℃. But the velocities of all experimental samples start to drop above a temperature point. For an example, this greatly reduce the speed of wave propagation in amphibolite and serpentinite above 760 ℃ and above 550 ℃ respectively, which may be due to dehydrate of amphibole and serpentine. P-wave anisotropy coefficients of those rocks range from 2% to 10% at 800MPa and 500 ℃. The anisotropies decrease with increasing pressure at room temperature, but hardly change as function of temperature at constant 800MPa or 600MPa pressure. The average velocity of the six crustal rocks is 6.28km/s under the condition of 800MPa and 550 ℃, which is consistent with the result of deep seismic sounding data. Based on this experimental result, we deduce there may exist a lot of felsic granulites and amphibolites at the depth of 15-25km underground. With increasing temperature and pressure, the deformation behavior of the rocks undergoes from localized brittle fracture, semi-brittle deformation (cataclastic flow or semi-brittle faulting, semi-brittle flow) to homogeneous crystal-plastic flow. This transition is associated with mechanical behavior and micro-mechanism. It is very important to understanding earthquake source mechanics, the strength of the lithosphere and the style of deformation. The experiments were conducted at temperature of 600-1000 ℃, confining pressure of 500MPa, and stain rates of 10~(-4)-10~(-6) S~(-1). For fine-grained natural amphibolite, the results of experiments show that brittle faulting is major failure mode at temperature <600 ℃, but crystal-plastic deformation is dominate at temperature >800 ℃, and there is a transition with increasing temperature from sembrittle faulting to cataclastic flow and sembrittle flow at temperature of 670-750 ℃. For medium-grained natural Felsic granulite, the results of experiments show that brittle faulting is major failure mode at temperature <500 ℃, but crystal-plastic deformation is dominate at temperature >700 ℃, and there is a transition with increasing temperature from semibrittle faulting to cataclastic flow and sembrittle flow at temperature of 500-600 ℃.

边坡动力响应分析及应用研究

As a marginal subject, dynamic responses of slopes is not only an important problem of engineering geology (Geotechnical problem), but also of other subjects such as seismology, geophysics, seismic engineering and engineering seismic and so on. Owning to the gulf between different subjects, it is arduous to study dynamic responses of slopes and the study is far from ripeness. Studying on the dynamic responses of slopes is very important in theories as well as practices. Supported by hundreds of bibliographies, this paper systemically details the development process of this subject, introduces main means to analyze this subject, and then gives brief remarks to each means respectively. Engineering geology qualitative analysis is the base of slopes dynamic responses study. Because of complexity of geological conditions, engineering geology qualitative analysis is very important in slopes stability study, especially to rock slopes with complex engineering geology conditions. Based on research fruits of forerunners, this paper summarizes factors influencing slopes dynamic stability into five aspects as geology background, stratums, rock mass structure, and topography as well as hydrogeology condition. Based on rock mass structure controlling theory, engineering geology model of the slope is grouped into two typical classes, one is model with obvious controlling discontinuities, which includes horizontal bedded slope, bedding slope, anti-dip slope, slide as well as slope with base rock and weathered crust; the other is model without obvious controlling discontinuities, which includes homogeneous soil slope, joint rock mass slope. Study on slope failure mechanism under dynamic force, the paper concludes that there are two effects will appear in slope during strong earthquake, one is earthquake inertia force, the other is ultra pore pressure buildup. The two effects lead to failure of the slope. To different types of slope failure, the intensity of two effects acting on the slope is different too. To plastic flow failure, pore pressure buildup is dominant; to falling rock failure and toppling failure, earthquake inertia force is dominant in general. This paper briefly introduces the principle of Lagrangian element method. Through a lot of numerical simulations with FLAC3D, the paper comprehensively studies dynamic responses of slopes, and finds that: if the slope is low, displacement, velocity and acceleration are linear enlarging with elevation increasing in vertical direction; if the slope is high enough, displacement, velocity and acceleration are not linear with elevation any more, on the other hand, they fluctuate with certain rhythm. At the same time, the rhythm appears in the horizontal direction in the certain area near surface of the slope. The distribution form of isoline of displacement, velocity and acceleration in the section of the slope is remarkably affected by the slope angle. In the certain area near the slope surface, isoline of displacement, velocity and acceleration is parallel to the surface of the slope, in the mean time, the strike direction of the extreraum area is parallel to the surface of the slope too. Beyond this area, the isoline direction and the strike direction of the extremum area turn to horizontal with invariable distance. But the rhythm appearing or not has nothing to with the slope angle. The paper defines the high slope effect and the low slope effect of slopes dynamic responses, discusses the threshold height H^t of the dynamic high slope effect, and finds that AW is proportional to square root of the dynamic elastic moduli El P , at the same time, it is proportional to period Tof the dynamic input. Thus, the discriminant of H^t is achieved. The discriminant can tell us that to a slope, if its height is larger than one fifth of the wavelength, its response regular will be the dynamic high slope effect; on the other hand, its response regular will be the dynamic low slope effect. Based on these, the discriminant of different slopes taking on same response under the same dynamic input is put forward in this paper. At the same time, the paper studies distribution law of the rhythm extremum point of displacement, velocity and acceleration, and finds that there exists relationship of N = int among the slope height H, the number of the rhythm extremum VHlhro) point N and ffthre- Furthermore, the paper points out that if N^l, the response of the slope will be dynamic high slope effect; \fN

中国大陆及邻区地震活动背景、环境触发与地震预测研究

This dissertation that includes most of the P. PH.D research work during 2001~2002 covers the large-scale distribution of continental earthquakes in mainland China, the mechanism and statistic features of grouped strong earthquakes related to the tidal triggering, some results in earthquake prediction with correlativity analysis methods, and the flushes from the two strong continental earthquakes in South Asia in 2001. Mainland China is the only continental sub-plate that is compressed by collision boundaries at its two sides, within which earthquakes are dispersive and distributed as seismic belts with different widths. The control capability of the continental block boundaries on the strong earthquakes and seismic hazards is calculated and analyzed in this dissertation. By mapping the distribution of the 31282 ML:3s2,0 earthquakes, I found that the depth of continental earthquakes depend on the tectonic zonings. The events on the boundaries of relatively integrated blocks are deep and those on the new-developed ruptures are shallow. The average depth of earthquakes in the West of China is about 5km deeper than that in the east. The western and southwestern brim of Tarim Basin generated the deepest earthquakes in mainland China. The statistic results from correlation between the grouped M7 earthquakes and the tidal stress show that the strong events were modulated by tidal stress in active periods. Taking Taiwan area as an example, the dependence of moderate events on the moon phase angles (￡>) is analyzed, which shows that the number of the earthquakes in Taiwan when D is 50° ,50° +90° and 50° +180° is more than 2 times of standard deviation over the average frequency at each degree, corresponding to the 4th, 12th and 19th solar day after the new moon. The probability of earthquake attack to the densely populated Taiwan island on the 4th solar day is about 4 times of that on other solar days. On the practice of earthquake prediction, I calculated and analyzed the temporal correlation of the earthquakes in Xinjinag area, Qinghai-Tibet area, west Yunnan area, North China area and those in their adjacent areas, and predicted at the end of 2000 that it would be a special time interval from 2001 to 2003, within which moderate to strong earthquakes would be more active in the west of China. What happened in 2001 partly validated the prediction. Within 10 months, there were 2 great continental earthquakes in south Asia, i.e., the M7.8 event in India on Jan 26 and M8.1 event in China on Nov. 14, 2001, which are the largest earthquake in the past 50 years both for India and China. No records for two great earthquakes in Asia within so short time interval. We should speculate the following aspects from the two incidences: The influence of the fallacious deployment of seismic stations on the fine location and focal mechanism determination of strong earthquakes must be affronted. It is very important to introduce comparative seismology research to seismic hazard analysis and earthquake prediction research. The improvement or changes in real-time prediction of strong earthquakes with precursors is urged. Methods need to be refreshed to protect environment and historical relics in earthquake-prone areas.

接收函数方法与中国及邻区上地幔间断面的研究

The receiver function method applied in researching the discontinuities in upper mantle was systematically studied in this paper. Using the theoretical receiver functions, the characteristics of P410S and P660S phases were analyzed, and the influencing factors for detection of these phases were discussed. The stability of receiver function was studied, and a new computational method of receiver function, RFSSMS (Receiver Function of Stack and Smooth of Multi seismic-records at a Single station), was put forward. We built initial reference velocity model for the media beneath each of 18 seismic stations respectively; then estimated the buried depths of 410-km and 660-km discontinuities(simply marked as '410' and '660') under the stations by using the arrive time differences of P410S and P660S with P. We developed a new receiver function inversion method -PGARFI (Peeling-Genetic Algorithm of Receiver Function Inversion), to obtain the whole crust and upper mantle velocity structure and the depths of discontinuities beneath a station. The major works and results could be summarized as follows: (1) By analysis of the theoretical receiver functions with different velocity models and different ray parameters, we obtain the knowledge: The amplitudes of P410S and P660S phases are decreasing with the increasing of epicentral distance A , and the arrival time differences of these phases with P are shorter as A is longer. The multiple refracted and/or reflected waves yielded on Moho and the discontinuities in the crust interfere the identification of P410S. If existing LVZ under the lithosphere, some multiple waves caused by LVZ will interfere the identification of P410S. The multiple waves produced by discontinuity lied near 120km depth will mix with P410s phase in some range of epicentral distance; and the multiple waves concerned with the discontinuity lied near 210km depth will interfere the identification of P660S. The epicentral distance for P4i0s identification is limited, the upper limit is 80° . The identification of P660S is not restricted by the epicenter distance obviously. The identification of P410S and P6gos in the theoretical receiver functions is interfered weakly from the seismic wave attenuation caused by the media absorption if the Q value in a reasonable range. (2) The stability of receiver function was studied by using synthetic seismograms with different kind of noise. The results show that on the condition of high signal-noise-ratio of seismic records, the high frequency background noise and the low frequency microseism noise do not influence the calculating result of receiver function. But the media "scattering noise" influence the stability of receiver function. When the scattering effect reach some level, the identification of P4iOs and P66os is difficult in single receiver function which is yielded from only one seismic record. We provided a new method to calculate receiver function, that is, with a group of earthquake records, stacking the R and Z components respectively in the frequency domain, and weighted smooth the stacked Z component, then compute the complex spectrum ratio of R to Z. This method can improve the stability of receiver function and protrude the P4i0s and P66os in the receiver function curves. (3) 263 receiver functions were provided from 1364 three component broadband seismograms recorded at 18 stations in China and adjacent areas for the tele-earthquakes. The observed arrival time differences of P410S and P660S with P were obtained in these receiver functions. The initial velocity model for every station was built according to the prior research results. The buried depths of '410' and '660' under a station were acquired by the way of adjusting the depths of these two discontinuities in the initial velocity model until the theoretical arrival time differences of P410S and P660S with P well conformed to the observed. The results show an obvious lateral heterogeneity of buried depths of ' 410' and (660' . The depth of '410' is shallower beneath BJI, XAN, LZH and ENH, but deeper under QIZ and CHTO, and the average is 403km . The average depth of '660' is 663km, deeper under MDJ and MAJO, but shallower under QIZ and HYB. (4) For inversing the whole crust and upper mantle velocity structure, a new inversion method -PGARFI (Peeling-Genetic Algorithm of Receiver Function Inversion) has- been developed here. The media beneath a station is divided into segments, then the velocity structure is inversed from receiver function from surface to deep successively. Using PGARFI, the multi reflection / refraction phases of shallower discontinuities are isolated from the first order refraction transform phase of deep discontinuity. The genetic algorithm with floating-point coding was used hi the inversion of every segment, and arithmetical crossover and non-uniform mutation technologies were employed in the genetic optimization. 10 independent inversions are completed for every segment, and 50 most excellent velocity models are selected according to the priority of fitness from all models produced in the inversion process. The final velocity structure of every segment is obtained from the weighted average of these 50 models. Before inversion, a wide range of velocity variation with depth and depth range of the main discontinuities are given according to priori knowledge. PGARFI was verified with numerical test and applied in the inversion of the velocity structure beneath HIA station down to 700km depth.

天然地震转换波偏移成像方法研究

Butovskaya, a scholar of Former Soviet Union, first determined the depth of basalt layer in Tashkent Zone by using converted waves on seismogram in 1952. From then on, more and more scholars developed the comprehensive research that imaged the earth interior structures by applying converted waves information. With the digitalization of earthquake observation, The inversion imaging of complete or partial waveform record can efficiently improve inversion quality and widen its usage scope, therefore great progress is made in converted wave imaging. This paper makes a certain study in converted wave imaging on that basis. Transmitted PP waves and converted PS waves are generated when a P-wave propagates through an interface separating two media with large impedance contracts. A PS converted wave is a seismic body wave, which result from the conversion of an incident parent P wave at a boundary within the crust to a refracted S wave. The thickness of a single crustal layer can theoretically be determined by observing, with three-componented seismometer at a single station, the difference in time of the arrival of the parent P wave and the arrival of the PS converted wave. For a multilayered media, PS converted wave arrivals corresponding to each of the layers can theoretically be observed, provided the station is sufficiently from the source of the parent P wave to allow initial penetration of the P wave beneath the deepest layer considered. To avoid the difficulty of picking up transmitted P-wave and converted wave phases, this paper proposed a converted wave migration method by estimating the travel time difference between PS converted wave and PP transmitted wave. To verify its validity, we apply the converted wave PS migration algorithm to synthetic data generated by three forward modeling. The migration results indicate that PS converted wave may be migrated to reconstruct the transmitting interface. This technique is helpful to investigate the deep earth structures by using earthquake data.

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

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.

东营凹陷沙三中亚段湖相浊积岩储层特征及油气成藏规律研究

Turbidity sandstone reservoirs have been an important field of hydrocarbon exploration and development in the basins all over the world, as well as in China. Lithologic pools are composed of turbidity sandstones and other sandstones are frequently found in the Jiyang Depression that is a Mesozoic-Cenozoic non-marine oil-bearing basin. The Dongying Sag lies in the sedimentary center of the basin. The subtle traps with turbidity reservoirs are generally difficult to be predicted and described by using current techniques. The studies on turbidity reservoirs plays thus an important theoretical and theoretical practical role in exploration and development in the Jiyang Depression. The attention is, in this thesis, focused on the petrologic properties and oil accumulating behaviors in lake turbidity sedimentary systems in the middle part of the third section of Shahejie Formation in the Dongying Sag, especially in Dongxin area, which lies on the central uplift of the Sag. The paper has disclosed the origin types of turbidity sandstones, distribution pattern and controlling factors of turbidity sandstones, and set up hydrocarbon accumulation patterns of the middle part of the third section of Shahejie Formation in Dongxin, based on nonmarine high resolution sequence stratigraphy, event sedimentology and new theories of hydrocarbon forming. By studying prediction method and technology of turbidity sandstone reservoirs, using precise geological model developing, new techniques of high resolution seismic inversion constrained by logging, the paper has forecast low permeability turbidity sandstone reservoirs and pointed out advantage exploration aims to progressive exploration and development. The paper has obtained mainly many productions and acknowledges as follows: 1.Turbidity sandstone reservoirs of the third section of Shahejie Formationin Dongying Sag are formed in such specifical geological background as rift and extension of basin. The inherited Dongying delta and transgression make up many turbidity distribution areas by overlaying and joining together. The hydrocarbon migrates from depression area to adjacent turbidity sandstone continuously. Accumulation area which is sufficient in oil is formed. 2.The paper has confirmed distinguishable sign of sequence boundary , established stratigraphic framework of Dongying Sag and realized isotime stratigraphic correlation. Es3 of Dongying delta is divided into eleven stages. Among them, the second period of the lower section in Es3, the sixth period of the middle section in Es3, the third period of the upper section in Es3 correspond to eleven sedimentary isotime surface in seismic profile, namely Es3 is classified into eleven Formations. 3.According to such the features of turbidity sandstone as deep in burial, small in area, strong in subtle property, overlaying and joining together and occurring in groups, management through fault and space variations of restriction quantum are realized and the forecast precision of turbidity sandstone by using precise geological model developing, new techniques of high resolution seismic inversion constrained by logging, based on the analysis of all kinds of interwell seismic inversion techniques. 4.According to the features of low permeable turbidity sandstone reservoirs, new method of log interpretation model is put forward. At the same time, distinguish technology of familiar low resistivity oil layer in the turbidity sandstone reservoirs is studied based on petrophysical laboratory work and "four properties" interrelationship between lithological physical Jogging and bearing hydrocarbon properties. Log interpretation model and reservoir index interpretation model of low resistivity oil layer are set up. So the log interpretation precision is improved. 5.The evolution law and its difference of the turbidity sandstone are embodies as follows: the source of sediments come from the south and east of the study area in the middle period of Es3. East source of sediments is pushed from west to east. However, the south source supply of sediments in the early and middle period of Es3 is in full, especially in Es3. subsequently, the supply is decreased gradually. Turbidity fan moves back toward the south and the size of fan is minished accordingly. The characteristic of turbidity sandstone in Dongying Sag is different in different structural positions. Dongxin in the middle-east of the central lift and Niuzhuang Sag He in Dongying delta front and prodelta deep lake subfacies. Although the turbidity sandstone of the two areas root in the Dongying delta sedimentary system, the sand body has different remarkably characteristic. 6.The sedimentary model of the turbiditys in study area have three types as follows: (1) collapse turbidity fan in respect of delta; (2) fault trench turbidity fan; (3) other types of microturbidity sandstone. Middle fan and outer fan, can be found mainly in sublacustrine fan. Middle fan includes braided channel microfacies, central microfacies and braided interchannel microfacies, which is main prospecting oil-bearing subfacies. The middle section of the third section of Shahejie Formation in study area (for example the central lift) can be divided into middle-lower and upper part. The middle-lower part is characteristic of turbidity fan. The upper part is sedimented mainly by delta-collapse fan. 7.The turbidity reservoirs of the middle part of the third section of Shahejie Formation in study area characterize by low maturity both in component and texture, strong in diagenesis and low in permeability. The reservoir can be classified into four types. Type III is the body of reservoir and comprises two types of H a and HI b. M a belongs to middle porosity - low permeability reservoir and distributes in the central lift. Hlb belongs to low porosity - low permeability and distributes in Haojia region. 8.A11 single sand body of lens turbidity reservoir of the middle part of the third section of Shahejie Formation in study area are surrounded by thick dark source rocks. The oil-water system is complex and behaves that every sandstone is single seal unit. The water body is 1/3-1-5 of the sand body. The edge water is not active. The gas exists in the top of reservoir in the form of mixed gas. For far-range turbidity fan with big scale channel, the area and volume of sand body is large and the gap is big in oil packing degree. There are lots of edge water and bottom water, and the latter increases rapidly during the course of development. 9.By exerting the modern hydrocarbon forming theories, the third section of Shahejie Formation in study area belongs to abnormally pressured fluid compartment. The lithological reservoir of the third section of Shahejie Formation is formed in the compartment. The reservoir-formed dynamic system belongs to lower self-source enclosed type. The result and the practice indicate that the form and accumulation of lithological oil reservoirs are controlled by the temperature and pressure of stratum, microfacies, thickness of sand body, fault and reservoir heterogeneity. 10. Based on studies above, the emphases focus on in south and north part of Dongying structure, west Dongxin region and south part Xinzhen structure in the application of production. The practice proves that the turbidity sandstone reservoirs in Ying 11 block and the fault-lithological reservoirs in Xin 133 block have been obtained significant breakthrough. The next target is still sandstone groups of the third section of Shahejie Formation in the bordering areas of Dongxin region for instance Xin 149 area, He 89 area, Ying 8 area etc.

西北太平洋某海域40ka来的地磁场相对强度研究

Paleointensity changes of geomagnetic field help us to understand the evolutionprocess of earth completely and provide further constraints for earth interior process and geodynamo model. Marine sediments are good carriers for relative paleointensity of geomagnetic field. But in most cases, deep sea sediments that conform with magnetic "uniformity" usually have low sedimentation rate about l-2cm/ka and lie under the Carbonate Compensate Depth with little carbonate content. Therefore, the number of relative paleointensity records with detailed oxygen stratigraphy is still rare. This thesis focus on four cores from east of Ryukyu Trench which have foraminiferal content of 5-30% and sedimentation rate of lOcm/ka and wish to get centennial -millennial changes of relative paleointensity.The sediments from east of Ryukyu Trench conform with magnetic "uniformity" and remanences of four cores all show single component with stable direction and faithfully record the magnetic field. The NRM301T1T/ARM and NRMsomT/ x are still affected by grain size and concentration changes although the sediments are "uniform" , indicating the uniformity might not enough for relative paleointensity. After renomalized by grain size parameter MDF, the intensities remove the effect of grain size changes to different degrees and show coherency in 1-1 Oka scale with results from ODP983/984. The characteristics of paleointensity of geomagnetic field arefrom 32-24kaBP, paleointensity of geomagnetic field is low;24-12kaBP, paleointensity of geomagnetic field is high and shows two peaks boundary with 19kaBP.3) 12-5.3ka, paleointensity is low. Then increase from 5,3kaBP until a small trough at 2.7kaBP, then increase till now.

汶川地震后青少年创伤后应激障碍症状及其相关因素研究

The present paper studies focus on the symptoms of the post-traumatic stress disorder in adolescents post-disaster. 482 students from 6 secondary schools and 17 primary schools in the extremely severe disaster areas in Mianzhu, Sichuan province and 785 students from 3 primary schools and 9 secondary schools in the severe disaster areas in Baoji, Shaanxi respectively were surveyed on the symptoms of the post-traumatic stress disorder and the extent of disaster exposure after the Wenchuan earthquake. Self-compiled background information questionnaire and CRIES were used for the investigation. In this study, we contrast the extent of disaster exposure in the two areas in order to explore the related factors about the post-traumatic stress disorder in adolescent post-disaster. The main results of this paper can be summarized as follows： 1. There are significant positive corrections between the post-traumatic stress disorder and the extent of disaster exposure(get trapped in the earthquake、relatives and friends have been injured in the earthquake、look at relatives and friends dying in the earthquake).The more exposed in the disaster, the more serious symptom of the post-traumatic stress disorder. The trauma exposure indicators (get trapped in the earthquake, relatives and friends have been injured in the earthquake、look at relatives and friends dying in the earthquake)were all significant predictors for PTSD severity. 2. There are significant sex difference in the extent（F=8.750, p <0.05） and the incidence rate of PTSD（χ =20.735, df=5,p =0.001）, the extent and the incidence rate of girls in Mianzhu is significantly higher than that of boys. 3. The age is also an influence factor of PTSD. The extent （F=7.246, p <0.001）and the incidence rate （χ =20.735, df=5,p =0.001）of PTSD get higher as adolescent in Mianzhu get older. 4. As the extremely severe disaster areas, the extent of disaster exposure of Mianzhu areas significantly higher than that of the severe disaster areas Baoji. However, there are not difference in the extent of PTSD between two areas（t=0.181，df=1265,p=0.857）, there are only significant difference in the incidence rate of PTSD between two areas（χ =8.766，df=1,p=0.003）, the incidence rate of PTSD in Mianzhu areas significantly higher than that of Baoji areas.