俯冲带深部流体特征——以西南天山变质俯冲杂岩为例

The subduction zone is an important site of the fluid activity and recycling of chemical elements. The fluid characteristic of deep subduction zones is a top scientific problem attracting the petrologists, geochemists and tectonists. In this dissertation, the characteristics of fluid activity within a deep subduction zone have been explored on the basis of the studies on the petrography, mineral chemistry, fluid inclusions, geochemistry and metamorphic P–T conditions of the omphacite-bearing high-pressure veins and related hosts from the low-temperature/high-pressure metamorphic belt in southwestern Tianshan, China. Multiple high-pressure veins are exposed in host eclogites and blueschists. The veins are composed predominantly of omphacite, garnet, quartz, and other minerals. Some veins contain cm-sized rutiles. In general, the vein can be divided into three types, the ‘in situ dehydration’ vein, the ‘external transport’ vein and the ‘composite’ vein. The omphacites within the veins and related host rocks contain lots of two-phase or three-phase primary fluid inclusions. The final melting temperature (Tfm) of fluid inclusions varies mainly from -0.6 to -4.3 °C, the homogeneous temperature (Th) varies from 185 to 251 °C, the salinity varies from 1.1 to 6.9 wt.% NaCl equivalent and the density varies from 0.81 to 0.9 g/cm3. The fluids were released under the conditions of T = 520–580°C and P = 15–19 kbar at blueschist facies to eclogite facies transition. The fluids include not only Li, Be, LILE, La, Pb-enriched and HFSE- and HREE-depleted aqueous fluids but also HFSE (Ti-Nb-Ta)-rich aqueous fluids. The complex composed of aluminosilicate polymers and F was the catalyst which had caused the Ti-Nb-Ta to be dissolved into the fluids. During the transport of the LILE-rich and HFSE- and HREE-poor fluids, they can exchange some chemical elements with country rocks and leach some trace elements in some extent. The rutile could be precipitated from the HFSE (Ti-Nb-Ta)-rich aqueous fluids when CO2 was added into the fluids. The host rocks could obtain some elements, such as Ca, Cs, Rb, Ba and Th, from the external fluids. The fluids with complex composition had been released within the deep subduction zone (>50 km) in Early Carboniferous during the subduction of the South Tianshan Ocean under the Yili–Central Tianshan Plate. The results obtained in this dissertation have made new progress compared with the published data (e.g. Tatsumi, 1989; Becker et al., 1999; Scambelluri and Philippot, 2001; Manning, 2004; Hermann et al., 2006; Spandler and Hermann, 2006).

复杂介质中地震波场数值模拟研究——基于广义正交多项式迭积微分算子法

In this paper, we propose a new numerical modeling method – Convolutional Forsyte Polynomial Differentiator (CFPD), aimed at simulating seismic wave propagation in complex media with high efficiency and accuracy individually owned by short-scheme finite differentiator and general convolutional polynomial method. By adjusting the operator length and optimizing the operator coefficient, both global and local informations can be easily incorporated into the wavefield which is important to invert the undersurface geological structure. The key issue in this paper is to introduce the convolutional differentiator based on Forsyte generalized orthogonal polynomial in mathematics into the spatial differentiation of the first velocity-stress equation. To match the high accuracy of the spatial differentiator, this method in the time coordinate adopts staggered grid finite difference instead of conventional finite difference to model seismic wave propagation in heterogeneous media. To attenuate the reflection artifacts caused by artificial boundary, Perfectly Matched Layer (PML) absorbing boundary is also being considered in the method to deal with boundary problem due to its advantage of automatically handling large-angle emission. The PML formula for acoustic equation and first-order velocity-stress equation are also derived in this paper. There is little difference to implement the PML boundary condition in all kind of wave equations, but in Biot media, special attenuation factors should be taken. Numerical results demonstrate that the PML boundary condition is better than Cerjan absorbing boundary condition which makes it more suitable to hand the artificial boundary reflection. Based on the theories of anisotropy, Biot two-phase media and viscous-elasticity, this paper constructs the constitutive relationship for viscous-elastic and two-phase media, and further derives the first-order velocity-stress equation for 3D viscous-elastic and two-phase media. Numerical modeling using CFPD method is carried out in the above-mentioned media. The results modeled in the viscous-elastic media and the anisotropic pore elastic media can better explain wave phenomena of the true earth media, and can also prove that CFPD is a useful numerical tool to study the wave propagation in complex media.

