华北及周边地区层析成像研究

Our study deals with the high resolution body wave tomography in North china and adjacent areas(30°N-43°N，100°E-130°E), where earthquakes occurred many times in history and has a very complicated geological structure. 6870 events recorded at 273 digital seismic stations from CDSN during 1996－2002 and stations settled by Seislab of IGCAS in Bohai Bay area, including 1382 local earthquakes and 5488 teleseismic earthquakes are used in this study. In the data we used, the average number of received stations is greater than 5, the error of picking up direct arrival time is 0.1－0.5s. Before the inversion, we use Checkerboard method to confirm the reliability of result of Local events; use Restoring Resolution Test to confirm the reliability of result of teleseismic events. We also analyzed the effect of different parameters in the inversion. Based the analysis above, the model used in this paper is divided into small blocks with a dimension of 0.33°in the latitude and longitude directions and 5km、15km、30km in depth, and initial velocity model. Using pseudobending method to calculate the ray traveling path, LSQR algorithm to inverse, finally, we got the body velocity images below 25km and above 480km in this area using Joint- inversion with local events and teleseismic events. We made the conclusion at last: (1)at top zone of the south of Sichuan Basin , there exits low velocity anomalies, below 40km is the high velocity zone extend to 300km; (2) Above the 40km of Ordos block exits low velocity zone, while below 40km until 240km, the high velocity anomalies are interlaced by low velocity anomalies. Below 300km, the anomalies are unclear any more; (3) On the whole, the velocity structure below 400km on the mantle transition zone of Eastern China area shows its changes from low velocity to high velocity.

井间地震速度层析成像方法与软件系统

Crosshole Seismic tomography has been broadly studied and applied in the fields of resource exploration and engineering exploration because of its special observing manner and better resolution than normal seismic exploration. This thesis will state the theory and method of Crosshole Seismic tomography. Basing on the previous studies，the thesis studied the initial velocity model，ray-tracing method, and developed the three-dimension tomography software. All the cells that a ray passes through are of the same velocities if the paths from transmitters to receivers are straight. The cells that the each ray passes through are recorded, and rays that pass through each cell are calculated. The ray average velocity which passes through a cell is set as the cell velocity. Analogously we can make a initial node velocity model because the velocity sum is calculated on the all cells which own to a certain node, and the cell number is summed about each nodes，the ratio of the velocity sum to the all cells number is set as the node velocity. The inversion result from the initial node velocity model is better than that of the average velocity model. Ray-bending and Shortest Path for Rays (SPR) have shortcomings and limitations respectively. Using crooked rays obtained from SPR rather than straight lines as the starting point can not only avoid ray bending converging to the local minimum travel time path, but also settle the no smooth ray problem obtained by SPR. The hybrid method costs much computation time, which is roughly equal to the time that SPR expends. The Delphi development tool based on the Object Pascal language standard has an advantage of object-oriented. TDTOM (Three Dimensions Tomography) was developed by using Delphi from the DOS version. Improvement on the part of inversion was made, which bring faster convergence velocity. TDTOM can be used to do velocity tomography from the first arrival travel time of the seismic wave, and it has the good qualities of friendly user interface and convenient operation. TDTOM is used to reconstruct the velocity image for a set of crosshole data from Karamay Oil Field. The geological explanation is then given by comparing the inversion effects of different ray-tracing methods. High velocity zones mean the cover of oil reservoir, and low velocity zones correspond to the reservoir or the steam flooding layer.

数据驱动型多次波衰减方法的研究

