Scanning tunnelling microscopy observation of cytochrome-c denaturation induced by bromopyrogal red on highly oriented pyrolytic graphite

The denaturation of cytochrome-e (cyt-c) induced by bromopyrogal red (BPR) was studied by scanning tunnelling microscopy (STM) on the electrochemically pretreated highly oriented pyrolytic graphite (HOPG) surface. STM images reveal that denatured cyt-c molecules exist in variable states including aggregates, globular compact, partially unfolded and combined with BPR molecule. The apparently low image contrast of denatured cyt-c observed in this experiment comparing to that of native cyt-c molecules, and the relative low image contrast of the unfolded part comparing with the compact globular part, are ascribed to the unfavourable tunnelling paths for the conformational variations of denatured cyt-c molecules. (C) 1997 Elsevier Science B.V.

A new topological index for the Changchun institute of applied chemistry C-13 NMR information system

A method to assign a single number representation for each atom (node) in a molecular graph, Atomic IDentification (AID) number, is proposed based on the counts of weighted paths terminated on that atom. Then, a new topological index, Molecular IDentification (MID) number is developed from AID. The MID is tested systematically, over half a million of structures are examined, and MID shows high discrimination for various structural isomers. Thus it can be used for documentation in the Changchun Institute of Chemistry C-13 NMR information system.

非完整自主车基于圆轨迹的道路避障

本文就四轮非完整自主车提出了一种基于圆轨迹的道路避障策略。先将道路上的障碍按照障碍距离自主车的远近划分层次,使一个层次的障碍能在自主车视场中全部出现。然后给出基于圆轨迹的避障算法,即自主车沿由自主车出发位姿和子目标点确定的圆弧轨迹走行。在此之前推导四轮非完整车的运动模型为提出避障策略的基础准备。尽量减小自主车在走行过程中运动状态的改变,基于圆轨迹避障策略能够很好地满足这一要求。最后引入代价函数,给出对于此方法的评价,体现了本方法的优越性。

一种多移动机器人避碰规划方法

本文采用集中预规划方法 ,通过调整机器人的运动速度实现多机器人避碰 ,所提算法的基本思想为 :将机器人的运动路径分段 ,然后按避碰要求对机器人通过各段的时间进行约束 ,从而将避碰问题转化为高维线性空间的优化问题 ,并进一步将其转化为线性方程的求解 ,使问题具有明确的解析解 .由于该方法的复杂度较高 ,在实现过程中采用了多种方法降低复杂度 ,简化计算 .本文给出了该算法的基本思路 ,有关定理及证明 ,算法的化简方法 ,最后给出了实验结果及分析 .

两轮驱动自动引导车辆定位停车研究

建立了两轮独立驱动自动引导车辆的运动学模型，计算了两轮以不同速度行驶时的运动轨迹。基于所建的运动学模型，提出了两轮驱动自动引导车的定位控制策略。通过试验证实，这种定位控制策略具有很快的停车速度和较高的定位精度。

浅谈网络电子期刊的利用

介绍了网络电子期刊的定义和特点,论述了网络电子期刊的获取途径,探讨了网络电子期刊对图书馆等信息服务机构的影响。

基于Multi-Agent的生产过程建模与仿真技术研究

某些流程行业由于采用按配方进行分组加工的模式组织生产,在排产时存在多条路径调度优化的问题,应用一般的优化算法对于现场在线调度难以给出满意结果,而基于Agent的过程仿真在解决离散、非线性系统模拟方面有显著的优势,本文采用Agent的方法对生产过程建模,然后对方案组内的备选方案进行仿真,通过对比各方案的仿真结果找到最优的方案作为执行方案,为现场的优化排产提供决策支持。

青藏高原及邻区面波频散反演--岩石圈-软流层体系结构研究

