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

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.

华南震旦纪-寒武纪两期成磷事件及其地球动力学意义

Phosphorus is an important biological and ecological element that to a certain degree constrains ecological environment and nutrient (including carbon) cycling. Marine sedimentary phosphorites are the principal phosphorus supply of the mankind. In the eastern to southern margins of the Yangtze Craton, South China, there are two phosphogenetic events at the Doushantuo stage of the Late Sinian and the Meishucun stage of the Early Cambrian respectively, corresponding two explosion events of life across the Precambrian\Cambrian boundary. Phosphorus ores from the Sinian and Cambrian phosphate in South China can be classified roughly into two categories, namely, grained and non-grained phosphorites. Grained phosphorites, hosted in dolostone type of phosphogenetic sequences and with larger industrial values, occur mainly in margins of the Upper Yangtze Platform, formed in shallow-water environments with high hydraulic energy and influenced by frequent sea-level change. Non-grained phosphorites, hosted principally in black-shale type of phosphogenetic sequences and with smaller industrial values, are distributed mainly in the Jiangnan region where deeper-water sub-basins with low hydraulic energy were prevailing at the time of phosphogenesis. Secular change ofδ~(13)C, δ~(18) O, ~(86)Sr/~(87)Sr values of carbonates from Sinian and Cambrian sequences were determined. A negative abnormal ofδ~(13)C, δ~(18)O values and positive abnormal of 86Sr/87Sr values from the fossiliferous section of the Lowest Cambrian Meishucun Formation implies life depopulation and following explosion of life across the PrecambriamCambrian boundary. Based on a lot of observations, this paper put forward a six-stage genetic model describing the whole formational process of industrial phosphorites: 1) Phosphorus was transported from continental weathering products and stored in the ocean; 2) dissolved phosphates in the seawater were enriched in specific deep seawater layer; 3) coastal upwelling currents took this phosphorus-rich seawater to a specific coastal area where phosphorus was captured by oceanic microbes; 4) clastic sediments in this upwelling area were enriched in phosphorus because of abundant phosphorus-rich organic matters and because of phosphorus absorption on grain surfaces; 5) during early diagenesis, the phosphorus enriched in the clastic sediments was released into interstitial water by decomposition and desorption, and then transported to the oxidation-reduction interface where authigenic phosphates were deposited and enriched; 6) such authigenic phosphate-rich layers were scoured, broken up, and winnowed in shallow-water environments resulting in phosphate enrichment. The Sinian-Cambrian phosphorites in South China are in many aspects comparable with coastal-upwelling phosphorites of younger geological ages, especially with phosphorites from modern coastal upwelling areas. That implies the similarities between the Sinian-Cambrian ocean and the modern ocean. Although Sinian-Cambrian oceanic life was much simpler than modern one, but similar oceanic planktons prevail, because oceanic planktons (particularly phytoplanktons) are crucial for phosphate enrichment related to coastal upwelling. It implies also a similar seawater-layering pattern between the Sinian-Cambrian ocean and the modern ocean. The two global phosphate-forming events and corresponding life-explosion events at the Sinian and Cambrian time probably resulted from dissolved-phosphate accumulation in seawater over a critical concentration during the Earth's evolution. Such an oceanic system with seawater phosphorus supersaturation is evidently unstable, and trends to return to normal state through phosphate deposition. Accordingly, this paper put forward a new conception of "normal state <=> phosphorus-supersaturation state" cycling of oceanic system. Such "normal state <=> phosphorus-supersaturation state" cycling was not only important for the three well-known global phosphate-forming events, also related to the critical moments of life evolution on the Earth. It might be of special significance. The favorable paleo-oceanic orientation in regard to coastal-upwelling phosphorite formation suggests a different orientation of the Yangtze Craton between the Sinian time and the present time (with a 135° clockwise difference), and a 25° anti-clockwise rotation of the Yangtze Craton from late Sinian to early Cambrian. During the Sinian-Cambrian time, the Yangtze Craton might be separated from the Cathaysia Block, but might be still associated with the North China Craton.