内蒙古白云鄂博超大型REE-Nb-Fe矿床脉状矿化时代与成矿流体特征

Bayan Obo REE-Nb-Fe ore deposit is the largest REE deposit in the world. Owing to its unique type and tremendous economic value, this deposit has widely attracted interests from geological researchers and vast amount of scientific data have been accumulated. However, its genesis, especially ore-forming age and REE sources, have been under dispute for a long time. On the basis of previous research works, this paper mainly conducts studies on the Early Paleozoic ore-forming event in the Bayan Obo deposit. The following results and conclusions can be suggested: Sm-Nd isotopic analytical results of bastnaesite, beloeilite, albite and fluorite samples from a coarse-crystalline ore lode present an isochron age of 436±35Ma. Besides, Rb-Sr isotope dating of the coarse-crystalline biotite lode that intruded into banded ores gives an isochron age of 459±39Ma. The two ages verify the exist of Early Paleozoic ore-forming event at Bayan Obo, which characterized by extensive netted mineralization of REE fluorocarbonates, aeschynite and monazite, accompanied by widely fluorite-riebeckite-aegirine-apatite alteration. Sr-Nd isotope composition of vein minerals is located between EMI and ancient lower crust component in the ISr(t)-εNd(t) correlation diagram, indicating that there is a crustal contamination during veined mineralization. A large area late Paleozoic granitoids are distributed in the southeast region of east open pit of the mine. The granitoids intruded directly into the ore-bearing dolomite, and produced intense skarnization. Moreover, at 650-660m of the drill core on 22 line and 1598m level flat in the south of East Open Pit, we firstly found skarnization rocks. Single grain and low background Rb-Sr isochrone dating on phlogopite in skarn gives 309±12Ma. Considering the intruded contacting relationship, the late Paleozoic granitoids, already extended to the under part of REE ore bodies, must be posterior to the latest intense REE mineralization, and is only a destructive tectonic and magmatic activity. Fluid inclusion types of fluorite in the Bayan Obo deposit consist of multiphase daughter mineral-bearing inclusion, two or three phase CO2-bearing inclusion and two phase aqueous inclusion. Petrography, laser Raman analysis and microthermometry study indicate that the fluids involving in REE-Nb-Fe mineralization at Bayan Obo might be mainly of H2O-CO2-NaCl-(F-REE) system. The presence of REE-carbonate as a daughter mineral in fluid inclusions shows that the original ore-forming fluids are rich in REE elements.

复杂介质边界元法数值模拟研究

