东喜马拉雅构造结新生代碰撞变形与岩石圈演化

The Eastern Himalayan Syntaxis (EHS) is one of the strongest deformation area along the Himalayan belt resulted from the collision between Indian plate and the Eurasian Plate since the 50～60Ma, and has sensitivity tracked and preserved the whole collisional processes. It should depend on the detail geological investigations to establish the deformational accommodate mode, and the uplift history, to elucidate the deep structure and the crust-mantle interaction of the Tibet Plateau of the EHS. The deep-seated (Main Mantle Thrusts) structures were exhumed in the EHS. The MMT juxtapose the Gangdese metamorphic basement and some relic of Gangdese mantle on the high Himalayan crystalline series. The Namjagbawa group which is 1200～1500Ma dated by U/Pb age of zircon and the Namla group which is 550Ma dated by U/Pb age of zircon is belong to High Himalayan crystalline series and Gangdese basement respectively. There is some ophiolitic relic along the MMT, such as metamorphic ocean mantle peridotite and metamorphic tholeiite of the upper part of ocean-crust. The metamorphic ocean mantle peridotites (spinel-orthopyroxene peridotite) show U type REE patterns. The ~(87)Sr/~(86)Sr ratios were, 0.709314～0.720788, and the ~(143)Nd/~(144)Nd ratios were 0.512073～0.512395, plotting in the forth quadrant on the ~(87)Sr/~(86)Sr-~(143)Nd/~(144)Nd isotope diagram. Some metamorphic basalt (garnet amphibolite) enclosures have been found in the HP garnet-kynite granulite. The garnet amphibolites can be divided two groups, the first group is deplete of LREE, and the second group is flat or rich LREE, and their ~(87)Sr/~(86)Sr, ~(143)Nd/~(144)Nd ratios were 0.70563～0.705381 and 0.512468～0.51263 respectively. Trace element and isotopic characteristics of the garnet amphibolites display that they formed in the E-MORB environment. Some phlogolite amphibole harzburgites, which exhibit extensive replacement by Phl, Amp, Tc and Dol etc, were exhumed along the MMT. The Phl-Amp harzburgites are rich in LREE and LILE, such as Rb, K etc, and depletes Eu (Eu~* = 0.36 ～ 0.68) and HFSE, such as Nb, Ta, Zr, Hf, P, Ti etc. The trace element indicate that the Phl-Amp harzburgites have island arc signature. Their ~(87)Sr/~(86)Sr are varied from 0.708912 to 0.879839, ~(143)Nd/~(144)Nd from 0.511993 to 0.512164, ε Nd from- 9.2 to - 12.6. Rb/Sr isochrone age of the phlogolite amphibole harzburgite shows the metasomatism took place at 41Ma, and the Amp ~(40)Ar/~(39)Ar cooling age indcate the Phl-Amp harzburgite raising at 16Ma. There is an intense crust shortening resulted from the thrust faults and folds in the Cayu block which is shortened more 120km than that of the Lasha block in 35～90Ma. With the NE corner of the India plate squash into the Gangdese arc, the sinistral Pai shear fault and the dextral Aniqiao shear fault on the both sides of the Great bent of Yalun Zangbu river come into active in 21～26Ma. On the other hand, the right-lateral Gongrigabu strike-slip faults come into activity at the same period, a lower age bound for the Gongrigabu strike-slip fault is estimated to be 23～24Ma from zircon of ion-probe U/Pb thermochronology. The Gongrigabu strike-slip faults connect with the Lhari strike-slip fault in the northwestern direction and with the Saganing strike-slip at the southeastern direction. Another important structure in the EHS is the Gangdese detachment fault system (GDS) which occurs between the sedimental cover and the metamorphic basement. The lower age of the GDS is to be 16Ma from the preliminary 40Ar/39Ar thermochronology of white mica. The GDS is thought to be related to the reverse of the subducted Indian crust and the fast uplift of the EHS. Structural and thermochronology investigation of the EHS suggest that the eastern Tibet and the western Yunnan rotated clockwise around the EHS in the period of 35～60Ma. Later, the large-scale strike-slip faults (RRD, Gaoligong and Saganing fault) prolongate into the EHS, and connect with the Guyu fault and Gongrigabu fault, which suggest that the Indianchia block escape along these faults. Two kind of magmatic rocks in the EHS have been investigated, one is the mantle-derived amphibole gabbro, dioposide diorite and amphibole diorite, another is crust origin biotit-garnet adamellite, biotit-garnet granodiorite and garnet-amphibole-biotite granite. The amphibole gabbro dioposite diorite and amphibole diorite are rich in LREE, and LILE, such as Ba, Rb, Th, K, Sr etc, depleted in HFSE, such as Nb, Ta, Zr, Hf, Ti etc. The ratio of ~(87)Sr/~(86)Sr are from 0.7044 to 0.7048, ~(143)Nd/~(144)Nd are from 0.5126 to 0.5127. The age of the mantle origin magamatic rocks, which result from the partial melt of the raising and decompression anthenosphere, is 8Ma by ~(40)Ar/~(39)Ar dating of amphibole from the diorite. The later crust origin biotite-garnet adamellite, biotite-garnet granodiorite and garnet-amphibole-biotite granite are characterized by aboudance in LREE, and strong depletion of Eu. The ratios of ~(87)Sr-~(86)Sr are from 0.795035 to 0.812028, ~(143)Nd/~(144)Nd from 0.51187 to 0.511901. The ~(40)Ar/~(39)Ar plateau age of the amphibole from the garnet-amphibole-biotite granite is 17.5±0.3Ma, and the isochrone age is 16.8±0.6Ma. Their geochemical characteristics show that the crust-derived magmatic rocks formed from partial melting of the lower curst in the post-collisional environment. A group of high-pressure kaynite-garnet granulites and enclave of high-pressure garnet-clinopyroxene grnulites and calc-silicate grnulites are outcroped along the MMT. The peak metamorphic condition of the high-pressure granulites yields T=800～960 ℃, P=1.4～1.8Gpa, corresponding the condition of 60km depth. The retrograde assemblages of the high-pressure grnulites occur at the condition of T=772.3～803.3 ℃, P=0.63～0.64Gpa. The age of the peak metamorphic assemblages are 45 ～ 69Ma indicated by the zircon U/Pb ion-plobe thermochronology, and the retrograde assemblage ages are 13～26Ma by U/Pb, ~(40)Ar/~(39)Ar thermochronology. The ITD paths of the high-pressure granulites show that they were generated during the tectonic thickening and more rapid tectonic exhumation caused by the subducting of the Indian plate and subsequent break-off of the subducted slab. A great deal of apatite, zircon and sphene fission-track ages, isotopic thermochronology of the rocks in the EHS show that its rapid raising processes of the EHS can be divided into three main periods. There are 35～60Ma, 13～25Ma, 0～3Ma. 3Ma is a turn in the course of raising in the EHS which is characterized by abruptly acceleration of uplifting. The uplift ratios are lower than 1mm .a~(-1) before 3Ma, and higher than 1mm .a~(-1) with a maximum ratio of 30mm .a~(-1) since 3Ma. The bottom (knick point) of the partial anneal belt is 3.8km above sea level in the EHS, and correspond to age of 3Ma determined by fission-track age of apatite. The average uplift ratio is about 1.4 mm .a~(-1) below the knick point. The EHS has raised 4.3km from the surface of 2.36km above sea level since 3Ma estimated by the fossil partial anneal belt of the EHS. We propose a two-stage subduction model (B+A model) basing on Structural, thermochronological, magmatical, metamorphic and geophysical investigations of the EHS. The first stage is the subduction of the Indian continental margin following after the subduction of the Tethys Ocean crust and subsequent collision with the Gangdese arc, and the second stage is the Indian crust injecting into the lower crust and upper mantle of the Tibet plateau. Slab break-off seems to be occurred between these two stages.

