Fast Update of Conditional Simulation Ensembles

Gaussian random field (GRF) conditional simulation is a key ingredient in many spatial statistics problems for computing Monte-Carlo estimators and quantifying uncertainties on non-linear functionals of GRFs conditional on data. Conditional simulations are known to often be computer intensive, especially when appealing to matrix decomposition approaches with a large number of simulation points. This work studies settings where conditioning observations are assimilated batch sequentially, with one point or a batch of points at each stage. Assuming that conditional simulations have been performed at a previous stage, the goal is to take advantage of already available sample paths and by-products to produce updated conditional simulations at mini- mal cost. Explicit formulae are provided, which allow updating an ensemble of sample paths conditioned on n ≥ 0 observations to an ensemble conditioned on n + q observations, for arbitrary q ≥ 1. Compared to direct approaches, the proposed formulae proveto substantially reduce computational complexity. Moreover, these formulae explicitly exhibit how the q new observations are updating the old sample paths. Detailed complexity calculations highlighting the benefits of this approach with respect to state-of-the-art algorithms are provided and are complemented by numerical experiments.

Reliability models for existing structures based on dynamic state estimation and data based asymptotic extreme value analysis

The problem of time variant reliability analysis of existing structures subjected to stationary random dynamic excitations is considered. The study assumes that samples of dynamic response of the structure, under the action of external excitations, have been measured at a set of sparse points on the structure. The utilization of these measurements m in updating reliability models, postulated prior to making any measurements, is considered. This is achieved by using dynamic state estimation methods which combine results from Markov process theory and Bayes' theorem. The uncertainties present in measurements as well as in the postulated model for the structural behaviour are accounted for. The samples of external excitations are taken to emanate from known stochastic models and allowance is made for ability (or lack of it) to measure the applied excitations. The future reliability of the structure is modeled using expected structural response conditioned on all the measurements made. This expected response is shown to have a time varying mean and a random component that can be treated as being weakly stationary. For linear systems, an approximate analytical solution for the problem of reliability model updating is obtained by combining theories of discrete Kalman filter and level crossing statistics. For the case of nonlinear systems, the problem is tackled by combining particle filtering strategies with data based extreme value analysis. In all these studies, the governing stochastic differential equations are discretized using the strong forms of Ito-Taylor's discretization schemes. The possibility of using conditional simulation strategies, when applied external actions are measured, is also considered. The proposed procedures are exemplifiedmby considering the reliability analysis of a few low-dimensional dynamical systems based on synthetically generated measurement data. The performance of the procedures developed is also assessed based on a limited amount of pertinent Monte Carlo simulations. (C) 2010 Elsevier Ltd. All rights reserved.

黄土沟壑区小流域土壤pH值的空间分布及条件模拟

土壤pH值是影响土壤养分有效性和化学物质在土壤中行为的主要因素,研究土壤pH值的空间分布特征对于土壤养分管理和土壤污染预测具有重要意义。该文用地统计学方法研究了环境因素复杂的黄土高原小流域土壤pH值空间分布特征。结果表明,黄土沟壑区小流域土壤pH值具有球形—指数套合模型的空间结构特征,其空间异质性主要来源于流域内土地利用和土壤侵蚀等随机因素。与有机质协同的Kriging法能较好地对土壤pH值进行估值,其估值范围小于实测数据,估值误差来源于复杂的环境因素。序贯高斯条件模拟的土壤pH值范围与实测数据接近,模拟的平均值低于实测数据,模拟误差来源于模拟过程中独特的Kriging算法及高斯特性。

大港油田孔南地区孔一段储层流动单元研究

