应用生烃动力学法研究川东上二叠统烃源岩生烃史

现今川东上二叠统烃源岩生烃史研究具有局限性.首先,没有对上二叠统Ⅰ-Ⅱ1型灰岩生油、Ⅲ型泥岩的生气史分开评价;其次,海相镜质组反射率通过拟合公式换算成镜质体反射率评价烃源岩成熟度具有局限性.针对上述不足,笔者通过热模拟实验,利用化学动力学方法,标定出两类源岩生油、生气的动力学参数,并结合川东地区的埋藏史及热史,模拟出源岩有机质的成烃转化率曲线.研究表明:上二叠统灰岩、泥岩有机质在距今200 Ma和190 Ma分别进入了生油、生气门限,而在距今170 Ma和140 Ma生烃结束.     

西藏冈底斯带石炭纪—二叠纪岛弧造山作用：火山岩及地球化学证据

系统研究了西藏冈底斯带石炭纪—二叠纪火山岩的时空分布、岩相学、元素及Sr、Nd、Pb 同位素地球化学和构造环境、源区性质,并与喜马拉雅带二叠纪火山岩进行了对比研究。冈底斯带石炭纪—二叠纪火山岩近东西向集中分布在冈底斯构造带中北部地带,空间上从东至西火山活动的强度和规模渐次减小,时间上从早至晚火山活动的强度和规模总体由弱到强。冈底斯带石炭纪—二叠纪火山岩形成于活动大陆边缘的岛弧构造环境,从早到晚岛弧造山作用经历了初始岛弧→早期岛弧→成熟岛弧的发展演变过程,火山岩浆来源于富集型地幔部分熔融作用,原始岩浆在形成和演化的过程中有俯冲洋壳及随带的深海沉积物和再循环进人地慢的地壳物质组分的强烈混染,明显不同于受地壳物质组分强烈混染的喜马拉雅带二叠纪陆缘裂陷型火山岩。综合研究冈底斯带及其邻区近年来的最新调查与研究成果,从北向南拟建了石炭纪—二叠纪冈底斯岛弧→雅鲁藏布江弧后裂谷盆地→喜马拉雅陆缘裂陷盆地的弧盆系时空结构演化模式,探讨了冈瓦纳大陆北缘石炭纪—二叠纪活动大陆边缘的岛弧造山作用与青藏高原古特提斯演化的耦合关系及其动力学机制,讨论了冈底斯带松多乡榴辉岩的形成过程。     

幔源CO_2释出机理、脱气模式及成藏机制研究进展

针对幔源CO2如何从地幔岩浆中脱出并进入沉积地层中形成CO2气藏聚集这一关键问题,总结了国内外研究进展和前缘方向。研究表明,地幔深部的碱性玄武岩浆和碱性岩浆才是深部流体和CO2等挥发份大量赋存、渗滤和释出的场所。浅成侵入岩、次火山岩和火山通道等是CO2释放和聚集的有利位置,岩浆期后和岩浆衰弱期的热液活动阶段是CO2大量释放和聚集的有利时期。幔源CO2进入沉积盆地中具有3种脱气模式,即沿岩石圈断裂直接脱气模式、热流底辟体脱气模式和壳内岩浆房-基底断裂组合脱气模式。CO2的固有物化性质决定其运移相态多样,具有运移和聚集过程同步的特征。只有在满足大量的化学消耗及地层水或原油的溶解和耗散之后才能形成CO2有效聚集。幔源CO2成藏和分布主要受岩浆气源体和气源断裂体系的控制。今后,在超临界CO2及其对油气运移聚集的作用、CO2与深大断裂及火山岩的关系、CO2脱气运移机制、CO与常规烃类油气的耦合差异成藏机制等方面仍需要进一步的研究和探索。     

松辽盆地宋芳屯油田扶杨油层孔隙结构及其影响因素

通过镜下观察,宋芳屯油田扶杨油层孔隙类型较多,以粒间孔为主。细管压力曲线上看,曲线较陡,分选差,孔隙结构特征参数总体上较差,与孔隙度、渗透率有一定的关系,大的孔隙对渗透率的贡献最大。扶杨油层的孔隙结构较复杂主要由于成岩作用影响,压实作用使孔隙空间显著缩小,胶结作用和自生粘土矿物生长使孔隙进一步缩小和复杂化,后期的溶蚀作用使孔隙空间增加。     

Multi-Scale Texture Modeling

Patterns of crystallographic preferred orientation are referred to as texture. The specific subject of texture analysis is the experimental determination and interpretation of the statistical distribution of orientations of crystals within a specimen of polycrystalline material, which could be metals or rocks. The objective is to relate an observed pattern of preferred orientation to its generating processes and vice versa. In geosciences, texture of minerals in rocks is used to infer constraints on their tectono-metamorphic history. Since most physical properties of crystals, such as elastic moduli, the coefficients of thermal expansion, or chemical resistance to etching depends on crystal symmetry and orientation, the presence of texture imparts directional properties to the polycrystalline material. A major issue of mathematical texture analysis is the resolution of the inverse problem to determine a reasonable orientation density function on SO(3) from measured pole intensities on , which relates to the inverse of the totally geodesic Radon transform. This communication introduces a wavelet approach into mathematical texture analysis. Wavelets on the two-dimensional sphere and on the rotational group SO(3) are discussed, and an algorithms for a wavelet decomposition on both domains following the ideas of Ta-Hsin Li is given. The relationship of these wavelets on both domains with respect to the totally geodesic Radon transform is investigated. In particular, it is shown that the Radon transform of these wavelets on SO(3) are again wavelets on . A novel algorithm for the inversion of experimental pole intensities to an orientation density function based on this relationship is developed.     