Experimental investigation on the cavity-based scramjet model

The present work focused on improving the engine performance with different fuel equivalence ratios and fuel injections. A scramjet model with strut/cavity integrated configurations was tested under Mach 5.8 flows. The results showed that the strut may sreve as an effective tool in a kerosene-fueled scramjet. The integration of strut/cavities also had great effect on stablizing the combustion in a wide range of fuel equivalence ratio. The one-sdimensional analysis method was used to analyze the main characteristics of the model. The two-stage fuel injection should have better performance in increasing the chemical reaction rate in the first cavity region.

Characteristics of a CW HF chemical-laser calculated from a simplified two-dimensional model

A two-dimensional simplified model of an HF chemical laser is introduced. Using an implicit finite difference scheme, the solution of two adjacent parallel streams with diffusion mixing and chemical reaction is generated. A contour of mixing and reaction boundary is obtained without presupposition. The distribution of the HF(v) concentrations, gas temperature and the optical small signal gain (alpha sub V, J) on the flowing plane (X, Y) are presented. Compared with the solution solved directly from a set of Navier-Stokes equations, the results of these two methods agree with each other qualitatively. The influences of the different velocity, temperature (T sub 0) and composition of the two streams on the small signal gain after the nozzle exit are investigated. It is interesting that for larger J with a fixed v, the peaks of alpha sub v-T sub 0 profiles move towards higher T sub 0. The computing method is simple and only a short computing time is needed.