微功耗海底地震仪的研制

Our motherland has large area of maritime space. Searching and developing ocean becomes more and more important. So Ocean Bottom Seismometer (OBS) as an absolutely necessary equipment can be used in many oceanic fields. OBS not only is an important instrument for discovering structure of lithosphere of ocean bottom, but also plays a main role of oceanic geophysical exploration. The paper introduces my relational work. The MCI micro-power broad frequency seismometer was developed independently. Its power dissipation is less than 300mW. It has some merits including miniature volumeN light mass and cheap price. It is an ideal device not only for the collection high-resolution natural seismic data, but also for the fields of seismic sounding and engineering seismology. Many new high technique were applied to develop this instrument including over-sampling A/D converter, high performance 32bit Micro Process Unit and Flash memory with smart-media interface. Base on the achievement, I have accomplished the showpiece of OBS, which is applied to the deepwater oil and gas geophysical exploration. Because of micro-power dissipation, the seismograph and the sonar releaser can be integrated into a sphere cabin. By this means, the instrument's frequency of resonance and frequency of couple are improved obviously. The data acquisition system of OBS is improved from MCI seismometer. The capacity of flash memory is enlarge from 1G bytes to 8G bytes. The advance MPU in data acquisition system is used to integrate other function modules such as sonar, GPS, compass and digital transmitter.

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

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.

青藏高原东部电性结构特征及地球动力学意义

As the largest and highest plateau on the Earth, the Tibetan Plateau has been a key location for understanding the processes of mountain building and plateau formation during India-Asia continent-continent collision. As the front-end of the collision, the geological structure of eastern Tibetan Plateau is very complex. It is ideal as a natural laboratory for investigating the formation and evolution of the Tibetan Plateau. Institute of Geophysics, Chinese Academy of Sciences (CAS) carried out MT survey from XiaZayii to Qingshuihe in the east part of the plateau in 1998. After error analysis and distortion analysis, the Non-linear Conjugate Gradient inversion(NLCG), Rapid Relaxation Inversin (RRI) and 2D OCCAM Inversion algorithms were used to invert the data. The three models obtained from 3 algorithms provided similar electrical structure and the NLCG model fit the observed data better than the other two models. According to the analysis of skin depth, the exploration depth of MT in Tibet is much more shallow than in stable continent. For example, the Schmucker depth at period 100s is less than 50km in Tibet, but more than 100km in Canadian Shield. There is a high conductivity layer at the depth of several kilometers beneath middle Qiangtang terrane, and almost 30 kilometers beneath northern Qiangtang terrane. The sensitivity analysis of the data predicates that the depth and resistivity of the crustal high conductivity layer are reliable. The MT results provide a high conductivity layer at 20~40km depth, where the seismic data show a low velocity zone. The experiments show that the rock will dehydrate and partially melt in the relative temperature and pressure. Fluids originated from dehydration and partial melting will seriously change rheological characteristics of rock. Therefore, This layer with low velocity and high conductivity layer in the crust is a weak layer. There is a low velocity path at the depth of 90-110 km beneath southeastern Tibetan Plateau and adjacent areas from seismology results. The analysis on the temperature and rheological property of the lithosphere show that the low velocity path is also weak. GPS measurements and the numerical simulation of the crust-mantle deformation show that the movement rate is different for different terranes. The regional strike derived from decomposition analysis for different frequency band and seismic anisotropy indicate that the crust and upper mantle move separately instead of as a whole. There are material flow in the eastern and southeastern Tibetan Plateau. Therefore, the faults, the crustal and upper mantle weak layers are three different boundaries for relatively movement. Those results support the "two layer wedge plates" geodynamic model on Tibetan formation and evolution.

黄骅坳陷高尘头地区储层非均质性及成藏研究

