UKF算法及其在GPS卫星轨道短期预报中的应用

分析了GPS卫星预报星历,在比较分析EKF和UKF优缺点的基础上,将UKF引入GPS卫星轨道预报研究中.数值模拟和结果分析表明,UKF方法预报更稳定,能有效地提高轨道预报精度和稳定性.     

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

马槽滩磷矿区矿体形态复杂,底板起伏不定,多年来勘查工程间距不断改变,因此用多种方法讨论了合理工程间距.探采对比和数理统计学方法求得控制资源量的合理工程间距分别为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.     

Pelagic records from the Equatorial Ninetyeast Ridge and significant environmental events during the past 3.5 Ma

This paper presents pelagic records of planktic foraminifera, as well as data of stable isotope stratigraphy and carbonate stratigraphy since 3.5 Ma B.P. from site ODP758 in the Ninetyeast Ridge of the Indian Ocean. Based on these data, manifestations and related mechanisms of major tectonic and environmental events such as the rapid uplift of the Himalaya Mountains, “middle Pleistocene climatic transition” and “mid-Brunhes dissolution event” in the region are discussed. According to the analysis and comparison of various indices and changes in terms of foraminifera assemblage, paleotemperature, paleosalinity and themocline from site ODP758, the authors deduce that the paleoclimatic changes might correlate with the mid-Pleistocene transition at 1.4–1.7 Ma B.P. The changes of CaCO₃, mass accumulation rates (MAR) of CaCO₃ and non- CaCO₃ MAR indicate that the loaded terrigenous sediments increased at 1.7 Ma, which is in agreement with the uplift history of the Qinghai-Tibet plateau as shown by the available data. The last two changes coincide with the uplift of the Qinghai-Tibet plateau, hence they are called “Qinghai-Tibet movement” (1.7 Ma), and the “Kunlun-Yellow River movement” (1.2–0.6 Ma). The changes of the CaCO₃ content, coarse fraction (> 150 μm) content and planktonic foraminifera biostratigraphy show that strong dissolution of abyssal CaCO₃ occurred in the study region during 0.5–0.4 Ma. The event was consistent with the “mid-Brunhes dissolution event” in the sedimentary records of the Atlantic Ocean, Pacific Ocean, Indian Ocean and Nansha sea area of the South China Sea. __________ Translated from Geological Bulletin of China, 2007, 26(12): 1627–1632 [译自: 地质通报]     

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.     

Dynamics of the mid-ocean ridge plate boundary: Recent observations and theory

The global mid-ocean ridge system is one of the most active plate boundaries on the earth and understanding the dynamic processes at this plate boundary is one of the most important problems in geodynamics. In this paper I present recent results of several aspects of mid-ocean ridge studies concerning the dynamics of oceanic lithosphere at these diverging plate boundaries. I show that the observed rift valley to no-rift valley transition (globally due to the increase of spreading rate or locally due to the crustal thickness variations and/or thermal anomalies) can be explained by the strong temperature dependence of the power law rheology of the oceanic lithosphere, and most importantly, by the difference in the rheological behavior of the oceanic crust from the underlying mantle. The effect of this weaker lower crust on ridge dynamics is mainly influenced by spreading rate and crustal thickness variations. The accumulated strain pattern from a recently developed lens model, based on recent seismic observations, was proposed as an appealing mechanism for the observed gabbro layering sequence in the Oman Ophiolite. It is now known that the mid-ocean ridges at all spreading rates are offset into individual spreading segments by both transform and nontransform discontinuities. The tectonics of ridge segmentation are also spreading-rate dependent: the slow-spreading Mid-Atlantic Ridge is characterized by distinct bulls-eye shaped gravity lows, suggesting large along-axis variations in melt production and crustal thickness, whereas the fast-spreading East-Pacific Rise is associated with much smaller along-axis variations. These spreading-rate dependent changes have been attributed to a fundamental differences in ridge segmentation mechanisms and mantle upwelling at mid-ocean ridges: the mantle upwelling may be intrinsically plume-like (3-D) beneath a slow-spreading ridge but more sheet-like (2-D) beneath a fast-spreading ridge.     

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

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

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.     

极性转换期间地球磁场形态学研究

本文对采自中国黄土高原西峰（35.7°N,107.6°E）和段家坡（34.2°N,109.2°E）两个剖面中黄土层L8和古土壤层S8的1281块定向古地磁样品做了详细的岩石磁学和古地磁学研究.证实了Matuyama-Brunhes（M-B）极性转换带位于L8的中下部.提出了下列观点:1.M-B极性转换过程与地球磁场方向变化相联系的持续时间为3600-4500a,而与地球磁场强度变化相联系的持续时间则为8000-9000a,即强度变化存在“超前和滞后” 效应;2.M-B转换场的形态是由三次快速倒转和一次不成功的倒转构成,或者说,转换场具有快速变换极性的振荡特征;3.M-B转换过程中地球磁场并不是以轴对称的非偶极子场为主,而是偶极子场至少与非偶极子场相当;4.中国黄土-古土壤沉积物所含磁性矿物的主要成分是磁铁矿,它是研究极性转换期间地球磁场详细结构的良好物质.