Performance of a Supersonic Model Combustor using Vaporized Kerosene Injection

Characteristics of supersonic combustion by injecting kerosene vapor into a Mach 2.5 crossflow at various preheat temperatures and pressures were investigated experimentally. A two-stage heating system has been designed and tested, which can prepare heated kerosene of 0.8 kg up to 820 K at pressure of 5.5 Mpa with minimum/negligible fuel coking. In order to simulate the thermophysical properties of kerosene over a wide range of thermodynamic conditions, a three-component surrogate that matches the compound class of the parent fuel was employed. The flow rate of kerosene vapor was calibrated using a sonic nozzle. Computed flow rates using the surrogate fuel are in agreement with the experimental data. Kerosene jets at various preheat temperatures injecting into both quiescent environment and Mach 2.5 crossflow were visualized. It was found that at injection pressure of 4 Mpa and preheat temperature of 550 K the kerosene jet was completely in vapor phase, while keeping almost the same penetration depth as compared to the liquid kerosene injection. Supersonic combustion tests were also carried out to compare the combustor performance for the cases of vaporized kerosene injection, liquid kerosene injection, and effervescent atomization with hydrogen barbotage, under the similar stagnation conditions. Experimental results demonstrated that the use of vaporized kerosene injection leads to better combustor performance. Further parametric study on vaporized kerosene injection in a supersonic model combustor is needed to assess the combustion efficiency as well as to identify the controlling mechanism for the overall combustion enhancement.

Effects of coherence on anisotropic electromagnetic Gaussian-Schell model beams on propagation

An analytical formula for the cross-spectral density matrix of the electric field of anisotropic electromagnetic Gaussian-Schell model beams propagating in free space is derived by using a tensor method. The effects of coherence on those beams are studied. It is shown that two anisotropic stochastic electromagnetic beams that propagate from the source plane z = 0 into the half-space z > 0 may have different beam shapes (i.e., spectral density) and states of polarization in the half-space, even though they have the same beam shape and states of polarization in the source plane. This fact is due to a difference in the coherence properties of the field in the source plane. (C) 2007 Optical Society of America.

Performance of supersonic model combustors with staged injections of supercritical aviation kerosene

Supersonic model combustors using two-stage injections of supercritical kerosene were experimentally investigated in both Mach 2.5 and 3.0 model combustors with stagnation temperatures of approximately 1,750 K. Supercritical kerosene of approximately 760 K was prepared and injected in the overall equivalence ratio range of 0.5-1.46. Two pairs of integrated injector/flameholder cavity modules in tandem were used to facilitate fuel-air mixing and stable combustion. For single-stage fuel injection at an upstream location, it was found that the boundary layer separation could propagate into the isolator with increasing fuel equivalence ratio due to excessive local heat release, which in turns changed the entry airflow conditions. Moving the fuel injection to a further downstream location could alleviate the problem, while it would result in a decrease in combustion efficiency due to shorter fuel residence time. With two-stage fuel injections the overall combustor performance was shown to be improved and kerosene injections at fuel rich conditions could be reached without the upstream propagation of the boundary layer separation into the isolator. Furthermore, effects of the entry Mach number and pilot hydrogen on combustion performance were also studied.