The Qinghai-Tibet Plateau lies in the place of the continent-continent collision between Indian and Eurasian plates. Because of their interaction the shallow and deep structures are very complicated. The force system forming the tectonic patterns and driving tectonic movements is effected together by the deep part of the lithosphere and the asthenosphere. It is important to study the 3-D velocity structures, the spheres and layers structures, material properties and states of the lithosphere and the asthenosphere for getting knowledge of their formation and evolution, dynamic process, layers coupling and exchange of material and energy. Based on the Rayleigh wave dispersion theory, we study the 3-D velocity structures, the depths of interfaces and thicknesses of different layers, including the crust, the lithosphere and the asthenosphere, the lithosphere-asthenosphere system in the Qinghai-Tibet Plateau and its adjacent areas. The following tasks include: (1)The digital seismic records of 221 seismic events have been collected, whose magnitudes are larger than 5.0 over the Qinghai-Tibet Plateau and its adjacent areas. These records come from 31 digital seismic stations of GSN , CDSN、NCDSN and part of Indian stations. After making instrument response calibration and filtering, group velocities of fundamental mode of Rayleigh waves are measured using the frequency-time analysis (FTAN) to get the observed dispersions. Furthermore, we strike cluster average for those similar ray paths. Finally, 819 dispersion curves (8-150s) are ready for dispersion inversion. (2)From these dispersion curves, pure dispersion data in 2°×2° cells of the areas (18°N-42°N, 70°E-106°E) are calculated by using function expansion method, proposed by Yanovskaya. The average initial model has been constructed by taking account of global AK135 model along with geodetic, geological, geophysical, receiving function and wide-angle reflection data. Then, initial S-wave velocity structures of the crust and upper mantle in the research areas have been obtained by using linear inversion (SVD) method. (3)Taking the results of the linear inversion as the initial model, we simultaneously invert the S wave velocities and thicknesses by using non-linear inversion (improved Simulated Annealing algorithm). Moreover, during the temperature dropping the variable-scale models are used. Comparing with the linear results, the spheres and layers by the non-linear inversion can be recognized better from the velocity value and offset. (4)The Moho discontinuity and top interface of the asthenosphere are recognized from the velocity value and offset of the layers. The thicknesses of the crust, lithosphere and asthenosphere are gained. These thicknesses are helpful to studying the structural differentia between the Qinghai-Tibet Plateau and its adjacent areas and among geologic units of the plateau. The results of the inversion will provide deep geophysical evidences for studying deep dynamical mechanism and exploring metal mineral resource and oil and gas resources. The following conclusions are reached by the distributions of the S wave velocities and thicknesses of the crust, lithosphere and asthenosphere, combining with previous researches. (1)The crust is very thick in the Qinghai-Tibet Plateau, varying from 60 km to 80 km. The lithospheric thickness in the Qinghai-Tibet Plateau is thinner (130-160 km) than its adjacent areas. Its asthenosphere is relatively thicker, varies from 150 km to 230 km, and the thickest area lies in the western Qiangtang. India located in south of Main Boundary thrust has a thinner crust (32-38 km), a thicker lithosphere of about 190 km and a rather thin asthenosphere of only 60 km. Sichuan and Tarim basins have the crust thickness less than 50km. Their lithospheres are thicker than the Qinghai-Tibet Plateau, and their asthenospheres are thinner. (2)The S-wave velocity variation pattern in the lithosphere-asthenosphere system has band-belted distribution along east-westward. These variations correlate with geology structures sketched by sutures and major faults. These sutures include Main Boundary thrust (MBT), Yarlung-Zangbo River suture (YZS), Bangong Lake-Nujiang suture (BNS), Jinshajiang suture (JSJS), Kunlun edge suture (KL). In the velocity maps of the upper and middle crust, these sutures can be sketched. In velocity maps of 250-300 km depth, MBT, BNS and JSJS can be sketched. In maps of the crustal thickness, the lithospheric thickness and the asthenospheric thickness, these sutures can be still sketched. In particular, MBT can be obviously resolved in these velocity maps and thickness maps. (3)Since the collision between India and Eurasian plate, the “loss” of surface material arising from crustal shortening is caused not only by crustal thickening but also by lateral extrusion material. The source of lateral extrusion lies in the Qiangtang block. These materials extrude along the JSJS and BNS with both rotation and dispersion in Daguaiwan. Finally, it extends toward southeast direction. (4)There is the crust-mantle transition zone of no distinct velocity jump in the lithosphere beneath the Qiangtang Terrane. It has thinner lithosphere and developed thicker asthenosphere. It implies that the crust-mantle transition zone of partial melting is connected with the developed asthenosphere. The underplating of asthenosphere may thin the lithosphere. This buoyancy might be the main mechanism and deep dynamics of the uplift of the Qinghai-Tibet hinterland. At the same time, the transport of hot material with low velocity intrudes into the upper mantle and the lower crust along cracks and faults forming the crust-mantle transition zone.

陕北地区石家湾-李家岔探区三叠系延长组石油成藏条件分析