Bayesian estimation of seismic hazard for two sites in Switzerland

Seismic hazard analysis is based on data and models, which both are imprecise and uncertain. Especially the interpretation of historical information into earthquake parameters, e.g. earthquake size and location, yields ambiguous and imprecise data. Models based on probability distributions have been developed in order to quantify and represent these uncertainties. Nevertheless, the majority of the procedures applied in seismic hazard assessment do not take into account these uncertainties, nor do they show the variance of the results. Therefore, a procedure based on Bayesian statistics was developed to estimate return periods for different ground motion intensities (MSK scale).Bayesian techniques provide a mathematical model to estimate the distribution of random variables in presence of uncertainties. The developed method estimates the probability distribution of the number of occurrences in a Poisson process described by the parameter . The input data are the historical occurrences of intensities for a particular site, represented by a discrete probability distribution for each earthquake. The calculation of these historical occurrences requires a careful preparation of all input parameters, i.e. a modelling of their uncertainties. The obtained results show that the variance of the recurrence rate is smaller in regions with higher seismic activity than in less active regions. It can also be demonstrated that long return periods cannot be estimated with confidence, because the time period of observation is too short. This indicates that the long return periods obtained by seismic source methods only reflects the delineated seismic sources and the chosen earthquake size distribution law.     

Some comparisons between mining-induced and laboratory earthquakes

Although laboratory stick-slip friction experiments have long been regarded as analogs to natural crustal earthquakes, the potential use of laboratory results for understanding the earthquake source mechanism has not been fully exploited because of essential difficulties in relating seismographic data to measurements made in the controlled laboratory environment. Mining-induced earthquakes, however, provide a means of calibrating the seismic data in terms of laboratory results because, in contrast to natural earthquakes, the causative forces as well as the hypocentral conditions are known. A comparison of stick-slip friction events in a large granite sample with mining-induced earthquakes in South Africa and Canada indicates both similarities and differences between the two phenomena. The physics of unstable fault slip appears to be largely the same for both types of events. For example, both laboratory and mining-induced earthquakes have very low seismic efficiencies where _a is the apparent stress and is the average stress acting on the fault plane to cause slip; nearly all of the energy released by faulting is consumed in overcoming friction. In more detail, the mining-induced earthquakes differ from the laboratory events in the behavior of as a function of seismic momentM ₀. Whereas for the laboratory events 0.06 independent ofM ₀, depends quite strongly onM ₀ for each set of induced earthquakes, with 0.06 serving, apparently, as an upper bound. It seems most likely that this observed scaling difference is due to variations in slip distribution over the fault plane. In the laboratory, a stick-slip event entails homogeneous slip over a fault of fixed area. For each set of induced earthquakes, the fault area appears to be approximately fixed but the slip is inhomogeneous due presumably to barriers (zones of no slip) distributed over the fault plane; at constant , larger events correspond to larger _a as a consequence of fewer barriers to slip. If the inequality _a / 0.06 has general validity, then measurements of _a =µE _a /M ₀, where is the modulus of rigidity andE _a is the seismically-radiated energy, can be used to infer the absolute level of deviatoric stress at the hypocenter.     

Micromechanics of the velocity and normal stress dependence of rock friction

Among the second-order effects on friction the most important are those of variable normal stress and of slip velocity. Velocity weakening, which is usually considered the source of the stick-slip instability in rock friction, has been observed in velocity stepping experiments with Westerly granite. The friction change, , was –0.01 to –0.008 for a tenfold velocity increase. Using normal closure measurements, we observed dilation upon each increase in sliding rate. We also observed, for the first time, time-dependent closure between surfaces during static loading. The dilation that occurred during the velocity stepping experiment was found to be that expected from the static time-dependent closure phenomenon. This change in closure was used to predict friction change with an elastic contact model. The calculated friction change which results from a change in contact area and asperity interlocking, is in good agreement with the observed velocity dependence of steady-state friction. Variable normal stress during sliding has two effects, first in creating new partial slip contacts and locking some existing fully sliding contacts and second in increasing interlocking, for instance when normal load is suddenly increased. As a result, a transient change in friction occurs upon a sudden change in normal load.     

热液成矿作用的地质地球化学动力学——以胶东牟平-乳山金矿带为例

通过对胶东牟平—乳山热液脉状金,铅锌矿床成矿动力学控制规律的研究,说明了建立热液成矿反应体系的方法,并得出以下认识:（1）区域断裂构造活动,通过影响成矿反应体系的热力学性质和条件,控制元素的富集和分散。（2）热液成矿过程中容矿断裂活动可划分为两种作用方式:脆性破裂和韧—脆性张开,构成热液成矿的两种构造动力学环境。（3）热液成矿反应体系是一种开放的动态的热力学体系。断裂的脆性破裂阶段,使体系处于强烈的过饱和状态,矿质在远离平衡状态下快速结晶,加剧物质分异;在脆—韧性张开阶段,矿质在接近平衡态体系中缓慢晶出。（4）热液演化晚期,金在残余溶液中富集;当断裂再次发生构造脉动震颤时,早期形成的块状黄铁矿矿石碎裂,富金溶液充填其中,形成含裂隙金和晶隙金的富矿石。成矿体系热力学演化与构造动力学条件有利的匹配控制着富矿石的形成。     

铁砚磨穿指蘸墨 不用扬鞭自奋蹄——志《山西地震》创刊20周年

回顾了《山西地震》创刊20年的历史,叙述了《山西地震》.编辑部在自身建设、执行科技期刊编排标准、为防震减灾事业默默奉献的情况,总结了经验、教训,指出了今后努力的方向.