Geochemistry of magmatic rocks of the Koksharovka alkaline-ultrabasic massif,primorye, and results of trace-element modeling

New geochemical data are presented on the magmatic rocks of the Late Jurassic Koksharovka alkaline-ultrabasic massif, which is associated with deposits of vermiculite, apatite, V-bearing titanomagnetite, and placer isoferroplatinum. The REE geochemistry and strontium, oxygen, and carbon isotopic composition of carbonatites and related ijolites and pyroxenites, together with geological observations, point to the magmatic origin of the Koksharovka carbonatites. The origin of associated magmatic rocks is discussed. Trace element modeling of partial melting of mantle sources was conducted to decipher the genesis of the melts of the Koksharovka carbonatites and host titanite-kaersutite pyroxenites.     

Changing grading of soil: effect on critical states

Examples of situations are presented where the grading of a soil changes during its lifetime either by crushing of particles leading to an increase of fine material or by slow transport of fine particles with seepage leading to a decrease of fine material. Such grading changes influence the basic constitutive properties of the soil, in particular properties such as critical states which are dependent on the available range of densities of packing. Discrete element modelling is used to show the dependence of critical state conditions on grading and the way in which the particle assembly seeks out new critical state conditions as the grading changes.     

Oil and gas resource assessment: The linkage between geological analysis and discovery process model results

Two widely-used techniques to estimate the volume of remaining oil and gas resources are discovery process modeling and geologic assessment. Both were used in a recent national assessment of oil and gas resources of the United States. Parallel estimates were obtained for 27 provinces. Geological-based estimates can typically see into areas not available to discovery process models (that is areas with little or no exploration history) and thus, on average, yield higher estimates. However, a linear relation does exist between the mean estimates obtained from these two methods. In addition, other variables were found in a multiple regression model that explained much of the difference. Thus, it is possible to perform discovery process modeling and adjust the estimates to yield results that might be expected from geological-based assessments.     

The propagation of a single shear crack: A displacement discontinuity model

Numerical studies using the displacement discontinuity method show that a single shear crack under compression propagates in its own direction, because such propagation results in the maximum release of strain energy. The methods of linear elastic fracture mechanics may not be used for such a closed crack, and the stress intensity factors are meaningless in that case. Laboratory observations of propagation by means of kinks at an angle of approximately 70° to the crack may be due to heterogeneities, to the effect of a preexisting crack, to end effects, to microcracking, or to some combination of these factors. Such kinks may thus be local phenomena which cannot release most of the strain energy, and are not incompatible with our numerical results which are based on a global energy balance.     

Interpolation by regularized spline with tension: I. Theory and implementation

Bivariate and trivariate functions for interpolation from scattered data are derived. They are constructed by explicit minimization of a general smoothness functional, and they include a tension parameter that controls the character of the interpolation function (e.g., for bivariate case the surface can be tuned from a membrane to a thin steel plate), Tension can be applied also in a chosen direction, for modeling of phenomena with a simple type of anisotropy. The functions have regular derivatives of all orders everywhere. This makes them suitable for analysis of surface geometry and for direct application in models where derivatives are necessary. For processing of large datasets (thousands of data points), which are now common in geosciences, a segmentation algorithm with a flexible size of overlapping neighborhood is presented. Simple examples demonstrating flexibility and accuracy of the functions are presented.On leave from the Department of Physical Geography and Cartography, Comenius University, Mlynská dolina, Bratislava, Czechoslovakia.On leave from the Institute of Physics, Dúbravská cesta 9, Bratislava, Czechoslovakia.     

内蒙古中西部宝音图群石榴角闪岩的地球化学、变质作用及年代学特征