In the practical seismic profile multiple reflections tend to impede the task of even the experienced interpreter in deducing information from the reflection data. Surface multiples are usually much stronger, more broadband, and more of a problem than internal multiples because the reflection coefficient at the water surface is much larger than the reflection coefficients found in the subsurface. For this reason most attempts to remove multiples from marine data focus on surface multiples, as will I. A surface-related multiple attenuation method can be formulated as an iterative procedure. In this essay a fully data-driven approach which is called MPI —multiple prediction through inversion (Wang, 2003) is applied to a real marine seismic data example. This is a pretty promising scheme for predicting a relative accurate multiple model by updating the multiple model iteratively, as we usually do in a linearized inverse problem. The prominent characteristic of MPI method lie in that it eliminate the need for an explicit surface operator which means it can model the multiple wavefield without any knowledge of surface and subsurface structures even a source signature. Another key feature of this scheme is that it can predict multiples not only in time but also in phase and in amplitude domain. According to the real data experiments it is shown that this scheme for multiple prediction can be made very efficient if a good initial estimate of the multiple-free data set can be provided in the first iteration. In the other core step which is multiple subtraction we use an expanded multi-channel matching filter to fulfil this aim. Compared to a normal multichannel matching filter where an original seismic trace is matched by a group of multiple-model traces, in EMCM filter a seismic trace is matched by not only a group of the ordinary multiple-model traces but also their adjoints generated mathematically. The adjoints of a multiple-model trace include its first derivative, its Hilbert transform and the derivative of the Hilbert transform. The third chapter of the thesis is the application for the real data using the previous methods we put forward from which we can obviously find the effectivity and prospect of the value in use. For this specific case I have done three group experiments to test the effectiveness of MPI method, compare different subtraction results with fixed filter length but different window length, invest the influence of the initial subtraction result for MPI method. In terms of the real data application, we do fine that the initial demultiple estimate take on a great deal of influence for the MPI method. Then two approaches are introduced to refine the intial demultiple estimate which are first arrival and masking filter respectively. In the last part some conclusions are drawn in terms of the previous results I have got.

离散裂缝诱导各向异性粘弹性地震波模拟与特征分析

The real media always attenuate and distort seismic waves as they propagate in the earth. This behavior can be modeled with a viscoelastic and anisotropic wave equation. The real media can be described as fractured media. In this thesis, we present a high-order staggered grid finite-difference scheme for 2-D viscoelastic wave propagation in a medium containing a large number of small finite length fractures. We use the effective medium approach to compute the anisotropic parameters in each grid cell. By comparing our synthetic seismogram by staggered-grid finite-difference with that by complex-ray parameter ray tracing method, we conclude that the high-order staggered-grid finite-difference technique can effectively used to simulate seismic propagation in viscoelastic-anisotropic media. Synthetic seismograms demonstrate that strong attenuation and significant frequency dispersion due to viscosity are important factors of reducing amplitude and delaying arrival time varying with incidence angle or offset. On the other hand, the amount of scattered energy not only provides an indicator of orientation of fracture sets, but can also provide information about the fracture spacing. Analysis of synthetic seismograms from dry- and fluid-filled fractures indicates that dry-filled fractures show more significant scattering on seismic wavefields than fluid-filled ones, and offset-variations in P-wave amplitude are observable. We also analyze seismic response of an anticlinal trap model that includes a gas-filled fractured reservoir with high attenuation, which attenuates and distorts the so-called bright spot.

叠前深度偏移速度模型的建立方法研究

Seismic exploration is the main tools of exploration for petroleum. as the society needs more petroleum and the level of exploration is going up, the exploration in the area of complex geology construction is the main task in oil industry, so the seismic prestack depth migration appeared, it has good ability for complex construction imaging. Its result depends on the velocity model strongly. So for seismic prestack depth migration has become the main research area. In this thesis the difference in seismic prestack depth migration between our country and the abroad has been analyzed in system. the tomographical method with no layer velocity model, the residual curve velocity analysical method based on velocity model and the deleting method in pre-processing have been developed. In the thesis, the tomographysical method in velocity analysis is been analyzed at first. It characterized with perfection in theory and diffculity in application. This method use the picked first arrivial, compare the difference between the picked first arrival and the calculated arrival in theory velocity model, and then anti-projected the difference along the ray path to get the new velocity model. This method only has the hypothesis of high frequency, no other hypothesis. So it is very effective and has high efficiency. But this method has default still. The picking of first arrival is difficult in the prestack data. The reasons are the ratio of signal to noise is very low and many other event cross each other in prestack data. These phenomenon appear strongly in the complex geology construction area. Based on these a new tomophysical methos in velocity analysis with no layer velocity model is been developed. The aim is to solve the picking problem. It do not need picking the event time contiunely. You can picking in random depending on the reliability. This methos not only need the pick time as the routine tomographysical mehtod, but also the slope of event. In this methos we use the high slope analysis method to improve the precision of picking. In addition we also make research on the residual curve velocity analysis and find that its application is not good and the efficiency is low. The reasons is that the hypothesis is rigid and it is a local optimizing method, it can solve seismic velocity problem in the area with laterical strong velocity variation. A new method is developed to improve the precision of velocity model building . So far the pattern of seismic prestack depth migration is the same as it aborad. Before the work of velocity building the original seismic data must been corrected on a datum plane, and then to make the prestack depth migration work. As we know the successful example is in Mexico bay. It characterized with the simple surface layer construction, the pre-precessing is very simple and its precision is very high. But in our country the main seismic work is in land, the surface layer is very complex, in some area the error of pre-precessing is big, it affect the velocity building. So based on this a new method is developed to delete the per-precessing error and improve the precision of velocity model building. Our main work is, (1) developing a effective tomographical velocity building method with no layer velocity model. (2) a new high resolution slope analysis method is developed. (3) developing a global optimized residual curve velocity buliding method based on velocity model. (4) a effective method of deleting the pre-precessing error is developing. All the method as listed above has been ceritified by the theorical calculation and the actual seismic data.