A major impetus to study the rough surface and complex structure in near surface model is because accuracy of seismic observation and geophysical prospecting can be improved. Wave theory study about fluid-satuated porous media has important significance for some scientific problems, such as explore underground resources, study of earth's internal structure, and structure response of multi-phase porous soil under dynamic and seismic effect. Seismic wave numerical modeling is one of the effective methods which understand seismic propagation rules in complex media. As a numerical simulation method, boundary element methods had been widely used in seismic wave field study. This paper mainly studies randomly rough surface scattering which used some approximation solutions based on boundary element method. In addition, I developed a boundary element solution for fluid saturated porous media. In this paper, we used boundary element methods which based on integral expression of wave equation to study the free rough surface scattering effects of Kirchhoff approximation method, Perturbation approximation method, Rytov approximation method and Born series approximation method. Gaussian spectrum model of randomly rough surfaces was chosen as the benchmark model. The approximation methods result were compared with exact results which obtained by boundary element methods, we study that the above approximation methods were applicable how rough surfaces and it is founded that this depends on and ( here is the wavenumber of the incident field, is the RMS height and is the surface correlation length ). In general, Kirchhoff approximation which ignores multiple scatterings between any two surface points has been considered valid for the large-scale roughness components. Perturbation theory based on Taylor series expansion is valid for the small-scale roughness components, as and are .Tests with the Gaussian topographies show that the Rytov approximation methods improves the Kirchhoff approximation in both amplitude and phase but at the cost of an extra treatment of transformation for the wave fields. The realistic methods for the multiscale surfaces come with the Born series approximation and the second-order Born series approximation might be sufficient to guarantee the accuracy of randomly rough surfaces. It could be an appropriate choice that a complex rough surface can be divided into large-, medium-, and small-scale roughness components with their scattering features be studied by the Kirchhoff or Rytov phase approximations, the Born series approximation, and the perturbation theory, respectively. For this purpose, it is important to select appropriate parameters that separate these different scale roughness components to guarantee the divided surfaces satisfy the physical assumptions of the used approximations, respectively. In addition, in this paper, the boundary element methods are used for solving the porous elastic wave propagation and carry out the numerical simulation. Based on the fluid-saturated porous model, this paper analyses and presents the dynamic equation of elastic wave propagation and boundary integral equation formulation of fluid saturated porous media in frequency domain. The fundamental solutions of the elastic wave equations are obtained according to the similarity between thermoelasticity and poroelasticity. At last, the numerical simulation of the elastic wave propagation in the two-phase isotropic media is carried out by using the boundary element method. The results show that a slow quasi P-wave can be seen in both solid and fluid wave-field synthetic seismograms. The boundary element method is effective and feasible.

内蒙古乌拉特中旗北部宝音图群研究

The Baoyintu Group, lies in Wulate-zhongqi, Inner Mongolia, is a set of medium-grade metamorphic rock series which undergoes complex deformations. It consists of pelite schist, greenschist, plagioclase amphibolite quartzite and marble. The pelite schist is the main rock type and contains the classic medium pressure metamorphic minerals. The author divided Baoyintu group into five assemblages, investigated the rock association and plotted geological section of each assemblage in this area. Based on the systemically study of structural geology, petrology, geochemistry and mineralogy, the author reconstructs the protolith, sedimentary environment and tectonic evolution, discusses the mesoscopic and microscopic structure, metamorphism, geochemistry characters and the correlation between porphyroblast growth and deformation-metamorphism. There are three phase deformations in the research area: the earliest one occurred as the Baoyintu group deformed and metamorphosed and the main structure pattern is tight fold within layers during the Dl, large scale reversed fold and two phase faults (Fl fault and F2 fault) during the D2, and superimposed fold and F3 fault during D3. The F3 trancate the Wenduermian group of Silurian. The second and third phase deformation are relate to the orogenic event of late period of early Paleozoic. According to the rock association ,characteristics of the rocks and research of geochemistry, we get some information of the sedimentary environment and tectonic evolution of Baoyintu group. The source rocks are a set of terrigenous deposits-volcanic formation which reflect the history of the tectonic setting: stable- active-restable. And there are two sedimentary cycle from first assemblage to fifth assemblage: from first assemblage to fourth assemblage is a course of progression and the fifth assemblage is a start of regression. We also get the information of the P-T-t path by studying petrographies and calculating temperature and pressure. The path is not similiat to any classic type. And the interpretation is different from the traditional opinion. The P-T-t path reflects the dynamic course of convergence and uplift, magma underplating, back-arc extension and convergence of continental margin. Applying the theory of deformation partitioning to this area, the author discuss the relationship between deformation and porphyroblast growth, and get the conclusion of the sequence of deformation and metamorphism. At the first time we measure the distribution of chemical composition within the porphyroblast by XRF, confirm the theory of deformation partitioning quantitative and get new understanding about growth phase of porphyroblast and growth mode of porphyroblast: porphyroblast grow in the manner as "rose flower", the growth is controlled by the deformation. The elements distribution in porphyroblast reflects the growth manner and indicate history of metamorphism and deformation. So, we can deduce the metmorphism and deformation from the elements distribution in porphyroblast.