Shijiawan –Lijiacha area, lying on the northeastern part of the Shanbei Slope of Ordos Basin, was selected as studying area. The previous explorations proved that the 2nd segment and 6th segment of the Yanchang Formation are the most important oil-bearing formations. It is indicated that the sedimentary facies and reservoir characteristics restricted the hydrocarbon accumulation regularity. Therefore, with petrology methodologies, such as outcrop observation, core description, geophysical logging interpretation, thin section determination, scanning electron microscope, as well as rock property analysis, the reservoirs was were systematically studied and characterized. The sedimentary micro-facies, seals, reservoir-seal combines, migration pathways and entrapping modes were taken into account. The author tempted to establish a base for further studies on reservoirs and on petroleum geology, and to provide some reliably geological evidences for later prospect activities. It was found that the sediments in the 2nd and 3rd segments of the Yanchang Formation in Shijiawan –Lijiacha area were deposited in braided rivers, and most sandy-bodies were identified as channel sandbars. The 4+5th and 6th segments were principally deposited in deltaic-plain environment, consisting of corresponding sub-facies such as distributary channels, natural levee, crevasse-splay and marsh. The skeleton sandy-bodies were identified as sandy sediments of distributary channels. The sand grains in reservoir in studied area possess generally low mineralogical maturity and moderate structural maturity, and the form of pores may be classified into intergranular types and dissolved types. Most reservoirs of Yanchang Formation in Shijiawan –Lijiacha area belong to extreme low-porosity low-permeability ones (type III), and the 2nd sediments belongs to low permeability one (type II) and the 6th segment belong to super low-permeability one(type Ⅳ). The reservoirs in the 2nd segment behave more heterogeneous than those in the 6th segment. The statistic analysis results show that, for 6th and 4+5th segments, the high quality reservoir-seal combines may be found everywhere in the studied area except in the northwest and the southwest parts; and for 1st and 2nd segments, in the northeast, central and southwest parts Petroleum migration happened in the duration of the Early Cretaceous period in both lateral and vertical directions. The migration paths were mainly constructed by permeable sandy-bodies. The superimposed channel sandy-bodies consist of the principal part of the system of carriers. the vertical fractures, that may travel through the seals between reservoirs, offered the vertical paths for migrating oil. It may be synthesized that oil coming from south kitchens migrated first laterally in carriers in the 6th segment. When arrived at the studied area, oil will migration laterally or/and vertical within both the sandy-bodies and fractures, in a climbing-stair way. The results demonstrate that the oil was entrapped in traps structure-lithology and/or lithology traps. In some cases, the hydrodynamic force may help to trap oil. Accumulation of oil in the area was mainly controlled by sedimentary facies, seals, structure, and heterogeneity of reservoir in the 2nd, 4+5th and 6th segments. Especially, the oil distributions in both the 2nd and 6th segments were obviously influenced by seals in the 4+5th segment. The existence of seals in 1st segment seems important for accumulation in the 2nd segment.

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

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.

旋转轴对称弱各向异性介质模型的射线追踪方法

Ray tracing is a rapid and effective method for wave field calculation. Not only in the field of seismic-wave theory, but also in the field of seismic inversion and migration imaging,the seismic ray tracing method has become one of the most important methods. In anisotropic media, group velocity and phase velocity have different propagation directions. The seismic wave propagates along the direction of group velocity , it does not depend on the direction of phase velocity. Ray angle is a complex function with respect to phase angle, it is difficult to measure and calculate. But most rocks are weak anisotropic, so the expression of phase velocity can be simplified greatly. Based on the approximate expression of phase velocity this thesis for rotating axisymmetric weak anisotropic media deduces an expression of the partial derivative of phase velocity and an expression of group velocity with the method of linear approximation. This paper uses the fourth order Runge-Kutta method together with the two-dimensional interpolation and linear interpolation to obtain the parameters of the physical locations. At last the paths of seismic wave in rotating axisymmetric weak anisotropic media are computed. According to the analysis of the computational results, it indicates that the method developed in this paper has strong adaptability, high computational efficiency and high accuracy for rotating axisymmetric weak anisotropic media.

波前重建法射线追踪及走时表压缩