释光测年实验室的建立与沉积物释光测年技术研究

The main research projects reported in this paper are the establishment of a luminescence (OSL/TL) dating laboratory in The Institute of Geology and Geophysics, CAS, and studies on OSL dating technique and protocol of sediments from North China. These projects have been suggested in order to fit in with the needs of research developments in environmental changes, in particular the aridity and desertification in North China. A new luminescence dating laboratory in which there are a Rise TL/OSL-DA-15B/C reader with Sr-90 beta source, a set of Little More Tape 9022 alpha and beta irradiators, three set of Daybreak 583 intelligent alpha counters and sample preparation system has been set up in the Institute in June 2001. The courses of the establishment of a new laboratory involved a series of technical works, besides making a suitable choice of the equipment, as follows: installing and testing TL/OSL reader, calibrating the dose rate of the beta and alpha sources in the irradiators with the standard sources, testing and calibrating the count rates of the thick source alpha counting in the alpha counters with a standard sample, and then dating of the know age samples to check and examine the OSL/TL dating system. All data obtained from above calibrations and tests show that the established OSL/TL system, including the used equipment in it, can be used to determine age of the geological and archaeological samples with an error of equivalent dose (De) of less than 5%. The OSL dates of several sediment samples obtained from the system are good agreement with those from the OSL dating laboratory in Hong Kong University and ~(14)C dates within 1 - 2 standard deviations. The studies on OSL dating technique and protocol of sediment samples being in progress involve the De determinations with single aliquot regeneration (SAR) (Murray and Wintle, 2000) of the coarse grain quartz from sand dune samples and comparison of the De determinations obtained from SAR with those measured by using multiple aliquot regeneration of loess fine grains. The preliminary results from these research works are shown as follows. The very low natural equivalent dose (De) of about 0.012 - 0.03 Gy, corresponding age of less than 10 years, for BLSL (blue light stimulated luminescence) of the coarse grain quartz from modern sand dune samples in Horqin sand fields has been determined with both the SAR and multiple aliquot regeneration (MAR) techniques. This imply that the BLSL signal zeroing of the quartz could be reached before burying of the sand in Horqin sand fields. The De values and ages of the coarse grain quartz measured with SAR protocol are in good agreement with those obtained from multiple aliquot technique for the modern sand dune samples, but the errors of De from the MAR is greater than those from the SAR. This may imply that the higher precision of age determination for younger sand dune samples could be achieved with the SAR of coarse grain quartz. The MAR combining with "Australian Slide method" may be a perfect choice for De measurements of loess fine grain samples on the basis of analysis of De values obtained from the SAR and from the MAR. The former can be employed to obtain a reliable age estimate of loess sample as older as approximately SO ka BR There is a great difference between De determinations from the (post-IR) OSL of the SAR (Roberts and Wintle, 2001) and those from independent or expected estimates for the older samples. However, the age estimates obtained from the (post-IR) OSL of the SAR are mostly closed to the independent age determinations for the younger (age less than 10 ka) fine grain samples. It may be suggested that the (post-IR) OSL of the SAR protocol of the fine grain fraction would be a suitable choice to dating of the younger samples, but may be unsuitable for the older samples.

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

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.

安多地区早中侏罗世岩浆岩岩石地球化学和同位素年代学研究及其对构造演化的制约

