Efficient generation of conditional simulations by chebyshev matrix polynomial approximations to the symmetric square root of the covariance matrix     

C. R. Dietrich  G. N. Newsam 《Mathematical Geology》1995,27(2):207-228

Consider the problem of generating a realization y₁ of a Gaussian random field on a dense grid of points ₁ conditioned on field observations y₂ collected on a sparse grid of points ₂. An approach to this is to generate first an unconditional realization y over the grid =_{1 2}, and then to produce y₁ by conditioning y on the data y₂. As standard methods for generating y, such as the turning bands, spectral or Cholesky approaches can have various limitations, it has been proposed by M. W. Davis to generate realizations from a matrix polynomial approximations to the square root of the covariance matrix. In this paper we describe how to generate a direct approximation to the conditional realization y₁, on ₁ using a variant of Davis' approach based on approximation by Chebyshev polynomials. The resulting algorithm is simple to implement, numerically stable, and bounds on the approximation error are readily available. Furthermore we show that the conditional realization y₁ can be generated directly with a lower order polynomial than the unconditional realization y, and that further reductions can be achieved by exploiting a nugget effect if one is present. A pseudocode version of the algorithm is provided that can be implemented using the fast Fourier transform if the field is stationary and the grid ₁ is rectangular. Finally, numerical illustrations are given of the algorithm's performance in generating various 2-D realizations of conditional processes on large sampling grids.  相似文献   

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

Akhil Datta-Gupta  Larry W. Lake  Gary A. Pope 《Mathematical Geology》1995,27(6):763-787

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.  相似文献   

Regional Seas and Sustainable Development Results of an International Workshop     

Josep Maria Bayona  Joan Albaigés 《GeoJournal》1995,36(4):416-416

Reports

Regional Seas and Sustainable Development Results of an International Workshop  相似文献   

Estimation for partially observed Markov processes     

M. E. Thompson  T. N. Kaseke 《Stochastic Hydrology and Hydraulics》1995,9(1):33-47

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.  相似文献   

合肥遥测台网地震速报系统     

刘东旺  黎捷 《地震地磁观测与研究》1995,16(6):27-34

论述了适用于遥测地震台网大震速报的一套计算机处理系统,系统设计以《全国遥测台网观测技术规范》和《近震分析》理论为依据,实现了数据输入,震相选择,数据处理,地震定位,震级计算,结果输出及贮存等多项功能,达到了快速确定地震三要素,提高速报速度和精度之目的,是一个较完善的,实用的速报系统。  相似文献   

Sedimentary Features of Shallow Ancient River Channels on the Northern Shelf of the South China Sea   总被引：1，自引：0，他引：1  

Kou Yangqi  Du Deli nd Marine Geological Investigation Party  MGMR  Guangzhou  Guangdong Jiang M inxi 《《地质学报》英文版》1995,69(1)

Quaternary buried ancient river channels are widespread in the shallow-level sediments of the northern shelf of the South China Sea. The sedimentary sequence mainly of fluvial deposits comprise an important component part of the low-stand system tract and transgressive system tract in the study region. The plannar variation and spatial association of the sedimentary features such as incised valley fillings, deltaic foreset wedges and block slides of shelf-marginal fans reflect the palaeogeographic environment during the fall of the regional sea level in the northern part of the South China Sea. Based on the high-resolution seismic reflection data and gelogical data from boreholes, the present paper makes an integrated interpretation of the Quaternary ancient river channels in the shallow sediments of the study area, studies the sedimentary features of the ancient channels such as their spatial distribution, seismic facies reflection indicators, sedimentary facies and sand -body types, and discusses thei  相似文献   

湘中“桃江式”锰矿含锰岩系的定量研究及其找矿意义          下载免费PDF全文

潘汉军  侯景儒  张廷勋 《地质找矿论丛》1995,10(2):16-27

本文从分析地质体数学特征的思想出发,来研究“桃江式”锰矿的含矿母岩－“含锰岩系”的数学特征。含锰岩系主要由五种岩石组成：粘土岩、黑色页岩、条带状页岩、含锰灰岩和碳酸锰矿层,其数学特征研究包括其统计特征分析、结构特征分析和空间特征分析。文最后讨论了这些数学特征在寻找同类锰矿中的意义。  相似文献   

扬子东南大陆边缘晚前寒武纪古海洋演化的稀土元素记录   总被引：10，自引：1，他引：9  

伊海生  彭军 《沉积学报》1995,13(4):131-137

本文采用中子活化方法分析了晚震旦～早寒武纪１９件样品的稀土元素组成,发现硅质岩石ＲＥＥ配分形式有着复杂的变化,初步可以划分为五种类型,包括“台地型”、“红海型”、“海水型”及两种特殊的配分形式。剖面上Ｃｅ异常由小到大的演变反映了古海洋底层水由缺氧环境向氧化条件的逐渐转化。正Ｅｕ异常的出现代表了热水沉积事件的存在。  相似文献   

    

Daniel Y. Gezari  Dana E. Backman 《Astrophysics and Space Science》1995,224(1-2):463

Mid-infrared imaging photometry of the Orion BN/KL infrared cluster at eight wavelengths between 5 and 20µm using a 58 × 62 pixel imaging array camera has revealed new compact sources and the large-scale structure of the region in diffraction-limited (1 arcsec) detail. Several new objects have been detected within a few arcsec of IRc2, widely thought to be the principal luminosity source for the entire BN/KL complex. Detailed color temperature and emission opacity images are derived from the 7.8, 12.4 and 20.0µm observations, and the 9.8µm image is used to derive an image of silicate dust extinction for the region. The color temperature, opacity, and extinction images show that IRc2 may not be the single dominant luminosity source for the BN/KL region; substantial contributions to the luminosity could be made by IRc7, BN, KL, and five new compact 10µm sources detected within a few arcseconds of IRc2. We suggest that a luminous, early-type star near IRc2, which is associated with the compact radio source I and the Orion SiO maser, is the dominant luminosity source in the BN/KL region, hidden from view by cool dust material with at least Av 60 mag of visible extinction.  相似文献   

A topological approach to understand a multiple-loop solar flare     

L. G. Bagalá  C. H. Mandrini  M. G. Rovira  P. Démoulin  J. C. Hénoux 《Solar physics》1995,161(1):103-121

We analyze the UV and X-ray data obtained by the SMM satellite for the flare starting at 02:36 UT on November 12, 1980 in AR 2779. From a detailed revision of the Ov emission, we find that the observations are compatible with energy being released in a zone above the magnetic inversion line of the AR intermediate bipole. This energy is then transported mainly by conduction towards the two distant kernels located in the AR main bipole. One of these kernels is first identified in this paper. Accelerated particles contribute to the energy transport only during the impulsive phase.We model the observed longitudinal magnetic field by means of a discrete number of subphotospheric magnetic poles, and derive the magnetic field overall topology. As in previous studies of chromospheric flares, the Ov kernels are located along the intersection of the computed separatrices with the photosphere. Especially where the field-line linkage changes discontinuously, these kernels can be linked in pairs by lines that extend along separatrices. Our results agree with the hypothesis of magnetic energy released by magnetic reconnection occurring on separatrices.  相似文献