As a fast and effective method for approximate calculation of seismic numerical simulation, ray tracing method, which has important theory and practical application value, in terms of seismic theory and seismic simulation, inversion, migration, imaging, simplified from seismic theory according to geometric seismic, means that the main energy of seismic wave field propagates along ray paths in condition of high-frequency asymptotic approximation. Calculation of ray paths and traveltimes is one of key steps in seismic simulation, inversion, migration, and imaging. Integrated triangular grids layout on wavefront with wavefront reconstruction ray tracing method, the thesis puts forward wavefront reconstruction ray tracing method based on triangular grids layout on wavefront, achieves accurate and fast calculation of ray paths and traveltimes. This method has stable and reasonable ray distribution, and overcomes problems caused by shadows in conventional ray tracing methods. The application of triangular grids layout on wavefront, keeps all the triangular grids stable, and makes the division of grids and interpolation of a new ray convenient. This technology reduces grids and memory, and then improves calculation efficiency. It enhances calculation accuracy by accurate and effective description and division on wavefront. Ray tracing traveltime table, which shares the character of 2-D or 3-D scatter data, has great amount of data points in process of seismic simulation, inversion, migration, and imaging. Therefore the traveltime table file will be frequently read, and the calculation efficiency is very low. Due to these reasons, reasonable traveltime table compression will be very necessary. This thesis proposes surface fitting and scattered data compression with B-spline function method, applies to 2-D and 3-D traveltime table compression. In order to compress 2-D (3-D) traveltime table, first we need construct a smallest rectangular (cuboidal) region with regular grids to cover all the traveltime data points, through the coordinate range of them in 2-D surface (3-D space). Then the value of finite regular grids, which are stored in memory, can be calculated using least square method. The traveltime table can be decompressed when necessary, according to liner interpolation method of 2-D (3-D) B-spline function. In the above calculation, the coefficient matrix is stored using sparse method and the liner system equations are solved using LU decomposition based on the multi-frontal method according to the sparse character of the least square method matrix. This method is practiced successfully in several models, and the cubic B-spline function can be the best basal function for surface fitting. It make the construction surface smooth, has stable and effective compression with high approximate accuracy using regular grids. In this way, through constructing reasonable regular grids to insure the calculation efficiency and accuracy of compression and surface fitting, we achieved the aim of traveltime table compression. This greatly improves calculation efficiency in process of seismic simulation, inversion, migration, and imaging.

印度与欧亚板块强烈相互作用区深部结构特征研究

The central-south Tibet is a part of the products of the continental plate collision between Eurasia and India. To study the deep structure of the study area is significant for understanding the dynamics of the continental-continental collision. A 3-D density model matched well with the observations in the central-south Tibet was proposed in this study. In addition, this study has also used numerical simulation method to prove that Quasi-Love (QL) wave is deduced by anisotropy variation but not by lateral heterogeneity. Meanwhile, anisotropy variation in the upper mantle of the Qiangtang terrane and Lhasa terrane is detected by the QL waves observed in recorded seismograms. Based on the gravity modeling, some results are summarized as follows: 1) Under the constrain of geometrical structure detected by seismic data, a 3-D density model and Moho interface are proposed by gravity inversion of the central-south Tibet. 2) The fact that the lower crustal densities are smaller than 3.2 g/cm3, suggests absence of eclogite or partial eclogitization due to delamination under the central-south Tibet. 3) Seismicity will be strong or weak in the most negative Bouguer gravity anomaly. So there is no a certain relationship between seismicity and Bouguer gravity anomaly. 4) Crustal composition are determined after temperature-pressure calibration of seismic P wave velocity. The composition of lower crust might be one or a mixture of: 1. amphibolite and greenschist facies basalt beneath the Qiangtang terrane; 2. gabbro-norite-troctolite and mafic granulite beneath the Lhasa terrane. Because the composition of the middle crust cannot be well constrained by the above data set, the data set published by Rudnick & Fountain (1995) is used for comparison. It indicated the composition of the middle crust is granulite facies and might be pelitic gneisses.Granulite facies used to be interpreted as residues of partial melting, which coincidences with the previous study on partial melting middle crust. Amphibolite facies are thought to be produced after delamination, when underplating works in the rebound of the lower crust and lithospheric mantle. From the seismology study, I have made several followed conclusions: 1) Through the numerical simulation experiment of surface wave propagating in heterogeneity media, we can find that amplitude and polarization of surface wave only change a little when considering heterogeneity. Furthermore, it is proved that QL waves, generated by surface wave scattering, are caused by lateral variation of anisotropy but not by heterogeneity. 2) QL waves are utilized to determine the variation of uppermost mantle anisotropy of the Tibetan plateau. QL waves are identified from the seismograms of the selected paths recorded by the CAD station. The location of azimuth anisotropy gradient is estimated from the group velocities of Rayleigh wave, Love wave and QL wave. It suggests that south-north lateral variation of azimuthal anisotropy locates in Tanggula mountain, and east-west lateral variation in the north of Gandese mountain with 85°E longitude and near the Jinsha river fault with 85°E longitude.

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

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.

