水平井含水率上升影响因素

利用塔里木油田塔中4油田(TZ4)底水油藏相关的地质、流体数据建立数值模型。在所建模型的基础上,应用数值模拟计算的累积产油、产水和产液量回归俞启泰水驱特征曲线,以反映水平井见水特征的参数b。以参数b为研究对象,采用正交试验的方法研究不同因素对b值的影响,筛选影响水平井见水特征的主要因素,认为原油粘度、油层厚度、非均质性及水平井水平段在油藏中的位置是影响水平井含水上升趋势的主要指标。最后建立主要因素与b值的关系式,结合俞启泰曲线b值图版,提出预测水平井见水规律的公式——图版法(F-b法)。应用F-b法对塔里木水平井含水率进行预测,并与其他相关方法及实际生产数据对比,认为F-b法可作为预测塔里木油田水平井含水率、估算可采储量的一种有效方法。     

火山沉积岩及其矿床研究中伽马测井曲线解释的误区

根据酸性凝灰岩和酸性火山事件粘土岩的自然伽马测井曲线,对当前流传的"沉积岩层的放射性强度(或放射性核素的含量)随泥质含量的增加而增高"的概念和用自然伽马值及经验公式求泥质含量提出质凝。酸性凝灰岩的自然伽马曲线有高异常响应,若解释为泥岩显然是误解,故将沉积岩伽马曲线高异常一律解释为泥岩是片面的。各类火山事件粘土岩的伽马值相差悬殊,但其泥质含量几乎相等,用它们的伽马值计算泥质含量误差甚大。最后对铝土矿层的伽马曲线稍加解释,指出核测井应用的远景。     

TM影像几何校正算法的精度比较

同大多数的线阵推扫式影像相比,TM影像作为双向扫描型的影像其几何校正方法有其特殊性和复杂性。针对这个问题,本文简要介绍了一般多项式、空间斜墨卡托投影、有理函数模型等几种几何近似校正算法,并利用广州地区的TM影像进行了各种试验分析和精度比较。结果表明:同其他方法相比,空间斜墨卡托投影是一种较好的算法,校正精度最高,能够达到一个像素左右,而且该方法不受地形起伏的影响,无需地面高程信息。     

抗差最小二乘配置在数字化地图几何纠正中的应用讨论

抗差估计具有较好的抗拒异常观测值及粗差的能力,而最小二乘配置又能较好地处理系统误差,本文结合两者的优点,利用抗差最小二乘配置对数字化地图进行几何纠正,其中对协方差函数采用抗差拟合,得到了较好的结果。实验证明在GIS数据处理的扫描数字化地图几何纠正中,抗差最小二乘配置在抗拒异常值和处理系统误差方面优于单纯的最小二乘估计和单纯的最小二乘配置方法。     

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.     

Using Sequential Self-Calibration Method to Identify Conductivity Distribution: Conditioning on Tracer Test data

An iterative inverse method, the sequential self-calibration method, is developed for mapping spatial distribution of a hydraulic conductivity field by conditioning on nonreactive tracer breakthrough curves. A streamline-based, semi-analytical simulator is adopted to simulate solute transport in a heterogeneous aquifer. The simulation is used as the forward modeling step. In this study, the hydraulic conductivity is assumed to be a deterministic or random variable. Within the framework of the streamline-based simulator, the efficient semi-analytical method is used to calculate sensitivity coefficients of the solute concentration with respect to the hydraulic conductivity variation. The calculated sensitivities account for spatial correlations between the solute concentration and parameters. The performance of the inverse method is assessed by two synthetic tracer tests conducted in an aquifer with a distinct spatial pattern of heterogeneity. The study results indicate that the developed iterative inverse method is able to identify and reproduce the large-scale heterogeneity pattern of the aquifer given appropriate observation wells in these synthetic cases.     

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

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

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.     

地震综合效应场函数及其在地震预报中的应用──（一）地震综合效应场函数及其计算方法

在系统地分析了目前各种测震学地震预报方法科学思路的基础上,认为测震学地震预报方法基本上可以分为两大类。一类是以已经发生的一些地震作为未来可能发生的地震的“因”,即由于已经发生的地震对区域应力场的影响,导致未来发生较强地震。这一类包括的预报方法较多,如空区、条带、ｂ值、地震迁移、相关地震等等及其由此衍生出来的各种方法。另一类是把已经发生的一些地震作为区域应力场增强的“果”,即已经发生的地震是区域应力场增强过程中的一种反映,而未来地震不一定是已经发生的地震所导致的结果。这一类包括“地震窗口”、小震群活动等方法。针对第一类方法,各种预报方法都是力图从地震三要素中提取未来地震的信息,而具体作法又都是利用地震三要素这个多维空间的某个剖面。为了从地震活动诸要素的多维空间提取综合信息,我们对每个地震加入了破裂面方位,构成了地震第四要素,并依据地震４要素建立了地震综合效应场函数。地震综合效应场函数概括了多种测震学地震预报方法的科学思路和预报经验,从而可以形成测震学的综合预报方法。     

Estimation of covariance matrix from geochemical data with observations below detection limits

Multivariate statistical analyses have been extensively applied to geochemical measurements to analyze and aid interpretation of the data. Estimation of the covariance matrix of multivariate observations is the first task in multivariate analysis. However, geochemical data for the rare elements, especially Ag, Au, and platinum-group elements, usually contain observations the below detection limits. In particular, Instrumental Neutron Activation Analysis (INAA) for the rare elements produces multilevel and possibly extremely high detection limits depending on the sample weight. Traditionally, in applying multivariate analysis to such incomplete data, the observations below detection limits are first substituted, for example, each observation below the detection limit is replaced by a certain percentage of that limit, and then the standard statistical computer packages or techniques are used to obtain the analysis of the data. If a number of samples with observations below detection limits is small, or the detection limits are relatively near zero, the results may be reasonable and most geological interpretations or conclusions are probably valid. In this paper, a new method is proposed to estimate the covariance matrix from a dataset containing observations below multilevel detection limits by using the marginal maximum likelihood estimation (MMLE) method. For each pair of variables, sayY andZ whose observations containing below detection limits, the proposed method consists of three steps: (i) for each variable separately obtaining the marginal MLE for the means and the variances, , , , and forY andZ: (ii) defining new variables by and and lettingA=C+D andB=C–D, and obtaining MLE for variances, and forA andB; (iii) estimating the correlation coefficient _YZ by and the covariance _YZ by . The procedure is illustrated by using a precious metal geochemical data set from the Fox River Sill, Manitoba, Canada.