A novel self-consistent-field lattice model for block copolymers

We develop a self-consistent-field lattice model for block copolymers and propose a novel and general method to solve the self-consistent-field equations. The approach involves describing the polymer chains in a lattice and employing a two-stage relaxation procedure to evolve a system as rapidly as possible to a free-energy minimum. In order to test the validity of this approach, we use the method to study the microphases of rod-coil diblock copolymers. In addition to the lamellar and cylindrical morphologies, micellar, perforated lamellar, gyroid, and zigzag structures have been identified without any prior assumption of the microphase symmetry. Furthermore, this approach can also give the possible orientation of the rods in different structures.

基于MRF模型的油气储层属性随机建模方法研究

Stochastic reservoir modeling is a technique used in reservoir describing. Through this technique, multiple data sources with different scales can be integrated into the reservoir model and its uncertainty can be conveyed to researchers and supervisors. Stochastic reservoir modeling, for its digital models, its changeable scales, its honoring known information and data and its conveying uncertainty in models, provides a mathematical framework or platform for researchers to integrate multiple data sources and information with different scales into their prediction models. As a fresher method, stochastic reservoir modeling is on the upswing. Based on related works, this paper, starting with Markov property in reservoir, illustrates how to constitute spatial models for catalogued variables and continuum variables by use of Markov random fields. In order to explore reservoir properties, researchers should study the properties of rocks embedded in reservoirs. Apart from methods used in laboratories, geophysical means and subsequent interpretations may be the main sources for information and data used in petroleum exploration and exploitation. How to build a model for flow simulations based on incomplete information is to predict the spatial distributions of different reservoir variables. Considering data source, digital extent and methods, reservoir modeling can be catalogued into four sorts: reservoir sedimentology based method, reservoir seismic prediction, kriging and stochastic reservoir modeling. The application of Markov chain models in the analogue of sedimentary strata is introduced in the third of the paper. The concept of Markov chain model, N-step transition probability matrix, stationary distribution, the estimation of transition probability matrix, the testing of Markov property, 2 means for organizing sections-method based on equal intervals and based on rock facies, embedded Markov matrix, semi-Markov chain model, hidden Markov chain model, etc, are presented in this part. Based on 1-D Markov chain model, conditional 1-D Markov chain model is discussed in the fourth part. By extending 1-D Markov chain model to 2-D, 3-D situations, conditional 2-D, 3-D Markov chain models are presented. This part also discusses the estimation of vertical transition probability, lateral transition probability and the initialization of the top boundary. Corresponding digital models are used to specify, or testify related discussions. The fifth part, based on the fourth part and the application of MRF in image analysis, discusses MRF based method to simulate the spatial distribution of catalogued reservoir variables. In the part, the probability of a special catalogued variable mass, the definition of energy function for catalogued variable mass as a Markov random field, Strauss model, estimation of components in energy function are presented. Corresponding digital models are used to specify, or testify, related discussions. As for the simulation of the spatial distribution of continuum reservoir variables, the sixth part mainly explores 2 methods. The first is pure GMRF based method. Related contents include GMRF model and its neighborhood, parameters estimation, and MCMC iteration method. A digital example illustrates the corresponding method. The second is two-stage models method. Based on the results of catalogued variables distribution simulation, this method, taking GMRF as the prior distribution for continuum variables, taking the relationship between catalogued variables such as rock facies, continuum variables such as porosity, permeability, fluid saturation, can bring a series of stochastic images for the spatial distribution of continuum variables. Integrating multiple data sources into the reservoir model is one of the merits of stochastic reservoir modeling. After discussing how to model spatial distributions of catalogued reservoir variables, continuum reservoir variables, the paper explores how to combine conceptual depositional models, well logs, cores, seismic attributes production history.

Two-dimensional simulation of an electron cyclotron resonance plasma source with self-consistent power deposition

