Browsing by Subject "covariance"
Now showing 1 - 2 of 2
Results Per Page
Sort Options
Item An analysis of Texas rainfall data and asymptotic properties of space-time covariance estimators(2009-06-02) Li, BoThis dissertation includes two parts. Part 1 develops a geostatistical method to calibrate Texas NexRad rainfall estimates using rain gauge measurements. Part 2 explores the asymptotic joint distribution of sample space-time covariance estimators. The following two paragraphs briefly summarize these two parts, respectively. Rainfall is one of the most important hydrologic model inputs and is considered a random process in time and space. Rain gauges generally provide good quality data; however, they are usually too sparse to capture the spatial variability. Radar estimates provide a better spatial representation of rainfall patterns, but they are subject to substantial biases. Our calibration of radar estimates, using gauge data, takes season, rainfall type and rainfall amount into account, and is accomplished via a combination of threshold estimation, bias reduction, regression techniques and geostatistical procedures. We explore a varying-coefficient model to adapt to the temporal variability of rainfall. The methods are illustrated using Texas rainfall data in 2003, which includes WAR-88D radar-reflectivity data and the corresponding rain gauge measurements. Simulation experiments are carried out to evaluate the accuracy of our methodology. The superiority of the proposed method lies in estimating total rainfall as well as point rainfall amount. We study the asymptotic joint distribution of sample space-time covariance esti-mators of stationary random fields. We do this without any marginal or joint distri-butional assumptions other than mild moment and mixing conditions. We consider several situations depending on whether the observations are regularly or irregularly spaced, and whether one part or the whole domain of interest is fixed or increasing. A simulation experiment illustrates the asymptotic joint normality and the asymp- totic covariance matrix of sample space-time covariance estimators as derived. An extension of this part develops a nonparametric test for full symmetry, separability, Taylor's hypothesis and isotropy of space-time covariances.Item Automatic history matching in Bayesian framework for field-scale applications(Texas A&M University, 2006-04-12) Mohamed Ibrahim Daoud, AhmedConditioning geologic models to production data and assessment of uncertainty is generally done in a Bayesian framework. The current Bayesian approach suffers from three major limitations that make it impractical for field-scale applications. These are: first, the CPU time scaling behavior of the Bayesian inverse problem using the modified Gauss-Newton algorithm with full covariance as regularization behaves quadratically with increasing model size; second, the sensitivity calculation using finite difference as the forward model depends upon the number of model parameters or the number of data points; and third, the high CPU time and memory required for covariance matrix calculation. Different attempts were used to alleviate the third limitation by using analytically-derived stencil, but these are limited to the exponential models only. We propose a fast and robust adaptation of the Bayesian formulation for inverse modeling that overcomes many of the current limitations. First, we use a commercial finite difference simulator, ECLIPSE, as a forward model, which is general and can account for complex physical behavior that dominates most field applications. Second, the production data misfit is represented by a single generalized travel time misfit per well, thus effectively reducing the number of data points into one per well and ensuring the matching of the entire production history. Third, we use both the adjoint method and streamline-based sensitivity method for sensitivity calculations. The adjoint method depends on the number of wells integrated, and generally is of an order of magnitude less than the number of data points or the model parameters. The streamline method is more efficient and faster as it requires only one simulation run per iteration regardless of the number of model parameters or the data points. Fourth, for solving the inverse problem, we utilize an iterative sparse matrix solver, LSQR, along with an approximation of the square root of the inverse of the covariance calculated using a numerically-derived stencil, which is broadly applicable to a wide class of covariance models. Our proposed approach is computationally efficient and, more importantly, the CPU time scales linearly with respect to model size. This makes automatic history matching and uncertainty assessment using a Bayesian framework more feasible for large-scale applications. We demonstrate the power and utility of our approach using synthetic cases and a field example. The field example is from Goldsmith San Andres Unit in West Texas, where we matched 20 years of production history and generated multiple realizations using the Randomized Maximum Likelihood method for uncertainty assessment. Both the adjoint method and the streamline-based sensitivity method are used to illustrate the broad applicability of our approach.