宝音图群分布于内蒙古中西部狼山、图古日格和达茂旗一带，主要由云母片岩、石英片岩、石英岩、角闪岩和大理岩组成。在狼山格尔敖包沟和图古日格西南地区发育十分典型的石榴角闪岩，均以似层状或透镜状产自云母片岩中，但两地具有不同的矿物组合，前者主要以石榴石+角闪石+斜长石+含钛矿物（金红石和钛铁矿）为特征，后者主要以石榴石+角闪石+绿帘石+含钛矿物（金红石和钛铁矿）为特征。地球化学研究显示两个地区的石榴角闪岩原岩具有相似的化学组成，均为亚碱性玄武岩的拉斑系列。宝音图群的石榴角闪岩具有右倾的稀土配分模式，REE总量为83.31×10^-6~125.9×10^-6，（La/Yb）_N比值为2.17~6.48，δEu=0.87~0.98。Ta、Nb、Ti没有明显的负异常，其配分型式类似E-MORB特点。微量元素构造判别图解表明这些石榴角闪岩的原岩产于板内拉张环境中。本文通过模拟格尔敖包沟样品LS01和图古日格西南地区样品LS35的P-T视剖面图研究其变质作用及矿物演化过程，结果显示这两个样品均经历了早期进变质，峰期以及峰后近等温减压的顺时针型P-T轨迹。样品LS01的峰期温压条件为~11kbar/~735℃，模拟得到的峰期矿物组合可能为石榴石+角闪石+斜长石+黑云母+透辉石+金红石+熔体，P-M（H₂O）视剖面图显示，由于峰后演化过程中存在饱和流体渗透，导致峰期透辉石无法保留。样品LS35的峰期温度压力条件为~8kbar/675℃，其峰期组合为石榴石+角闪石+绿帘石+金红石。两个地区的石榴角闪岩的温度压力条件明显不同，推测宝音图群中呈透镜体或似层状产出的变质基性岩，存在着递增型的变质作用。锆石LA-ICP-MS U-Pb年代学的研究结果显示石榴角闪岩样品LS01的变质锆石年龄为394±8Ma。由此推测这些石榴角闪岩的原岩可能形成于华北北缘在新元古代发育的裂陷盆地，在泥盆纪中期，宝音图群作为华北克拉通的一部分，被卷入到与古亚洲洋闭合有关的造山过程，发育中压型递增变质作用。     

Characterizing heterogeneous permeable media with spatial statistics and tracer data using sequential simulated annealing

Characterizing heterogeneous permeable media using flow and transport data typically requires solution of an inverse problem. Such inverse problems are intensive computationally and may involve iterative procedures requiring many forward simulations of the flow and transport problem. Previous attempts have been limited mostly to flow data such as pressure transient (interference) tests using multiple observation wells. This paper discusses an approach to generating stochastic permeability fields conditioned to geologic data in the form of a vertical variogram derived from cores and logs as well as fluid flow and transport data, such as tracer concentration history, by sequential application of simulated annealing (SA). Thus, the method incorporates elements of geostatistics within the framework of inverse modeling. For tracer-transport calculations, we have used a semianalytic transit-time algorithm which is fast, accurate, and free of numerical dispersion. For steady velocity fields, we introduce a transit-time function which demonstrates the relative importance of data from different sources. The approach is illustrated by application to a set of spatial permeability measurements and tracer data from an experiment in the Antolini Sandstone, an eolian outcrop from northern Arizona. The results clearly reveal the importance of tracer data in reproducing the correlated features (channels) of the permeability field and the scale effects of heterogeneity.     

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.     

弹性波的三维有限元模拟

本文讨论了弹性波瞬态传播问题的三维有限元计算方法及当前存在的实际困难.针对要求计算机内存大和计算时间长的问题,采取了改进措施.由于采用了集中质量矩阵和修正的中心差分时间积分显格式相结合的方法,可以使计算机内存和计算时间大为减少;由于采用结点定位法,最适合用于目前发展的并行计算机系统,可使计算速度大大增快;还采用了有效激发震源法,有效激发区是随时间步进的增加而逐步增大,这不仅能节省计算时间,而且使波场的传播过程一目了然,本文计算了由两种介质组成的三维楔形问题,得到若干典型剖面的瞬时波场图及随时间变化的合成地震图.