Using a refined two-dimensional hybrid-model with self-consistent microwave absorption, we have investigated the change of plasma parameters such as plasma density and ionization rate with the operating conditions. The dependence of the ion current density and ion energy and angle distribution function at the substrate surface vs. the radial position, pressure and microwave power were discussed. Results of our simulation can be compared qualitatively with many experimental measurements.

A characteristics study on the performance of a 2-stage light gas gun

In order to obtain an overall and systematic understanding of the performance of a two-stage light gas gun (TLGG), a numerical code to simulate the process occurring in a gun shot is advanced based on the quasi-one-dimensional unsteady equations of motion with the real gas effect,;friction and heat transfer taken into account in a characteristic formulation for both driver and propellant gas. Comparisons of projectile velocities and projectile pressures along the barrel with experimental results from JET (Joint European Tons) and with computational data got by the Lagrangian method indicate that this code can provide results with good accuracy over a wide range of gun geometry and loading conditions.

Two-dimensional simulation of inductive plasma sources with self-consistent power deposition

A time averaged two-dimensional fluid model including an electromagnetic module with self-consistent power deposition was developed to simulate the transport of a low pressure radio frequency inductively coupled plasma source. Comparsions with experiment and previous simulation results show, that the fluid model is feasible in a certain range of gas pressure. In addition, the effects of gas pressure and power input have been discussed.

Effects of stochastic extension in ideal microcrack system

The effects of stochastic extension on the statistical evolution of the ideal microcrack system are discussed. First, a general theoretical formulation and an expression for the transition probability of extension process are presented, then the features of evolution in stochastic model are demonstrated by several numerical results and compared with that in deterministic model.

Thermal cracking of aviation kerosene for scramjet applications

Thermal cracking of China No.3 aviation kerosene was studied experimentally and analytically under supercritical conditions relevant to regenerative cooling system for Mach-6 scramjet applications. A two-stage heated tube system with cracked products collection/analysis was used and it can achieve a fuel temperature range of 700-1100 K, a pressure range of 3.5-4.5 MPa and a residence time of approximately 0.5-1.3 s. Compositions of the cracked gaseous products and mass flow rate of the kerosene flow at varied temperatures and pressures were obtained experimentally. A one-step lumped model was developed with the cracked mixtures grouped into three categories: unreacted kerosene, gaseous products and residuals including liquid products and carbon deposits. Based on the model, fuel conversion on the mass basis, the reaction rate and the residence time were estimated as functions of temperature. Meanwhile, a sonic nozzle was used for the control of the mass flow rate of the cracked kerosene, and correlation of the mass flow rate gives a good agreement with the measurements.

2-DELTA-FUNCTION GENERALIZED PLASMA POLE MODEL FOR 1ST PRINCIPLE CALCULATIONS OF QUASI-PARTICLE ENERGIES IN MANY-ELECTRON SYSTEMS

To evaluate the dynamical effects of the screened interaction in the calculations of quasiparticle energies in many-electron systems a two-delta-function generalized plasma pole model (GPP) is introduced to simulate the dynamical dielectric function. The usual single delta-function GPP model has the drawback of over simplifications and for the crystals without the center of symmetry is inappropriate to describe the finite frequency behavior for dielectric function matrices. The discrete frequency summation method requires too much computation to achieve converged results since ab initio calculations of dielectric function matrices are to be carried out for many different frequencies. The two-delta GPP model is an optimization of the two approaches. We analyze the two-delta GPP model and propose a method to determine from the first principle calculations the amplitudes and effective frequencies of these delta-functions. Analytical solutions are found for the second order equations for the parameter matrices entering the model. This enables realistic applications of the method to the first principle quasiparticle calculations and makes the calculations truly adjustable parameter free.

Two-dimensional numerical simulation of the channel electron in an In0.52Al0.48As/In0.53Ga0.48As HEMT

A two-dimensional quantum model based on the solution of Schrodinger and Poisson equations is first presented for In0.52Al0.48As/In0.53Ga0.47As/InP HEMT. According to the model, the two-dimensional distributions of electron density and transverse electric field in the channel of InAlAs/InGaAs HEMT are discussed.