从字形到语义：跨语言Stroop效应的ERP研究

According to the influential dual-route model of reading (Coltheart, Rastle et al. 2001), there are two routes to access the meaning of visual words: one directly by orthography (orthography-semantic) and the other indirectly via the phonology (phonology-semantic). Because of the dramatic difference between written Chinese and alphabetical languages, it is still on debate whether Chinese readers have the same semantic activation processes as readers of alphabetical languages. In this study, the semantic activation processes in alphabetical German and logographic Chinese were compared. Since the N450 for incongruent color words in the Stroop tasks was induced by the semantic conflict between the meaning of the incongruent color words and color naming, this component could be taken as an index for semantic activation of incongruent color words in Stroop tasks. Two cross-script Stroop experiments were adopted to investigate the semantic activation processes in Chinese and German. The first experiment focused on the the role of phonology, while the second one focused on the realative importance of orthography. Cultural differences in cognitive processing between individuals in western and eastern countries have been found (Nisbett & Miyamoto, 2005). In order to exclude potential differences in basic cognitive processes like visual discrimination capabilities during reading, a visual Oddball experiment with non-lexical materials was conducted with all participants. However, as indicated by the P300 elicited by deviant stimuli in both groups, no group difference was observed. In the first Stroop experiments, color words (e.g., “green”), color-word associates (e.g., “grass”), and homophones of color words were used. These words were embedded into color patches with either congruent color (e.g. word “green” in green color patch) or incongruent colors (e.g. word “green” in either red or yellow or blue color patch). The key point is to observe whether homophones in both languages could induce similar behavioral and ERP Stroop effects to that induced by color words. It was also interesting to observe to which extent the N450 was related to the semantic conflicts. Nineteen Chinese adult readers and twenty German adult readers were asked to respond to the back color of these words in the Stroop experiment in their native languages by pressing the corresponding keys. In the behavioral data, incongruent conditions (incongruent color words, incongruent color-word associates, incongruent homophones) had significantly longer reaction times as compared to corresponding congruent conditions. All incongruent conditions in the Geman group elicited an N450 in the 400 to 500 ms time window. In the Chinese group, the N450 in the same time window was also observed for the incongruent color words and incongruent color-word associates. These results indicated that the N450 was very sensitive to semantic conflict-even words with semantic association to colors (e.g. “grass”) could elicite similar N450. However, the N450 was absent for incongruent homophones of color words in the Chinese group. Instead, in a later time window (600-800 ms), incongruent homophones elicited a positivity over left posterior regions as compared to congruent homophones. Similar positivity was also observed for color words in the 700 to 1000 ms time window in the Chinese group and 600 to 1000 ms time window for incongruent color words and homophones in the Geman group. These results indicate that phonology plays an important role in Geman semantic activation processes, but not in Chinese. In the second Stroop experiment, color words and pseudowords which had similiar visual shape to color words in both languages were used as materials. Another group of eighteen Chinese and twenty Germans were involved in the Stroop experiment in their native languages.The ERPs were recorded during their performance. In the behavioral data, strong and comparable Stroop effects (as counted by substract the reaction times in the congruent conditions from reaction times in the incongruent conditions) were observed. In the ERP data, both incongruent color words and incongruent pseudowords elicited an N450 over the whole brain scalp in both groups. These results indicated that orthography played an equally important role in semantic activation processes in both languages. The results of the two Stroop experiments support the view that the semantic activation process in Chiense readers differs significantly from that in German readers. The former rely mainly on the direct route (orthography-semantic), while the latter use both direct route and incirect route (phonology-semantic). These findings also indicate that the characteritics of different languages shape the semantic activation processes.