Gaochentou region is located in the southwest direction of Gaochentou village in Huanghua city of Hebei province. In regionally structural position, It lies in Qikou sag In the middle part of Huanghua depression, which belongs to the east part of the south Dagang structure zone in the middle part of Huanghua depression. Its' very beneficial at regional structure in Gaochentou , and It becomes the advantage area for oil and gas gathered and preserved, Sandstone reservoir of Dongying Formation is main bearing bed .Dongying Formation in Gaochentou region of Huanghua depression is consisted of set of mudstone and sandstone interbeds by deposited delta fades . Dongying Formation can be divided into 3 members from above to below: the first member of Dongying Formation (FMDF), the second member of Dongying Formation (SMDF), and third member of Dongying Formation (TMDF). The lithology of the upper part of FMDF was consisted of mostly middle-grained and fine-grained sandstone, and it is small for the oil-bearing area of the sand bodies .The lithology of the lower part is coarse-grained sandstone bodies which are well connected between sandstone bodies of wells, and the lower part was main bed of oil production in Dongying Formation; SMDF and TMDF are consisted of larger scale set of mudstone, in which the sandbodies are lenticular and pinch out quickly, and the lithology was mostly fine sandstone and silt stone, in which there are little oil and gas .Because the reservoirs in this area are largely influenced by the factors such as lithology, fault and others, and the reservoirs have the strong,heterogeneity , there exists the problem of oil-down and water-up for vertical distribution of oil and gas bearing. It is not very clearly for the three dimension distribution of sandstone , and the geology researchs is not enough. So, it can't satisfy the need of further development and production for Gaochentou oilfield.Having the key problem of oil-down and water-up and the mechanism of the reservoir for Gaochentou area, There are as follow study works, the first, is study of the high-resolution correlation of sequence stratigraphy and sedimentary microfacies. Dongying Formation was divided into three parasequence sets and each parasequence set was divided into different amount of parasequences. FMDF, as the main oil and gas producing bed, can be divided into seven parasequences. Oil and gas are discovered in six parasequences except the seventh. On the basis of study of sedimentary microfacies, the sediments of Dongying Formation are considered deposited mainly in delta front subfacies. The microfacies types of Dongying Formation are sub-water distirbutary channel, sub-water natural bank, inter distributary channel bay, distributary channel mouth dam, and delta front mat sand.Seismic facies analysis and logging-constrained inversion technique were applied by Author for transverse prediction of sandstone reservoir. Having 4 modes of interwell single sandbodies correlation technique, Author have described distribution characteristics of sandbodies, and established geological reservoir model of Gaochentou reservoir.Author presented that the reservoirs characteristic have very strong heterogeneity ,and In the section of sandstone interlayed with mudstone,the folium sandstone interlayed with each other, and the wedge shaped sandbodies pinched out in the mudstone. So the pinch-out up sandstone trap and lenticular sandstone trap are easily formed. They are most small scale overlying pinches out in the place of slope. This article applies the concept of deep basin oil to resolve reasonably the problem of which the oil is below the water in Gaochentou area. Combined with the study of sedimentary facies, reservoir and other aspects, the mechanism and patterns of deep basin oil are studied on the basis of characteristics in Gaochentou area.On the basis of the above study, the mechanism of the oil and gas' migration and accumulation in isotropic sandstone and heterogeneous sandstone are thoroughly analyzed through experiments on physical modeling. Experiments on physical modeling show that the discrepancy between sand layers with different permeability and thickness has important influence on the direction, path, and injection layer of oil's migration. At the beginning of the injection of oil and gas in high permeability sand layer, the pressure is low, the migration resistance is small, and the oil and gas are more easily displacing the water in sand. So it can act as good transformation layer or reservoir. But at the beginning of the injection of oil and gas in sand layer with low permeability, the pressure is high, the migration resistance is big, and the oil and gas are more difficultly displacing the water in sand. So it can only act as bad or worse transformation layer or reservoir. Even if it cannot act as transformation layer or reservoir, it can act as water layer or dry layer. The discrepancy between sand layers on permeability and thickness can make discrepancy in injection of oil and gas between different layers. Consequently it leads to small amount of oil and gas injection in sand layers with low permeability. Ultimately it affects the oil's accumulation and distribution in different sand layers.At Last, combining analysis of the structure and pool forming condition, The thesis has established models of reservoir formation to predict the advantage distribution of oil and gas bearing , and put forward the prospective target It is not only of theoretical signification for explosion and importance, but also has realistic value in guiding the progressive petroleum exploration and exploitation.