The main reservoir type in the south of Dagang Oilfield is alluvial reservoir. In this paper, the reservoir structure model and the distribution of connected body and flow barrier were built on base of the study of high-resolution sequential stratigraphic skeleton and fine sedimentary microfacies on level of single sandbody. Utilizing the static and dynamic data synthetically and carrying out the comparision of the classification method for reservoir flow unit in different reservoir, the criterion, which can be used to classify the flow unit in first section of Kongdian formation of Kongnan area, was defined. The qualitative method of well-to-well correlation and the quantitative method of conditional simulation using multiple data are adopted to disclose the oil and water moving regulation in different flow unit and the distribution rule of remaining oil by physical simulation measure. A set of flow unit study method was formed that is suit for the Dagang Oilfield on account of the remaining oil production according to the flow unit. Several outstanding progresses was obtained in the following aspects:It is considered that the reservoir structure of Zao V iow oil group- Zao Vup4 layerand are jigsaw-puzzled reservoir, while ZaoVup3-ZaoVupi layers are labyrinth reservoir,which are studied on base of high-resolution sequential stratigraphic skeleton on the levelof single sandbody in first section of Kongdian formation of Kongnan area and accordingto the study of fine sedimentary microfacies and fault sealeing.When classifying the flow unit, only permeability is the basic parameter using thestatic and dynamic data and, and also different parameters should be chose or deleted, suchas porosity, effective thickness, fluid viscosity and so on, because of the weak or stronginterlayer heterogeneous and the difference of interlayer crude oil character.The method of building predicting-model of flow unit was proposed. This methodis according to the theories of reservoir sedimentology and high-resolution sequencestratigraphic and adopts the quantitative method of well-to well correlation and the quantitative method of stochastic simulation using integrateddense well data. Finally the 3-D predicting-model of flow unit and the interlay er distribution model in flow unit were built which are for alluvial fan and fan delta fades in first section of Kongdian formation of Kongnan area, and nine genetic model of flow unit of alluvial environment that spread in the space were proposed.(4) Difference of reservoir microscopic pore configuration in various flow units and difference of flow capability and oil displacement effect were demonstrated through the physical experiments such as nuclear magnetic resonance (NMR), constant rate mercury penetration, flow simulation and so on. The distribution of remaining oil in this area was predicted combining the dynamic data and numerical modeling based on the flow unit. Remaining oil production measure was brought up by the clue of flow unit during the medium and late course of the oilfield development.

How Porosity and Permeability Vary Spatially With Grain Size, Sorting, Cement Volume, and Mineral Dissolution In Fluvial Triassic Sandstones: The Value of Geostatistics and Local Regression

Although it is well known that sandstone porosity and permeability are controlled by a range of parameters such as grain size and sorting, amount, type, and location of diagenetic cements, extent and type of compaction, and the generation of intergranular and intragranular secondary porosity, it is less constrained how these controlling parameters link up in rock volumes (within and between beds) and how they spatially interact to determine porosity and permeability. To address these unknowns, this study examined Triassic fluvial sandstone outcrops from the UK using field logging, probe permeametry of 200 points, and sampling at 100 points on a gridded rock surface. These field observations were supplemented by laser particle-size analysis, thin-section point-count analysis of primary and diagenetic mineralogy, quantitiative XRD mineral analysis, and SEM/EDAX analysis of all 100 samples. These data were analyzed using global regression, variography, kriging, conditional simulation, and geographically weighted regression to examine the spatial relationships between porosity and permeability and their potential controls. The results of bivariate analysis (global regression) of the entire outcrop dataset indicate only a weak correlation between both permeability porosity and their diagenetic and depositional controls and provide very limited information on the role of primary textural structures such as grain size and sorting. Subdividing the dataset further by bedding unit revealed details of more local controls on porosity and permeability. An alternative geostatistical approach combined with a local modelling technique (geographically weighted regression; GWR) subsequently was used to examine the spatial variability of porosity and permeability and their controls. The use of GWR does not require prior knowledge of divisions between bedding units, but the results from GWR broadly concur with results of regression analysis by bedding unit and provide much greater clarity of how porosity and permeability and their controls vary laterally and vertically. The close relationship between depositional lithofacies in each bed, diagenesis, and permeability, porosity demonstrates that each influences the other, and in turn how understanding of reservoir properties is enhanced by integration of paleoenvironmental reconstruction, stratigraphy, mineralogy, and geostatistics.

