Numerical simulation of a stratigraphic model

A multi-lithology diffusive stratigraphic model is considered, which simulates at large scales in space and time the infill of sedimentary basins governed by the interaction between tectonics displacements, eustatic variations, sediment supply, and sediment transport laws. The model accounts for the mass conservation of each sediment lithology resulting in a mixed parabolic, hyperbolic system of partial differential equations (PDEs) for the lithology concentrations and the sediment thickness. It also takes into account a limit on the rock alteration velocity modeled as a unilaterality constraint. To obtain a robust, fast, and accurate simulation, fully and semi-implicit finite volume discre tization schemes are derived for which the existence of stable solutions is proved. Then, the set of nonlinear equations is solved using a Newton algorithm adapted to the unilaterality constraint, and preconditioning strategies are defined for the solution of the linear system at each Newton iteration. They are based on an algebraic approximate decoupling of the sediment thickness and the concentration variables as well as on a proper preconditioning of each variable. These algorithms are studied and compared in terms of robustness, scalability, and efficiency on two real basin test cases.     

基于ICCP算法的重力辅助惯性导航

迭代最近等值线算法(ICCP)是一种重要的匹配导航算法,文中首先介绍ICCP算法的基本原理,随后在0.2′×0.2′重力异常数据库的基础上,利用ICCP算法进行仿真计算得到最佳匹配位置。最后为了验证匹配位置是否可用于修正惯导误差,提出将匹配位置误差作为观测量,用卡尔曼滤波对惯导系统误差进行最优估计。由最后的仿真结果可以看出,ICCP算法可有效抑制惯导纬度误差的增长,且最大纬度误差不超过2,′以匹配位置误差作为观测量可以用来估计惯导方位误差角。     

顾及DEM误差自相关的坡度计算模型精度分析

基于DEM的坡度计算,其误差来源于DEM误差、DEM结构和坡度计算模型。在顾及DEM误差自相关的前提下,对四种DEM坡度计算模型进行了分析和评价。研究表明,三阶不带权差分能给出较高的坡度计算精度;在局部窗口中,格网点数量越多,坡度计算越准确;等权比不等权的坡度计算模型更准确;DEM误差自相关结构形式对坡度计算无影响。进一步的理论分析和试验分析还表明：DEM误差自相关性的存在,不仅能够改善地形分析的精度,也能改善DEM自身精度。     

基于规则库和网格算法的土地利用现状图自动数字注记

图形自动注记可极大地提高制图的工作效率。针对土地利用现状图数字注记的特点,本文从地图认知角度出发,建立注记自动配置的规则库,提出面状要素自动注记的一种新的算法——网格法。该算法最大的优点是容易实现、效率高,可方便地解决注记压盖冲突等问题。通过开发一个程序实例,对图斑要素进行了自动数字注记,实验表明95％以上的图形注记效果都较为理想。     

Sedimentary characteristics of terrigenous debris at site MD05-2905 in the northeastern part of the South China Sea since 36 ka and evolution of the East Asian monsoon

Sediments with high sedimentation rate at site MD05-2905 in the Northeastern slope of the South China Sea provide unique materials for a high-resolution study on the paleoenvironment. Based on precise dating of AMS ¹⁴C, grain size analysis of terrigenous debris at core MD05-2905 was conducted after organic matter, biological carbonate and biogenic opal were removed. The results show that 15.5–63.5 μm coarse grain size ingredients may indicate East Asian winter monsoon changes and that 2–9 μm fine grain size ingredients may be used as a proxy of evolution of the East Asian summer monsoon. The results of grain size analysis, which suggest East Asian monsoon intensity, reveal that a winter monsoon dominated the glacial regime and a summer monsoon dominated the Holocene regime. It was also shown that the summer monsoon increased gradually, experienced several abrupt changes and reached a culmination in the early Holocene (11200–8500 a B.P.) since 36 ka. Controlled by precession periodicity, it may be related with the amount of solar radiation at the highest stage, which needs further study. __________ Translated from Advances in Earth Science, 2007, 22(10): 1012–1018 [译自: 地球科学进展]     

临港新城地质特征及其对新城规划与建设的影响分析