利用宽频带台阵资料研究华北地区地壳及上地幔间断面结构

Intense tectonic renovation has occurred in the eastern continent of china since Mesozoic, as evidenced by the high heat flow, widespread magma extrusion and volcanic activities, and development of large sedimentary basins. To explain the cause and mechanism for the tectonic process in this period, some researchers have put forward various models, such as mantle plume, subduction of the Pacific slab, Yangtze Block-North China Block collision, etc. Their seismological evidence, however, is still scarce..During the period from 2000 to 2003, large temporary seismic arrays were established in North China by the Institute of the Geology and Geophysics, Chinese Academy of Sciences. Total 129 portable seismic stations were linearly emplaced across the western and eastern boundaries of the Bohai Bay Basin, and accumulated a large amount of high-quality data. Moreover, abundant data were also collected at the capital digital seismic network established in the ninth five-year period of national economic and social development. These provide an unprecedented opportunity for us to study the deep structure and associated geodynamic mechanism of lithospheric processes in North China using seismological techniques.Seismology is a kind of observation-based science. The development of seismic observations greatly promotes the improvement of seismologic theory and methodology. At the beginning of this thesis, I review the history of seismic observation progress, and present some routine processing techniques used in the array seismology. I also introduce two popular seismic imaging methods (receiver function method and seismic tomography).Receiver function method has been widely used to study the crustal and upper mantle structures, and many relevant research results have been published. In this thesis I elaborate the theory of this method, including the basic concept of receiver functions and the methodology for data pre-processing, stacking and migration. I also address some problems often encountered in practical applications of receiver function imaging.By using the teleseismic data collected at the temporary seismic arrays in North China, in particular, the traveltime information of P-to-S conversion and multiple reverberations of the Moho discontinuity, I obtain the distributions of the crustal thickness and the poisson ratio at the northwest boundary area of the Bohai Bay Basin and discuss the geological implications of the results.Through detailed intestigations on the crustal structural feature around the middle part of the Tanlu fault, considerable disparity in poisson ratios is found in the western and eastern sides of the Tanlu fault. Moreover, an obvious Moho offset is coincidently observed at the same surface location. A reasonable density model for the Tanlu fault area is also derived by simulating the observed gravity variations. Both receiver function study and gravity anomaly modeling suggest that the crustal difference between the western and eastern sides of the Tanlu fault is mainly resulted from their different compositions.With common conversion point imaging of receiver functions, I estimate the depths of the upper and lower boundaries of the mantle transition zone, i.e., the 410 and 660 km discontinuities, beneath most part of the North China continent The thickness of the transition zone (TTZ) in the study area is calculated by subtracting the depth of .410 km discontinuity from that of the 660km discontinuity. The resultant TTZ is 10-15 km larger in the east than in the west of the study area. Phase transitions at the 410 km and the 660 km discontinuities are known to have different Clapeyron slopes. Therefore, the TTZ is sensitive to the temperature changes in the transition zone. Previous studies have shown that the TTZ would be smaller in the mantle plume areas and become larger when the remnants of subducted slabs are present The hypothesis of mantle plume cannot give a reasonable interpretation to the observed TTZ beneath North China, Instead, the receiver function imaging results favor a dynamic model that correlates the thermal structure of the mantle transition zone and associated upper mantle dynamics of North China to the Pacific plate subduction process.

