基于物方影像匹配的DEM聚类变化检测

提出了DEM的变化检测问题,采用了基于密度的聚类分析方法对匹配所得的初始变化区域进行分析,实现了地表高程变化区域的检测.     

吐哈盆地丘东凝析气藏中侏罗统储层流动单元划分

对吐哈盆地丘东凝析气藏中侏罗统储层的流动单元的划分方法和参数选取进行探讨。针对丘东凝析低渗透气藏的储层和流体特点,纵向上对研究区中侏罗统地层进行小层精细划分,平面上应用SPSS统计分析软件对厚度、孔隙度、泥质含量、流动带指数等4个参数进行聚类分析和判别分析,将研究区中侏罗统储层砂体划分为4类流动单元。结果表明各流动单元的类型与储层物性、沉积微相具有很好的对应关系,能够客观地反映气藏储层的地质特征。虽然丘东凝析低渗透气藏中流动单元的划分和油藏中流动单元的划分没有太大差别,但是参数的选取却与油藏中的存在着明显的不同。     

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

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

马槽滩磷矿区的矿体形态特征及合理勘查工程间距

马槽滩磷矿区矿体形态复杂,底板起伏不定,多年来勘查工程间距不断改变,因此用多种方法讨论了合理工程间距.探采对比和数理统计学方法求得控制资源量的合理工程间距分别为200m×150m和200m×100m;控制半变异函数的变程亦为200m×100m,均与现行规范相吻合。但沿走向和倾向的变异相同,因此建议用正方形网格进行勘查更为合理,间距为150m×150m。     

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

抗差估计具有较好的抗拒异常观测值及粗差的能力,而最小二乘配置又能较好地处理系统误差,本文结合两者的优点,利用抗差最小二乘配置对数字化地图进行几何纠正,其中对协方差函数采用抗差拟合,得到了较好的结果。实验证明在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.     

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.