Anduo area is located in the Central Tibet, the middle segment of the Bangonghu-Nujiang suture. Anduo Block is the northern part of Lhasa terrane. The relationships among the different geological bodies were determined during the 1: 250000 regional geological surveying. Petrography, petrologic geochemistry, isotopic geochemistry and geochronology of igneous rocks from the suture and granitoids from Anduo Block were analyzed systematically as a whole for the first time. Then, their tectonic setting and history are discussed.Anduo ophiolitic melange consists of metamorphic peridotites, cumulates, plagiogranites, sheeted dykes swarm, pillow lava and radiolarian cherts. The concentration of Cr and Ni in the metamorphic peridotites is very high, with Mg# about 0.94 ~ 0.97, higher 87Sr/86Sr and Pb isotopic ratios, and lower 143Nd/i44Nd ratio. LREE is enriched relative to HREE and positive Eu anomaly is very clear. The REE distribution curve is U shape. Nb and Ta anomalies from cumulate gabbro and sheeted dyke swarm are not clear, while that are slightly negative from pillow lava. Plagiogranite belongs to strong calc-alkaline series with high Si, middle Al, low Fe, Mg and low K contents. Eu anomaly (~ 1.23) from plagiogranites is slightly positive. The character of all components of ophiolite is similar to that of the MORB, while to some extent the ophiolite was influenced by crustal material. Anduo ophiolite formed in a mature back-arc basin. Additionally, intermediate acidity volcanic rocks within Anduo phiolite melange are island arc calc-alkline rocks related to ocean subduction.The early-middle Jurassic plutonic rocks are tonalite, granodiorite bearing-phenocryst, magaporphyritic hornblende monzogranite, magaporphyritic monzogranite, monzogranite bearing-phenocryst and syenogranite in turn. They belong to calc-alkaline series which developed from middle K to high K series temporally. REE distribution curves of all plutonic rocks are similar and parallel to each other. SREE and negative Eu anomaly values decrease. In the multi-element spider diagram, the curves of different plutons are similar to each other, but troughs of Nb, Sr, P and Ti from young plutons become more evident. This suggests that thereare some closely petrogenetic affinities among plutonic rocks which make up amagma plutonism cycle of the early-middle Jurassic. Magma source is mainly crustal,but abundant mafic microgranular enclaves within granitoids indicate that crastalmagma should be mixed with mantle-derived magma and the mantle-derived magmadecreased subsequently. Tonalite has features of I-type granite, magaporphyriticmonzogranite is transition type, and monzogranite bearing-phenocryst is S-typegranite. The characteristic of granitoids from Anduo Block suggest that the formingtectonic setting is active continental margin.Reliable zircon U-Pb SHRIMP ages are obtained in the study area firstly. Plagiogranite from the Anduo ophiolite of the Bangonghu-Nujiang suture is 175.1 Ma, and granitoids from Anduo Block is 172.6-185.4 Ma. Additionally, plagioclase from the plagiogranite dates a 40Ar/39Ar age of 144 Ma, while biotite and hornblend from granitoids of Anduo Block give a 163-165 Ma.Similar cooling ages of plagiogranite from the Anduo ophiolitic melange and granitoids from Anduo Block and the spatial distribution of the ophiolitic rocks between Anduo, Naqu, and Shainzha area suggest that bilateral subduction of the Bangonghu-Nujiang oceanic basin took place in the early-middle Jurassic. During this subduction, Anduo ophiolitic rocks were related to north subduction of the Bangonghu-Nujiang oceanic basin and Anduo back-arc basin spreading, while granitoids from Anduo Block were related to south subduction.

环渤海地区地壳和上地幔地震层析成像研究