双聚焦偏移速度估计与可分表示叠前弹性反演

In exploration seismology, the geologic target of oil and gas reservoir in complex medium request the high accuracy image of the structure and lithology of the medium. So the study of the prestack image and the elastic inversion of seismic wave in the complex medium come to the leading edge. The seismic response measured at the surface carries two fundamental pieces of information: the propagation effects of the medium and the reflections from the different layer boundaries in the medium. The propagation represent the low-wavenumber component of the medium, it is so-called the trend or macro layering, whereas the reflections represent the high-wavenumber component of the medium, it is called the detailed or fine layering. The result of migration velocity analysis is the resolution of the low-wavenumber component of the medium, but the prestack elastic inversion provided the resolution of the high-wavvenumber component the medium. In the dissertation, the two aspects about the migration velocity estimation and the elastic inversion have been studied.Firstly, any migration velocity analysis methods must include two basic elements: the criterion that tell us how to know whether the model parameters are correct and the updating that tell us how to update the model parameters when they are incorrect, which are effected on the properties and efficiency of the velocity estimation method. In the dissertation, a migration velocity analysis method based on the CFP technology has been presented in which the strategy of the top-down layer stripping approach are adapted to avoid the difficult of the selecting reduce .The proposed method has a advantage that the travel time errors obtained from the DTS panel are defined directly in time which is the difference with the method based on common image gather in which the residual curvature measured in depth should be converted to travel time errors.In the proposed migration velocity analysis method, the four aspects have been improved as follow:? The new parameterization of velocity model is provided in which the boundaries of layers are interpolated with the cubic spline of the control location and the velocity with a layer may change along with lateral position but the value is calculated as a segmented linear function of the velocity of the lateral control points. The proposed parameterization is suitable to updating procedure.? The analytical formulas to represent the travel time errors and the model parameters updates in the t-p domain are derived under local lateral homogeneous. The velocity estimations are iteratively computed as parametric inversion. The zero differential time shift in the DTS panel for each layer show the convergence of the velocity estimation.? The method of building initial model using the priori information is provided to improve the efficiency of velocity analysis. In the proposed method, Picking interesting events in the stacked section to define the boundaries of the layers and the results of conventional velocity analysis are used to define the velocity value of the layers? An interactive integrate software environment with the migration velocity analysis and prestack migration is built.The proposed method is firstly used to the synthetic data. The results of velocity estimation show both properties and efficiency of the velocity estimation are very good.The proposed method is also used to the field data which is the marine data set. In this example, the prestack and poststack depth migration of the data are completed using the different velocity models built with different method. The comparison between them shows that the model from the proposed method is better and improves obviously the quality of migration.In terms of the theoretical method of expressing a multi-variable function by products of single-variable functions which is suggested by Song Jian (2001), the separable expression of one-way wave operator has been studied. A optimization approximation with separable expression of the one-way wave operator is presented which easily deal with the lateral change of velocity in space and wave number domain respectively and has good approach accuracy. A new prestack depth migration algorithm based on the optimization approximation separable expression is developed and used to testing the results of velocity estimation.Secondly, according to the theory of the seismic wave reflection and transmission, the change of the amplitude via the incident angle is related to the elasticity of medium in the subsurface two-side. In the conventional inversion with poststack datum, only the information of the reflection operator at the zero incident angles can be used. If the more robust resolutions are requested, the amplitudes of all incident angles should be used.A natural separable expression of the reflection/transmission operator is represented, which is the sum of the products of two group functions. One group function vary with phase space whereas other group function is related to elastic parameters of the medium and geological structure.By employing the natural separable expression of the reflection/transmission operator, the method of seismic wave modeling with the one-way wave equation is developed to model the primary reflected waves, it is adapt to a certain extent heterogeneous media and confirms the accuracy of AVA of the reflections when the incident angle is less than 45'. The computational efficiency of the scheme is greatly high.The natural separable expression of the reflection/transmission operator is also used to construct prestack elastic inversion algorithm. Being different from the AVO analysis and inversion in which the angle gathers formed during the prstack migration are used, the proposed algorithm construct a linear equations during the prestack migration by the separable expression of the reflection/transmission operator. The unknowns of the linear equations are related to the elasticity of the medium, so the resolutions of them provided the elastic information of the medium.The proposed method of inversion is the same as AVO inversion in , the difference between them is only the method processing the amplitude via the incident angle and computational domain.

