DCT域遥感影像融合算法

提出一种基于离散余弦变换(Discrete Cosine Transform,DCT)与IHS变换的多光谱与全色遥感影像融合方法及其改进算法。本文方法依据DCT系数的特点在DCT域进行遥感影像融合,适合压缩格式的遥感影像快速融合。目视效果与客观评价表明,相比常用遥感影像融合方法,本文方法能在提高空间分辨率与保持光谱特性之间得到更好的折衷。     

新鲜煤和氧化煤自燃倾向性的FTIR对比分析

利用傅立叶变换红外光谱（FTIR）分析技术,对乌达矿区低温氧化和新鲜煤样进行了红外光谱测试,通过两类煤样的特征官能团吸光度的对比分析．揭示了煤在氧化前后特征官能团的变化规律：氧化后的煤样中脂肪族和芳烃CH减少,酸酐C=O从无到有,芳烃骨架C=C基本不变,进而研究了其对煤自燃倾向性的影响。结果表明脂肪族和羟基含量越高,煤自燃倾向性越大。     

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.     

中强地震前后能量分布的自仿射特性

根据简化的地震波传播理论,“漫化”地震能量在空间的分配,减弱能量空间分布概率的极端情形及消除统计域为“空”的情形,由此构造类似于二维ｆＢｍ的地震能量空间分布概率曲面。与起伏不平的地形相似,由于基在水平方向及垂直方向的变化不是等比例的,因而这一分形曲面可能是统计自仿射的。     

Weak-contrast approximation of the elastic scattering matrix in anisotropic media

The plane-wave reflection and transmission coefficients at a plane interface between two anisotropic media constitute the elements of the elastic scattering matrix. For a 1-D anisotropic medium the eigenvector decomposition of the system matrix of the transformed elasto-dynamic equations is used to derive a general expression for the scattering matrix. Depending on the normalization of the eigenvectors, the expressions give scattering coefficients for amplitudes or for vertical energy flux.Computing the vertical slownesses and the corresponding polarizations, the eigenvector matrix and its inverse can be found. We give a simple formula for the inverse, regardless of the normalization of the eigenvectors. When the eigenvectors are normalized with respect to amplitudes of displacement (or velocity), the calculation of the scattering matrix for amplitudes is simplified.When the relative changes in all parameters are small, a weak-contrast approximation of the scattering matrix, based on the exactly determined polarization vectors in an average medium, is obtained. The same approximation is also derived directly from the transformed elasto-dynamic equations for a smooth vertically inhomogeneous medium, proving the consistency of the approximation.For monoclinic media, with the mirror symmetry plane parallel to the interface, the approximative scattering matrix is given in terms of analytic expressions for the non-normalized eigenvectors and vertical slownesses. For transversely isotropic media with a vertical axis of symmetry (VTI) and isotropic media, explicit solutions for the weak-contrast approximations of the scattering matrices have been obtained. The scattering matrix for amplitudes for isotropic media is well known. The scattering matrix for vertical energy flux may have applications in AVO analysis and inversion due to the reciprocity of the reflection coefficients for converted waves.Numerical examples for monoclinic and VTI media provide good agreement between the approximative and the exact reflection matrices. It is, however, expected that the approximations cannot be used when the symmetry properties of the two media are very different. This is because the approximation relies on a small relative contrast between the eigenvectors in the two media.Presented at the Workshop Meeting on Seismic Waves in Laterally Inhomogeneous Media, Castle of Trest, Czech Republic, May 22–27, 1995.     

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.     

The seismic hazard assessment of the Dead Sea rift, Jordan

The Dead Sea fault system and its branching faults represent one of the most tectonically active regions in the Middle East. The aim of this study is to highlight the degree of hazards related to the earthquake activities associated with the Dead Sea rift, in terms of speculating the possible future earthquakes. The present investigation mainly is based on available data and vertical crustal modeling of Jordan and the Dead Sea model for the Dead Sea basin with particular emphasis of the recent earthquake activities, which occurred on December 31st, 2003 (Mc = 3.7), February 11th, 2004 (strongest Mc = 4.9 R), and March 15th, 2004 (Mc = 4). The present research examines the location of the strong events and correlates them with the various tectonic elements in the area. The source mechanism of the main shock and the aftershock events is also examined. The analyses were based on the available short period seismogram data, which was recorded at the Natural Resources Authority of Jordan, Seismological Observatory. The seismic energy appears to have migrated from the south to the north during the period from December 31st up to March 12th, where the released seismic energy showed a migration character to the southern block of the eastern side of the Dead Sea, which led the seismic event to occur on March 15th.     

模糊理论在公路边坡稳定性分析中的应用

边坡工程是一个复杂的系统工程。边坡的稳定与否直接关系到其防护工作的程度。因此,如何有效、合理地判断边坡的稳定性对整个工程的经济与社会效益极为重要。影响边坡稳定的因素众多,判断其稳定与否应考虑各因素的综合效果。论文运用模糊理论中的多因素综合评判模型对某高速公路边坡的稳定性进行研究。通过全面考虑影响因素,建立相应的评判因素集及单因素评判矩阵,运用专家评判法和判断矩阵法确定出各因素的权重。最后通过多层次模糊评判得到该边坡的稳定性评判结果。研究证明,该法在边坡稳定性分析评判中是切实可行的。     

Spectral Simulation and Conditioning to Local Continuity Trends in Geologic Modeling

Spectral simulation has gained application in building geologic models due to the advantage of better honoring the spatial continuity of petrophysical properties, such as reservoir porosity and shale volume. Distinct from sequential simulation methods, spectral simulation is a global algorithm in the sense that a global density spectrum is calculated once and the inverse Fourier transform is performed on the Fourier coefficient also only once to generate a simulation realization. The generated realizations honor the spatial continuity structure globally over the whole field instead of only within a search neighborhood, as with sequential simulation algorithms. However, the disadvantage of global spectral simulation is that it traditionally cannot account for the local information such as the local continuity trends, which are often observed in reservoirs and hence are important to be accounted for in geologic models. This disadvantage has limited wider application of spectral simulation in building geologic models. In this paper, we present ways of conditioning geologic models to the relevant local information. To account for the local continuity trends, we first scale different frequency components of the original model with local-amplitude spectrum ratios that are specific to the local trend. The sum of these scaled frequency components renders a new model that displays the desired local continuity trend. The implementation details of this new method are discussed and examples are provided to illustrate the algorithm.     

地震信号的相位恢复及误差分析

常见的几种地震子波有:零相位子波、最小相位子波、最大相位子波和混合相位子波。以不同类型的信号为激发子波,对不同地层模型合成地震道。仅用原始合成地震道记录的振幅谱来恢复原始记录,讨论不同方法下,由于激发子波的不同,对比分析原始记录与恢复记录的差异,及其恢复效果。