Corrected kriging update formulae for batch-sequential data assimilation

Recently, a lot of effort has been spent in the efficient computation of kriging predictors when observations are assimilated sequentially. In particular, kriging update formulae enabling significant computational savings were derived. Taking advantage of the previous kriging mean and variance computations helps avoiding a costly matrix inversion when adding one observation to the TeX already available ones. In addition to traditional update formulae taking into account a single new observation, Emery (2009) proposed formulae for the batch-sequential case, i.e. when TeX new observations are simultaneously assimilated. However, the kriging variance and covariance formulae given in Emery (2009) for the batch-sequential case are not correct. In this paper, we fix this issue and establish correct expressions for updated kriging variances and covariances when assimilating observations in parallel. An application in sequential conditional simulation finally shows that coupling update and residual substitution approaches may enable significant speed-ups.

Análise de incertezas do modelo de teores associado aos investimentos de pesquisa de longo prazo

Os empreendimentos de mineração comumente demandam grande quantidade de investimentos financeiros e, na maioria das vezes, longos períodos de implantação, o que os torna altamente sujeitos a diversas fontes de incertezas. Tais incertezas comumente tendem a diminuir conforme a evolução do projeto. O objetivo deste estudo é correlacionar as incertezas associadas ao modelo de teores de cobre do depósito Sequeirinho com o volume de investimentos realizados ao longo de distintas fases da pesquisa geológica. Este depósito insere-se no contexto do Complexo de Mineração Sossego, localizado no município de Canaã dos Carajás (PA). Primeiramente, foram realizadas 100 simulações para cada domínio litológico em cada campanha de sondagem (pré-1998, 1999, 2000, 2002 e 2003) a partir do método de simulação sequencial gaussiana condicionada aos dados amostrais, totalizando 1.400 possíveis cenários. Para a avaliação das incertezas foram calculados três índices: variância condicional, coeficiente de variação condicional e intervalo interquartil. Por fim, a avaliação dos investimentos foi elaborada a partir dos custos estimados para o desenvolvimento de sondagens e análises químicas. Desde a campanha pré-1998, houve uma tendência de os teores médios do depósito aproximarem-se dos prováveis valores reais observados nas fases finais da pesquisa. No ano de 2000 ocorreu o maior investimento (cerca de 28 milhões de Reais) e a redução das incertezas atingiu o patamar de 15%. Os investimentos desenvolvidos em sondagens posteriores à campanha de 2000 foram da ordem de 9 milhões de Reais (cerca de 12 mil metros de sondagem), porém, não foram constatadas reduções significativas das incertezas. Este investimento seria melhor aproveitado caso fosse redirecionado a novas áreas de prospecção. Além do montante financeiro necessário para a redução das incertezas, foco deste estudo, as variações na interpretação geológica e a locação dos furos de sondagem são variáveis importantes na análise de incertezas associadas aos investimentos em pesquisa geológica.

Aplicação de métodos geoestatísticos no estudo das distribuições espaciais de condutividade hidráulica em áreas do Sistema Aquífero Guarani (SAG) e do Sistema Aquífero Bauru (SAB) no Estado de São Paulo