吉林金川西大甸子泥炭沉积中火山物质的发现及其意义附：泥炭中胡敏酸的稳定碳同位素值（δ<'13>C）作为古气候代用指标可行性研究

泥炭以发育广泛、连续沉积好、沉积速率大（可达1mm/a甚至更大）、记录丰富等特点,正成为过去全球变化研究的良好综合地质档案.参考冰芯、黄土、深海沉积物等中矿物固体的研究思路和方法,希望能对泥炭中可能存在的火山物质和其它风粉尘进行研究.在对吉林金川西大甸子泥炭样品的处理过程中,摸索出一套泥炭样品中固体物质（如火山物质）初步的处理、分析方法;并对金川西大甸子干玛珥泥炭沉积中部一段泥砂状物质进行鉴定,首次发现火山玻璃,判断这些泥砂物为火山灰砂.目前研究表明,泥炭纤维素的δ<'13>C序列能较好地反映古气候变化,那么泥炭中胡敏酸的δ<'13>C序列可否也能反映古气候变化,该文对此作了初步研究.用碱溶酸沉淀方法从泥炭中提取出所谓的胡敏酸,用燃烧法将胡敏酸、纤维素中的碳转化为CO<,2>气体,测定它们的稳定碳同位素组成（δ<'13>C）.与能较好地反映气候变化的纤维素δ<'13>C序列相比,胡敏酸的δ<'13>C序列反映气候变化不明显.

磷矿的形成与Rodinia超大陆裂解、生物爆发的关系-以贵州瓮安、开阳、织金磷矿床为例

