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

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.     

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.     

Correlations in aftershock and seismicity patterns

Correlations in space and time play a fundamental role in earthquake processes. One direct manifestation of the effects of correlations is the occurrence of aftershocks due to the stress transfer in the vicinity of a main shock. Less obvious and more speculative changes in correlations may occur in the background seismicity before large earthquakes. Using statistical physics it is possible to introduce a measure of spatial correlations through a correlation length. This quantity characterizes how local fluctuations can influence the occurrence of earthquakes over distances comparable with the correlation length. In this work, the physical basis of spatial correlations of earthquakes is discussed in the context of critical phenomena and the percolation problem. The method of two-point correlation function is applied to the seismicity of California. Well defined variations in time of the correlation length are found for aftershock sequences and background seismicity. The scaling properties of our obtained distributions are analyzed with respect to changes in several scaling parameters such as lower magnitude cutoff of earthquakes, the maximum time interval between earthquakes, and the spatial size of the area considered. This scaling behavior can be described in a unified manner by utilizing the multifractal fit. Utilizing the percolation approach the time evolution of clusters of earthquakes is studied with the correlation length defined in terms of the radius of gyration of clusters. This method is applied to the seismicity of California.     

A way to synchronize models with seismic faults for earthquake forecasting: Insights from a simple stochastic model

Numerical models are starting to be used for determining the future behaviour of seismic faults and fault networks. Their final goal would be to forecast future large earthquakes. In order to use them for this task, it is necessary to synchronize each model with the current status of the actual fault or fault network it simulates (just as, for example, meteorologists synchronize their models with the atmosphere by incorporating current atmospheric data in them). However, lithospheric dynamics is largely unobservable: important parameters cannot (or can rarely) be measured in Nature. Earthquakes, though, provide indirect but measurable clues of the stress and strain status in the lithosphere, which should be helpful for the synchronization of the models.The rupture area is one of the measurable parameters of earthquakes. Here we explore how it can be used to at least synchronize fault models between themselves and forecast synthetic earthquakes. Our purpose here is to forecast synthetic earthquakes in a simple but stochastic (random) fault model. By imposing the rupture area of the synthetic earthquakes of this model on other models, the latter become partially synchronized with the first one. We use these partially synchronized models to successfully forecast most of the largest earthquakes generated by the first model. This forecasting strategy outperforms others that only take into account the earthquake series. Our results suggest that probably a good way to synchronize more detailed models with real faults is to force them to reproduce the sequence of previous earthquake ruptures on the faults. This hypothesis could be tested in the future with more detailed models and actual seismic data.     

T -distributed Random Fields: A Parametric Model for Heavy-tailedWell-log Data1

Histograms of observations from spatial phenomena are often found to be more heavy-tailed than Gaussian distributions, which makes the Gaussian random field model unsuited. A T-distributed random field model with heavy-tailed marginal probability density functions is defined. The model is a generalization of the familiar Student-T distribution, and it may be given a Bayesian interpretation. The increased variability appears cross-realizations, contrary to in-realizations, since all realizations are Gaussian-like with varying variance between realizations. The T-distributed random field model is analytically tractable and the conditional model is developed, which provides algorithms for conditional simulation and prediction, so-called T-kriging. The model compares favourably with most previously defined random field models. The Gaussian random field model appears as a special, limiting case of the T-distributed random field model. The model is particularly useful whenever multiple, sparsely sampled realizations of the random field are available, and is clearly favourable to the Gaussian model in this case. The properties of the T-distributed random field model is demonstrated on well log observations from the Gullfaks field in the North Sea. The predictions correspond to traditional kriging predictions, while the associated prediction variances are more representative, as they are layer specific and include uncertainty caused by using variance estimates.     

Ecological characteristics of Vrishabhavathy River in Bangalore (India)

River Vrishabhavathy, a tributary of Cauvery River was studied for 12 physico-chemical parameters at four sites over a distance of 50 km for a period of 2 years (1999–2001) at monthly intervals. Water was faintly alkaline, with pH showing negative correlation with temperature. The dissolved oxygen content increased downstream with negative correlation to biological oxygen demand (BOD), chemical oxygen demand (COD) and turbidity. Bicarbonate alkalinity was very low compared with carbonate alkalinity. Carbonate alkalinity, total hardness, total dissolved solids, total suspended solids, electrical conductivity, BOD and COD decreased downstream, with an upward trend in the middle reaches due to the introduction of raw sewage. The seasonal and yearly trends are also discussed. The river is a sewer collector undergoing self-purification.     

VML在WebGIS图形显示中的应用研究

针对目前被动式WebGIS地图显示的局限,提出了引入VML的解决办法,使客户端容易实现矢量地图显示、地物动态变化显示、专题制图显示等,丰富了WebGIS的地图表现方式。     

早涝碱咸综合治理与生态环境良性循环

华北平原东部淡水资源短缺，旱涝碱成灾害限制了农业生产的可持续发展。海河的治理，解决了排洪排涝排咸出路。春季开发利用地下水包括微咸水和半咸水抗旱灌溉。夏季利用伏雨洗盐排咸，增大降雨入渗，减少径流流失，防治渍涝灾害，把降雨转化为地下水资源。秋冬引蓄河水，回灌地下水补源。以土壤与潜水的地层空间作为调节大气降水、土壤水、地下水、地表水的地下水库，以调控地下水埋深在临界动态为指标，最大限度地把时空分布不均的天然降雨转化为可持续利用的水资源。地表水地下水联合运用，促使水资源采补平衡，降雨灌溉淋洗脱盐强于干旱蒸发积盐过程，地下水淡化强于矿化过程。实现旱涝碱咸综合治理，水土资源可持续利用，经济社会可持续发展，生态环境良性循环。     

一种基于核学习的储集层渗透率预测新方法

基于核学习的支持向量机,是一种采用结构风险最小化原则代替传统经验风险最小化原则的新型统计学习方法,具有完备的理论基础。这里提出了核学习技术在储集层非均质特性描述中渗透率参数预测的新用途。在复杂地层中,基于支持向量机的智能和自适应模式识别能力而建立了常规测井多参数信息输入的渗透率预测模型,然后对实际油田储集层渗透率进行了预测。与常规线性回归模型预测结果相对比,所提出的方法更易于使用,很少受不确定因素的影响,并具有较强的信息整合能力以及更高的预测准确性和可信度。