In Circum-Bohai region (112°~124°E, 34°~42°N ), there exists rich gas-petroleum while inner-plate seismic activity is robust. Although the tectonic structure of this region is very complicated, plenty of geological, geophysical and geochemical researches have been carried out.In this paper, guided by the ideas of "One, Two, Three and Many" and "The depth controls the shallow, the regional constrains the local", I fully take advantage of previous results so as to establish a general image of this region. After collecting the arrival-time of P-wave phases of local events and tele-seismic events recorded by the stations within this region from 1966 to 2004, I process all these data and build an initial model. Then, a tomography image of crust and upper-mantle of this region is obtained. With reference to previous results, we compare the image of various depths and five cross-profiles traverse this region along different direction. And finally, a discussion and conclusion is made.The principle contents is listed as below: 1) in the first chapter, the purpose and meaning of this thesis, the advance in seismic tomography, and the research contents and blue-print is stated; 2) in the second chapter, I introduce the regional geological setting of Circum-Bohai region, describe the tectonic and evolutionary characteristics of principle tectonic units, including Bohai Bay Basin, Yanshan Fold Zone, Taihangshan Uplifted Zone, Jiao-Niao Uplifted Zone and Luxi Uplifted Zone, and primary deep faults; 3) In the third chapter, the previous geophysical researches, i.e., gravity and geomagnetic characters, geothermal flow, seismic activity, physical character of rocks, deep seismic sounding, and previous seismic tomography, are discussed; 4) in the fourth chapter, the fundamental theory and approach of seismic tomography is introduced; 5) in the fifth chapter, the technology and approaches used in this thesis, including collecting and pre-processing of data, the establishment of initial velocity model and relocation of all events; 6) in the sixth chapter, I discuss and analyze the tomography image of various depth and five cross-sections; 7)in the seventh chapter, I make a conclusion of the results, state the existing problems and possible solutions.

Mantle-derived gaseous components in ore-forming fluids of the Xiangshan uranium

The Xiangshan U deposit, the largest hydrothermal U deposit in China, is hosted in late Jurassic felsic volcanic rocks although the U mineralization post dates the volcanics by at least 20 Ma. The mineralization coincides with intrusion of local mantle-derived mafic dykes formed during Cretaceous crustal extension in South China. Ore-forming fluids are rich in CO2, and U in the fluid is thought to have been dissolved in the form of UO2 (CO3)22− and UO2 (CO3) 34− complexes. This paper provides He and Ar isotope data of fluid inclusions in pyrites and C isotope data of calcites associated with U mineralization (pitchblende) in the Xiangshan U deposit. He isotopic compositions range between 0.1 and 2.0Ra (where Ra is the 3He/4He ratio of air=1.39×10−6) and correlates with 40Ar/36Ar; although there is potential for significant 3He production via 6Li(n,α)3H(β)3He reactions in a U deposit (due to abundant neutrons), nucleogenic production cannot account for either the 3He concentration in these fluids, nor the correlations between He and Ar isotopic compositions. It is more likely that the high 3He/4He ratios represent trapped mantle-derived gases. A mantle origin for the volatiles of Xiangshan is consistent with the δ13C values of calcites, which vary from −3.5‰ to −7.7‰, overlapping the range of mantle CO2. The He, Ar and CO2 characteristics of the ore-forming fluids responsible for the deposit are consistent with mixing between 3He- and CO2-rich mantle-derived fluids and CO2-poor meteoric fluids. These fluids were likely produced during Cretaceous extension and dyke intrusion which permitted mantle-derived CO2 to migrate upward and remobilize U from the acid volcanic source rocks, resulting in the formation of the U deposit. Subsequent decay of U within the fluid inclusions has reduced the 3He/4He ratio, and variations in U/3He result in the range in 3He/4He observed with U/3He ratios in the range 5–17×103 likely corresponding to U concentrations in the fluids b0.2 ppm.

北京市中学生考试焦虑状况及其相关因素分析