油气运聚动力学的模拟分析——以准噶尔盆地为例

The petroleum migration, happening in the geologic past, is the very important and complex dynamic processes in the petroleum systems. It plays a linking role among all static factors in a system. The accumulation is in fact the result of the petroleum migration. For the petroleum geology, the dynamics research of the petroleum migration refers to the mechanism and process research, as well as the use of the quantitative methods. In this thesis, combining with the qualitative analysis and quantitative modeling, the author manages to discuss theoretically some key problems dealing with migration processes, which have not been solved yet, and to apply the studied results in petroleum system analysis in actual basins. The basin analysis offers the base of the numerical modeling for geological phenomena occurring in sedimentary basins, that consists of the sedimentary facies analysis, the section reconstructing technique, eroded thickness estimating, etc. The methods to construct the geologic model, which is needed in the research of oil and gas migration and accumulation, are discussed. The basin analysis offers also the possibility for the latter modeling works to get and select the parameters, such as stratum's thickness, age, stratigraphy etc. Modeling works were done by using two basin modeling softwares: Basin_Mod and TPC_Mod. The role of compaction during the secondary migration and the heterogeneity of migrating paths within the clastic carrier are modeled. And the conclusions were applied in the migration studies in the Jungaer Basin, lying on the Northwest part of the China. To construct a reliable migration model, the author studied the characteristics of the sedimentation, the pore fluid pressure evolution, as well as the distribution and the evolution of fluid potential, following the tectonic evolution of the Jungaer Basin. The geochemical prospecting results were used to evidence and to calibrate the migration processes: the oil-source correlation, the distribution of the properties of oil, gas and water. Finally, two important petroleum systems, Permian one and Jurassic one were studied and identified, according, principally, to the studies on the petroleum migration within the Jungaer Basin. Since the oil, as well as the gas, moves mainly in separate phase during the secondary migration, their migrating behaviors would be determined by the dynamics conditions of migration, including the driving forces and pathways. Based on such a consideration, the further understandings may be acquired: the roles played by permeable carriers and low-permeable source rock would be very different in compaction, overpressure generation, petroleum migration, and so on. With the numerical method, the effect of the compaction on the secondary migration was analyzed and the results show that the pressure gradient and the flux resulted from compaction are so small that could be neglected by comparing to the buoyancy of oil. The main secondary migration driving forces are therefore buoyancy and capillary within a hydrostatic system. Modeling with the commercial software-Basin_Mod, the migration pathways of petroleum in clastic carriers seem to be inhomogeneous, controlled by heterogeneity of the driving force, which in turn resulted from the topography of seals, the fabrics and the capillary pressure of the clastic carriers. Furthermore, the direct and indirect methods to study fault-sealing properties in the course of migration were systemically summarized. They may be characterized directly by lithological juxtaposition, clay smear and diagenesis, and indirectly the comparing the pressures and fluid properties in the walls at two apartments of a fault. In Jungaer Basin, the abnormal pressures are found in the formations beneath Badaowan or Baijantan Formation. The occurrence of the overpressure seems controlled by the stratigraphy. The rapid sedimentation, tectonic pressuring, clay sealing, chemical diagensis were considered as the principal pressuring mechanisms. The evolution of fluid pressure is influenced differently at different parts of the basin by the tectonic stresses. So the basin appears different pressure evolution cycles from each part to another during the geological history. By coupling the results of thermal evolution, pressure evolution and organic matter maturation, the area and the period of primary migration were acquired and used to determine the secondary migration time and range. The primary migration in Fengcheng Formation happened from latter Triassic to early Jurassic in the main depressions. The main period of lower-Wuerhe Formation was at latter Jurassic in Changji, Shawan and Pen-1-jing-xi Depression, and at the end of early Cretaceous in Mahu Depression. The primary migration in Badaowan and Sangonghe Formation is at the end of early-Cretaceous in Changji Depression. After then, the fluid potential of oil is calculated at the key time determined from area and time of the primary migration. Generally, fluid potential of oil is high in the depressions and low at the uplifts. Synthetically, it is recognized that the petroleum migration in the Jungaer Basin is very complex, that leads us to classify the evolution of petroleum systems in Northwestern China as a primary stage and a reformed one. The remigration of accumulated petroleum, caused by the reformation of the basin, results in the generation of multiple petroleum systems. The faults and unconformities are usually the linkers among the original petroleum systems. The Permian petroleum system in Jungaer Basin is such a multiple petroleum system. However, the Jurassic petroleum system stays still in its primary stage, thought the strong influences of the new tectonic activities.