结合临港新城三维地质调查所取得的成果,分析了该城区地质特征及对新城建设的影响,为新城规划建设提供科学依据和服务。     

Aspherical correction for free vibration of elastically anisotropic solid by means of theXYZ algorithm

By applying the perturbation theory to theXYZ algorithm (a kind of variational method), the difference f in free vibration frequencies between sphere and ellipsoid was approximated as , where i and _i (i = x,y andz) (i=x, y andz) are aspherical coefficients and asphericities of the ellipsoid, respectively. We developed an analytic method to compute the aspherical coefficients719-4 by using theXYZ algorithm. A numerical example was given for an ellipsoidal olivine, and an attempt was made to estimate the asphericities of the specimen by a least-squares method, based on the relationship between frequency shift and asphericity.     

Estimation for partially observed Markov processes

Many stochastic process models for environmental data sets assume a process of relatively simple structure which is in some sense partially observed. That is, there is an underlying process (X_n, n 0) or (X_t, t 0) for which the parameters are of interest and physically meaningful, and an observable process (Y_n, n 0) or (Y_t, t 0) which depends on the X process but not otherwise on those parameters. Examples are wide ranging: the Y process may be the X process with missing observations; the Y process may be the X process observed with a noise component; the X process might constitute a random environment for the Y process, as with hidden Markov models; the Y process might be a lower dimensional function or reduction of the X process. In principle, maximum likelihood estimation for the X process parameters can be carried out by some form of the EM algorithm applied to the Y process data. In the paper we review some current methods for exact and approximate maximum likelihood estimation. We illustrate some of the issues by considering how to estimate the parameters of a stochastic Nash cascade model for runoff. In the case of k reservoirs, the outputs of these reservoirs form a k dimensional vector Markov process, of which only the kth coordinate process is observed, usually at a discrete sample of time points.     

The petroleum generation potential and effective oil window of humic coals related to coal composition and age

A worldwide data set of more than 500 humic coals from the major coal-forming geological periods has been used to analyse the evolution in the remaining (Hydrogen Index, HI) and total (Quality Index, QI) generation potentials with increasing thermal maturity and the ‘effective oil window’ (‘oil expulsion window’). All samples describe HI and QI bands that are broad at low maturities and that gradually narrow with increasing maturity. The oil generation potential is completely exhausted at a vitrinite reflectance of 2.0–2.2%R_o or T_max of 500–510 °C. The initial large variation in the generation potential is related to the original depositional conditions, particularly the degree of marine influence and the formation of hydrogen-enriched vitrinite, as suggested by increased sulphur and hydrogen contents. During initial thermal maturation the HI increases to a maximum value, HI_max. Similarly, QI increases to a maximum value, QI_max. This increase in HI and QI is related to the formation of an additional generation potential in the coal structure. The decline in QI with further maturation is indicating onset of initial oil expulsion, which precedes efficient expulsion. Liquid petroleum generation from humic coals is thus a complex, three-phase process: (i) onset of petroleum generation, (ii) petroleum build-up in the coal, and (iii) initial oil expulsion followed by efficient oil expulsion (corresponding to the effective oil window). Efficient oil expulsion is indicated by a decline in the Bitumen Index (BI) when plotted against vitrinite reflectance or T_max. This means that in humic coals the vitrinite reflectance or T_max values at which onset of petroleum generation occurs cannot be used to establish the start of the effective oil window. The start of the effective oil window occurs within the vitrinite reflectance range 0.85–1.05%R_o or T_max range 440–455 °C and the oil window extends to 1.5–2.0%R_o or 470–510 °C. For general use, an effective oil window is proposed to occur from 0.85 to 1.7%R_o or from 440 to 490 °C. Specific ranges for HI_max and the effective oil window can be defined for Cenozoic, Jurassic, Permian, and Carboniferous coals. Cenozoic coals reach the highest HI_max values (220–370 mg HC/g TOC), and for the most oil-prone Cenozoic coals the effective oil window may possibly range from 0.65 to 2.0%R_o or 430 to 510 °C. In contrast, the most oil-prone Jurassic, Permian and Carboniferous coals reach the expulsion threshold at a vitrinite reflectance of 0.85–0.9%R_o or T_max of 440–445 °C.