In order to analyze the test anxiety of Beijing's high school students, relevant factors affecting test anxiety and how those factors have effect on the test anxiety of those high school students, a cross-sectional study had been conducted among 2,089 students randomly selected from 15 high schools in Beijing, using the Test Anxiety Scale（TAS）, Parenting Style Scale (PSS) and Academic Aggregate Graphic (AAG). Five months later, 194 students from one of those high schools were retested with TAS with students remaining anonymous. In the retest, the classes and birth dates of those students were strictly matched. This paper is composed of three parts to examine the issue. The first part examines the test anxiety among high school students in Beijing and establishes the model; the second part conducts a comparative study of issues related to test anxiety and academic performance; the third part examines the factors affecting test anxiety and establishes the model. Results showed that: 1. The reliability and validity of TAS are satisfied and can meet survey requirement. The 25-item version of TAS turned out to have equal or even better performance compared with the original 37-item version. 2. Incidence of test anxiety in high schools of Beijing: 57.9 percent of samples have an overall score at or higher 15. 55 percent of male student samples have a score at or higher than 15 while that for female student samples is 61 percent. A score of 20 refers to fairly serious test anxiety and 31.9 percent of samples have a score over 20. 28.1 percent of the male student samples have a score of over 20 while that for female student samples is 35.9 percent. 3. The effect of grade and sex and the interaction between the two factors are statistically significant. Female students have higher test anxiety than male students and the level of test anxiety varies from grade to grade. 4. Samples are divided into two groups, one with test anxiety and the other without. The academic performance gap between the two groups is very significant. There is a significant negative correlation between academic score and test anxiety. 5. There is a negative correlation between test anxiety and mastery goal orientation and a positive correlation between test anxiety and performance goal orientation. Students with higher academic self-efficiency have lower test anxiety. The more one thinks study is valuable, the lower his/her test anxiety is. Those whose parents communicate smoothly with them have lower test anxiety. Those whose parents have an inconsistent communication style have higher test anxiety. 6. Achievement goal is mediate variable for the effect of lack of values on test anxiety. 7. The indirect effect of introduction of achievement goal and parenting style, including communication, monitoring and communication consistence, on test anxiety is significant. Key words: Test anxiety，reliability, validity, self-efficiency，parenting style，achievement goal orientation, academic values, academic score

Volume-Time Curve: An Alternative for Endotracheal Tube Cuff Management

BACKGROUND: Despite numerous studies on endotracheal tube cuff pressure (CP) management, the literature has yet to establish a technique capable of adequately tilling the cuff with an appropriate volume of air while generating low CP in a less subjective way. the purpose of this prospective study was to evaluate and compare the CP levels and air volume required to fill the endotracheal tubes cuff using 2 different techniques (volume-time curve versus minimal occlusive volume) in the immediate postoperative period after coronary artery bypass grafting. METHODS: A total of 267 subjects were analyzed. After the surgery, the lungs were ventilated using pressure controlled continuous mandatory ventilation, and the same ventilatory parameters were adjusted. Upon arrival in the ICU, the cuff was completely deflated and re-inflated, and at this point the volume of air to fill the cuff was adjusted using one of 2 randomly selected techniques: volume-time curve and minimal occlusive volume. We measured the volume of air injected into the cuff, the CP, and the expired tidal volume of the mechanical ventilation after the application of each technique. RESULTS: the volume-time curve technique demonstrated a significantly lower CP and a lower volume of air injected into the cuff, compared to the minimal occlusive volume technique (P < .001). No significant difference was observed in the expired tidal volume between the 2 techniques (P = .052). However, when the subjects were submitted to the minimal occlusive volume technique, 17% (n = 47) experienced air leakage as observed by the volume-time graph. CONCLUSIONS: the volume-time curve technique was associated with a lower CP and a lower volume of air injected into the cuff, when compared to the minimal occlusive volume technique in the immediate postoperative period after coronary artery bypass grafting. Therefore, the volume-time curve may be a more reliable alternative for endotracheal tube cuff management.

Time-dependent Hyperstar algorithm for robust vehicle navigation in time-dependent stochastic road networks.

The vehicle navigation problem studied in Bell (2009) is revisited and a time-dependent reverse Hyperstar algorithm is presented. This minimises the expected time of arrival at the destination, and all intermediate nodes, where expectation is based on a pessimistic (or risk-averse) view of unknown link delays. This may also be regarded as a hyperpath version of the Chabini and Lan (2002) algorithm, which itself is a time-dependent A* algorithm. Links are assigned undelayed travel times and maximum delays, both of which are potentially functions of the time of arrival at the respective link. The driver seeks probabilities for link use that minimise his/her maximum exposure to delay on the approach to each node, leading to the determination of the pessimistic expected time of arrival. Since the context considered is vehicle navigation where the driver is not making repeated trips, the probability of link use may be interpreted as a measure of link attractiveness, so a link with a zero probability of use is unattractive while a link with a probability of use equal to one will have no attractive alternatives. A solution algorithm is presented and proven to solve the problem provided the node potentials are feasible and a FIFO condition applies for undelayed link travel times. The paper concludes with a numerical example.