A condutividade hidráulica (K) é um dos parâmetros controladores da magnitude da velocidade da água subterrânea, e consequentemente, é um dos mais importantes parâmetros que afetam o fluxo subterrâneo e o transporte de solutos, sendo de suma importância o conhecimento da distribuição de K. Esse trabalho visa estimar valores de condutividade hidráulica em duas áreas distintas, uma no Sistema Aquífero Guarani (SAG) e outra no Sistema Aquífero Bauru (SAB) por meio de três técnicas geoestatísticas: krigagem ordinária, cokrigagem e simulação condicional por bandas rotativas. Para aumentar a base de dados de valores de K, há um tratamento estatístico dos dados conhecidos. O método de interpolação matemática (krigagem ordinária) e o estocástico (simulação condicional por bandas rotativas) são aplicados para estimar os valores de K diretamente, enquanto que os métodos de krigagem ordinária combinada com regressão linear e cokrigagem permitem incorporar valores de capacidade específica (Q/s) como variável secundária. Adicionalmente, a cada método geoestatístico foi aplicada a técnica de desagrupamento por célula para comparar a sua capacidade de melhorar a performance dos métodos, o que pode ser avaliado por meio da validação cruzada. Os resultados dessas abordagens geoestatísticas indicam que os métodos de simulação condicional por bandas rotativas com a técnica de desagrupamento e de krigagem ordinária combinada com regressão linear sem a técnica de desagrupamento são os mais adequados para as áreas do SAG (rho=0.55) e do SAB (rho=0.44), respectivamente. O tratamento estatístico e a técnica de desagrupamento usados nesse trabalho revelaram-se úteis ferramentas auxiliares para os métodos geoestatísticos.

Exchanging uncertainty: interoperable geostatistics?

Traditionally, geostatistical algorithms are contained within specialist GIS and spatial statistics software. Such packages are often expensive, with relatively complex user interfaces and steep learning curves, and cannot be easily integrated into more complex process chains. In contrast, Service Oriented Architectures (SOAs) promote interoperability and loose coupling within distributed systems, typically using XML (eXtensible Markup Language) and Web services. Web services provide a mechanism for a user to discover and consume a particular process, often as part of a larger process chain, with minimal knowledge of how it works. Wrapping current geostatistical algorithms with a Web service layer would thus increase their accessibility, but raises several complex issues. This paper discusses a solution to providing interoperable, automatic geostatistical processing through the use of Web services, developed in the INTAMAP project (INTeroperability and Automated MAPping). The project builds upon Open Geospatial Consortium standards for describing observations, typically used within sensor webs, and employs Geography Markup Language (GML) to describe the spatial aspect of the problem domain. Thus the interpolation service is extremely flexible, being able to support a range of observation types, and can cope with issues such as change of support and differing error characteristics of sensors (by utilising descriptions of the observation process provided by SensorML). XML is accepted as the de facto standard for describing Web services, due to its expressive capabilities which allow automatic discovery and consumption by ‘naive’ users. Any XML schema employed must therefore be capable of describing every aspect of a service and its processes. However, no schema currently exists that can define the complex uncertainties and modelling choices that are often present within geostatistical analysis. We show a solution to this problem, developing a family of XML schemata to enable the description of a full range of uncertainty types. These types will range from simple statistics, such as the kriging mean and variances, through to a range of probability distributions and non-parametric models, such as realisations from a conditional simulation. By employing these schemata within a Web Processing Service (WPS) we show a prototype moving towards a truly interoperable geostatistical software architecture.

3D modelling of the transmissivity of granitic rocks surrounding the old Quinta do Bispo mine | Modelação 3D da transmissividade do maciço granítico envolvente à antiga mina da Quinta do Bispo

Apresenta-se uma metodologia para caracterizar a transmissividade dos Granitos Hercínicos e Metasedimentos do Complexo Xisto-Grauváquico do maciço envolvente e subjacente à antiga área mineira de urânio da Quinta do Bispo. Inicia-se com a modelação das litologias e grau de alteração a que se segue a simulação condicional da densidade de fracturação. No final, a densidade de fracturação é convertida num modelo 3D de transmissividade por relação com os resultados dos ensaios de bombagem. The purpose of this work is to present a methodology for characterizing the transmissivity of the Hercynian granites and complex schist–greywacke metasediment rocks surrounding and underlying the old Quinta do Bispo uranium mining site. The methodology encompasses modelling of lithologies and weathering levels, followed by a conditional simulation of fracture density. Fracture density is then converted into a 3D model of transmissivity via a relationship with pumping tests.