Dynamic data integration for structural modeling: model screening approach using a distance-based model parameterization

This paper proposes a novel history-matching method where reservoir structure is inverted from dynamic fluid flow response. The proposed workflow consists of searching for models that match production history from a large set of prior structural model realizations. This prior set represents the reservoir structural uncertainty because of interpretation uncertainty on seismic sections. To make such a search effective, we introduce a parameter space defined with a “similarity distance” for accommodating this large set of realizations. The inverse solutions are found using a stochastic search method. Realistic reservoir examples are presented to prove the applicability of the proposed method.     

Uncertainty as the overlap of alternate conditional distributions

An important task in modern geostatistics is the assessment and quantification of resource and reserve uncertainty. This uncertainty is valuable support information for many management decisions. Uncertainty at specific locations and uncertainty in the global resource is of interest. There are many different methods to build models of uncertainty, including Kriging, Cokriging, and Inverse Distance. Each method leads to different results. A method is proposed to combine local uncertainties predicted by different models to obtain a combined measure of uncertainty that combines good features of each alternative. The new estimator is the overlap of alternate conditional distributions.     

西藏札达盆地控盆断裂有限元数值模拟

札达盆地是喜马拉雅构造带中的一个山间断陷盆地,其演化过程与盆地两侧的控盆构造密切相关。对控盆断裂的构造应力场进行模拟计算,将有助于进一步深化对本区构造控盆的认识。因此,在对盆地构造地质进行详细调查的基础上,结合本区的深部地质与地球物理资料,对札达盆地控盆断裂的构造应力场进行了模拟。计算结果表明,札达盆地的演化明显受盆地两侧边界断裂的控制,札达盆地是在整体南北向挤压应力的作用下,不同块体差异隆升作用的结果。其南侧的控盆断裂为北倾的正断层,北侧的控盆断裂为南倾的逆断层,二者共同形成了南降北升的翘板式断陷盆地运动过程,是喜马拉雅地块在陆内汇聚挤压构造环境中构造应力场调整的一种方式。     

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

系统研究了西藏冈底斯带石炭纪—二叠纪火山岩的时空分布、岩相学、元素及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.     

Influence of dimethylsulfide produced by marine algae on the global climate

DMS (dimethylsulfide), a breakdown product of cellular solutes of many species of macroalgae andphytoplankton plays an important role in regulating global climate and counteracting partly the "greenhouse" effect.In this paper, the advance and prospects of DMS study are reviewed and discussed with respectto DMS sample storage, measurement and importance in regulating global climate and the acidity ofrain and aerosol.     

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

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