贵州省分布有两期大型磷矿：震旦纪陡山沱期和早寒武纪梅树村时期。在这两期磷矿中分别发育有瓮安生物群和小壳动物梅树村生物群，这两个生物群代表了生命早期演化过程中的两次爆发事件。由此推测磷矿的形成可能与生物的爆发存在某种联系。另外，已有研究表明，磷矿层下伏的南沱冰碳岩和上覆的灯影帽碳酸盐的形成都与晚元古代末的Rodini。超大陆裂解相关，那么沉积于它们之间的磷矿也可能与Rodinia超大陆裂解相关联。研究磷矿形成、Rodinia超大陆裂解及生物爆发三者之间的关系对揭示磷矿的物质来源、成矿环境、构造条件、成因机制等有重要的理论意义；对解决全球范围晚元古代地层、几次生物爆发事件和．冰川事件的对比问题等有极其重要的理论和实际意义。本文以晚元古代陡山沱期（瓮安、开阳）和早寒武纪（织金）磷矿床为研究对象，借鉴岩相学、古生物学、元素地球化学、同位素地球化学、有机地球化学及Sm-Nd、Rb-Sr同位素定年等方法，讨论磷矿的形成与Rodini。超大陆裂解、生物爆发之间的关系。论文主要获得了以下几个方面的认识：研究了三个矿床的主要地质特征，表明寒武纪织金磷矿床的磷品位低于陡山沱期磷矿床。瓮安、织金磷矿床形成于氧化带环境，而开阳磷矿床形成于亚氧化带环境。确定了三个矿床的成矿时代。同位素年代学得出陡山沱组磷矿床形成年龄为583±19Ma（Sm-Nd法）、588±8.6Ma和582±7.7Ma（Rb-Sr法），寒武纪磷矿床的形成年龄为542±23Ma（Sm-Nd法）和541±12Ma（Rb-Sr法）。震旦纪陡山沱期和寒武纪梅树村期磷矿床的形成时代与华南地块从Rodinia超大陆中裂解漂移出来的时间（600Ma～550Ma、540Ma）相耦合，为Rodinia超大陆裂解事件与磷矿床的形成存在密切成因联系提供了年代学方面的证据。对三个磷矿床进行的大量地质地球化学研究表明：Rodinia超大陆裂解为磷矿的形成提供了热液源和慢源成矿物质，并造成了大规模的低温热液成矿作用。研究了发育于陡山沱组和寒武纪含磷地层中的瓮安生物群和梅树村小壳生物群爆发的原因，表明Rodini。超大陆裂解造成了两次生物大爆发。超大陆裂解引发地球化学异常，导致了生物的大发展。讨论了地层中生物与磷的关系。地质证据、伴生微量元素证据、磷灰石矿物地球化学证据及有机地球化学证据表明，生物有机质直接或间接地参与了成磷作用。其中生物直接参与成磷作用主要发生在成矿过程的磷质吸取阶段，同时也见于地球化学富集阶段还原型细菌活跃的成岩带中；有机质间接成磷作用主要发生在磷块岩成矿过程的地球化学富集阶段和物理富集阶段。总结了磷矿的形成与Rodinia超大陆裂解、生物爆发之间的关系，Rodinia超大陆裂解为磷矿的形成提供了热液源和慢源的成矿物质；同时引发了生物的大爆发，而生物大爆发则促进了成矿物质的超常聚集和大规模成矿。

Gas phase catalytic partial oxidation of toluene in a microchannel reactor

Gas phase partial oxidation of toluene over V/Ti oxide catalysts has been successfully performed in a microchannel reactor, which provides very good mass and heat transfer conditions. With the elimination of hot spots, which are known as the most negative factors for partial oxidation of hydrocarbons, steady and uniform reaction conditions can be achieved in the catalyst bed by using, the microreactor. Since the best performance of the catalysts might be exploited, the selectivity of partial oxidation products of toluene has remarkably increased compared to the traditional packed fixed-bed reactor, even without the bother of modifying the catalysts, diluting the reactants or catalysts with inert contents to avoid hot spots or improve the diffusion and mixing. Furthermore, in virtue of its inherent safety features, when using pure oxygen as oxidant, the reactions were handled safety within the explosion limits in the microreactor. With TiO2 carried V2O5 as catalysts, the total selectivity of benzaldehyde and benzoic acid reaches around 60%, and the toluene conversion is about 10%. The conversion can go up without violent decline of selectivity, unlike most fixed bed reactors. Space time yield of 3.12 kg h(-1) L-1 calculated on the basis of the channel volume has been achieved. The influence of operating conditions has been investigated in detail in the microreactor. (c) 2005 Elsevier B.V. All rights reserved.

Cluster-assistant generation of multiply charged atomic ions in nanosecond laser ionization of seeded methyl iodide beam

The photoionization of methyl iodide beam seeded in argon and helium is studied by time-of-flight mass spectrometry using a 25 ns, 532 nm Nd-YAG laser with intensities in the range of 2 x 10(10)-2 x 10(11) W/cm(2). Multiply charged ions Of Iq+ (q = 2-3) and C2+ with tens of eV kinetic energies have been observed when laser interacts with the middle part of the pulsed molecular beam, whose peak profiles are independent on the laser polarization directions. Strong evidences show that these ions are coming from the Coulomb explosion of multiply charged CH3I clusters, and laser induced inverse bremsstrahlung absorption of caged electrons plays a key role in the formation of multiply charged ions. (C) 2004 Elsevier B.V. All rights reserved.