共查询到20条相似文献，搜索用时 437 毫秒
1.
A probabilistic parametrization for geological uncertainty estimation using the ensemble Kalman filter (EnKF) 总被引：2，自引：0，他引：2
In the past years, many applications of history-matching methods in general and ensemble Kalman filter in particular have been proposed, especially in order to estimate fields that provide uncertainty in the stochastic process defined by the dynamical system of hydrocarbon recovery. Such fields can be permeability fields or porosity fields, but can also fields defined by the rock type (facies fields). The estimation of the boundaries of the geologic facies with ensemble Kalman filter (EnKF) was made, in different papers, with the aid of Gaussian random fields, which were truncated using various schemes and introduced in a history-matching process. In this paper, we estimate, in the frame of the EnKF process, the locations of three facies types that occur into a reservoir domain, with the property that each two could have a contact. The geological simulation model is a form of the general truncated plurigaussian method. The difference with other approaches consists in how the truncation scheme is introduced and in the observation operator of the facies types at the well locations. The projection from the continuous space of the Gaussian fields into the discrete space of the facies fields is realized through in an intermediary space (space with probabilities). This space connects the observation operator of the facies types at the well locations with the geological simulation model. We will test the model using a 2D reservoir which is connected with the EnKF method as a data assimilation technique. We will use different geostatistical properties for the Gaussian fields and different levels of the uncertainty introduced in the model parameters and also in the construction of the Gaussian fields. 相似文献
2.
Application of the EnKF and Localization to Automatic History Matching of Facies Distribution and Production Data 总被引：3，自引：3，他引：0
The performance of the ensemble Kalman filter (EnKF) for continuous updating of facies location and boundaries in a reservoir model based on production and facies data for a 3D synthetic problem is presented. The occurrence of the different facies types is treated as a random process and the initial distribution was obtained by truncating a bi-Gaussian random field. Because facies data are highly non-Gaussian, re-parameterization was necessary in order to use the EnKF algorithm for data assimilation; two Gaussian random fields are updated in lieu of the static facies parameters. The problem of history matching applied to facies is difficult due to (1) constraints to facies observations at wells are occasionally violated when productions data are assimilated; (2) excessive reduction of variance seems to be a bigger problem with facies than with Gaussian random permeability and porosity fields; and (3) the relationship between facies variables and data is so highly non-linear that the final facies field does not always honor early production data well. Consequently three issues are investigated in this work. Is it possible to iteratively enforce facies constraints when updates due to production data have caused them to be violated? Can localization of adjustments be used for facies to prevent collapse of the variance during the data-assimilation period? Is a forecast from the final state better than a forecast from time zero using the final parameter fields?To investigate these issues, a 3D reservoir simulation model is coupled with the EnKF technique for data assimilation. One approach to enforcing the facies constraint is continuous iteration on all available data, which may lead to inconsistent model states, incorrect weighting of the production data and incorrect adjustment of the state vector. A sequential EnKF where the dynamic and static data are assimilated sequentially is presented and this approach seems to have solved the highlighted problems above. When the ensemble size is small compared to the number of independent data, the localized adjustment of the state vector is a very important technique that may be used to mitigate loss of rank in the ensemble. Implementing a distance-based localization of the facies adjustment appears to mitigate the problem of variance deficiency in the ensembles by ensuring that sufficient variability in the ensemble is maintained throughout the data assimilation period. Finally, when data are assimilated without localization, the prediction results appear to be independent of the starting point. When localization is applied, it is better to predict from the start using the final parameter field rather than continue from the final state. 相似文献
3.
In this paper we present an extension of the ensemble Kalman filter (EnKF) specifically designed for multimodal systems. EnKF
data assimilation scheme is less accurate when it is used to approximate systems with multimodal distribution such as reservoir
facies models. The algorithm is based on the assumption that both prior and posterior distribution can be approximated by
Gaussian mixture and it is validated by the introduction of the concept of finite ensemble representation. The effectiveness
of the approach is shown with two applications. The first example is based on Lorenz model. In the second example, the proposed
methodology combined with a localization technique is used to update a 2D reservoir facies models. Both applications give
evidence of an improved performance of the proposed method respect to the EnKF. 相似文献
4.
M. Panzeri E. L. Della Rossa L. Dovera M. Riva A. Guadagnini 《Computational Geosciences》2016,20(3):637-653
We present a methodology that allows conditioning the spatial distribution of geological and petrophysical properties of reservoir model realizations on available production data. The approach is fully consistent with modern concepts depicting natural reservoirs as composite media where the distribution of both lithological units (or facies) and associated attributes are modeled as stochastic processes of space. We represent the uncertain spatial distribution of the facies through a Markov mesh (MM) model, which allows describing complex and detailed facies geometries in a rigorous Bayesian framework. The latter is then embedded within a history matching workflow based on an iterative form of the ensemble Kalman filter (EnKF). We test the proposed methodology by way of a synthetic study characterized by the presence of two distinct facies. We analyze the accuracy and computational efficiency of our algorithm and its ability with respect to the standard EnKF to properly estimate model parameters and assess future reservoir production. We show the feasibility of integrating MM in a data assimilation scheme. Our methodology is conducive to a set of updated model realizations characterized by a realistic spatial distribution of facies and their log permeabilities. Model realizations updated through our proposed algorithm correctly capture the production dynamics. 相似文献
5.
Andreas S. Stordal Randi Valestrand Hans Arnfinn Karlsen Geir N?vdal Hans Julius Skaug 《Computational Geosciences》2012,16(2):467-482
Over the last years, the ensemble Kalman filter (EnKF) has become a very popular tool for history matching petroleum reservoirs.
EnKF is an alternative to more traditional history matching techniques as it is computationally fast and easy to implement.
Instead of seeking one best model estimate, EnKF is a Monte Carlo method that represents the solution with an ensemble of
state vectors. Lately, several ensemble-based methods have been proposed to improve upon the solution produced by EnKF. In
this paper, we compare EnKF with one of the most recently proposed methods, the adaptive Gaussian mixture filter (AGM), on
a 2D synthetic reservoir and the Punq-S3 test case. AGM was introduced to loosen up the requirement of a Gaussian prior distribution
as implicitly formulated in EnKF. By combining ideas from particle filters with EnKF, AGM extends the low-rank kernel particle
Kalman filter. The simulation study shows that while both methods match the historical data well, AGM is better at preserving
the geostatistics of the prior distribution. Further, AGM also produces estimated fields that have a higher empirical correlation
with the reference field than the corresponding fields obtained with EnKF. 相似文献
6.
7.
Traditional ensemble-based history matching method, such as the ensemble Kalman filter and iterative ensemble filters, usually update reservoir parameter fields using numerical grid-based parameterization. Although a parameter constraint term in the objective function for deriving these methods exists, it is difficult to preserve the geological continuity of the parameter field in the updating process of these methods; this is especially the case in the estimation of statistically anisotropic fields (such as a statistically anisotropic Gaussian field and facies field with elongated facies) with uncertainties about the anisotropy direction. In this work, we propose a Karhunen-Loeve expansion-based global parameterization technique that is combined with the ensemble-based history matching method for inverse modeling of statistically anisotropic fields. By using the Karhunen-Loeve expansion, a Gaussian random field can be parameterized by a group of independent Gaussian random variables. For a facies field, we combine the Karhunen-Loeve expansion and the level set technique to perform the parameterization; that is, for each facies, we use a Gaussian random field and a level set algorithm to parameterize it, and the Gaussian random field is further parameterized by the Karhunen-Loeve expansion. We treat the independent Gaussian random variables in the Karhunen-Loeve expansion as the model parameters. When the anisotropy direction of the statistically anisotropic field is uncertain, we also treat it as a model parameter for updating. After model parameterization, we use the ensemble randomized maximum likelihood filter to perform history matching. Because of the nature of the Karhunen-Loeve expansion, the geostatistical characteristics of the parameter field can be preserved in the updating process. Synthetic cases are set up to test the performance of the proposed method. Numerical results show that the proposed method is suitable for estimating statistically anisotropic fields. 相似文献
8.
集合卡曼滤波由于易于使用而被广泛地应用到陆面数据同化研究中,它是建立在模型为线性、误差为正态分布的假设上,而实际土壤湿度方程是高度非线性的,并且当土壤过干或过湿时会发生样本偏斜.为了全面评估它在同化表层土壤湿度观测来反演土壤湿度廓线的性能,特引入不需要上述假设的采样重要性重采样粒子滤波,比较非线性和偏斜性对同化算法的影响.结果显示:不管是小样本还是大样本,集合卡曼滤波都能快速、准确地逼近样本均值,而粒子滤波只有在大样本时才能缓慢地趋近;此外,集合卡曼滤波的粒子边缘概率密度及其偏度和峰度与粒子滤波完全不同,前者粒子虽不完全满足正态分布,但始终为单峰状态,而后者粒子随同化推进经历了单峰到双峰再到单峰的变化. 相似文献
9.
10.
Using a small ensemble size in the ensemble Kalman filter methodology is efficient for updating numerical reservoir models
but can result in poor updates following spurious correlations between observations and model variables. The most common approach
for reducing the effect of spurious correlations on model updates is multiplication of the estimated covariance by a tapering
function that eliminates all correlations beyond a prespecified distance. Distance-dependent tapering is not always appropriate,
however. In this paper, we describe efficient methods for discriminating between the real and the spurious correlations in
the Kalman gain matrix by using the bootstrap method to assess the confidence level of each element from the Kalman gain matrix.
The new method is tested on a small linear problem, and on a water flooding reservoir history matching problem. For the water
flooding example, a small ensemble size of 30 was used to compute the Kalman gain in both the screened EnKF and standard EnKF
methods. The new method resulted in significantly smaller root mean squared errors of the estimated model parameters and greater
variability in the final updated ensemble. 相似文献
11.
In this work, we present an efficient matrix-free ensemble Kalman filter (EnKF) algorithm for the assimilation of large data
sets. The EnKF has increasingly become an essential tool for data assimilation of numerical models. It is an attractive assimilation
method because it can evolve the model covariance matrix for a non-linear model, through the use of an ensemble of model states,
and it is easy to implement for any numerical model. Nevertheless, the computational cost of the EnKF can increase significantly
for cases involving the assimilation of large data sets. As more data become available for assimilation, a potential bottleneck
in most EnKF algorithms involves the operation of the Kalman gain matrix. To reduce the complexity and cost of assimilating
large data sets, a matrix-free EnKF algorithm is proposed. The algorithm uses an efficient matrix-free linear solver, based
on the Sherman–Morrison formulas, to solve the implicit linear system within the Kalman gain matrix and compute the analysis.
Numerical experiments with a two-dimensional shallow water model on the sphere are presented, where results show the matrix-free
implementation outperforming an singular value decomposition-based implementation in computational time. 相似文献
12.
Reza Tavakoli Gergina Pencheva Mary F. Wheeler Benjamin Ganis 《Computational Geosciences》2013,17(1):83-97
We present a parallel framework for history matching and uncertainty characterization based on the Kalman filter update equation for the application of reservoir simulation. The main advantages of ensemble-based data assimilation methods are that they can handle large-scale numerical models with a high degree of nonlinearity and large amount of data, making them perfectly suited for coupling with a reservoir simulator. However, the sequential implementation is computationally expensive as the methods require relatively high number of reservoir simulation runs. Therefore, the main focus of this work is to develop a parallel data assimilation framework with minimum changes into the reservoir simulator source code. In this framework, multiple concurrent realizations are computed on several partitions of a parallel machine. These realizations are further subdivided among different processors, and communication is performed at data assimilation times. Although this parallel framework is general and can be used for different ensemble techniques, we discuss the methodology and compare results of two algorithms, the ensemble Kalman filter (EnKF) and the ensemble smoother (ES). Computational results show that the absolute runtime is greatly reduced using a parallel implementation versus a serial one. In particular, a parallel efficiency of about 35 % is obtained for the EnKF, and an efficiency of more than 50 % is obtained for the ES. 相似文献
13.
State estimation in high dimensional systems remains a challenging part of real time analysis. The ensemble Kalman filter addresses this challenge by using Gaussian approximations constructed from a number of samples. This method has been a large success in many applications. Unfortunately, for some cases, Gaussian approximations are no longer valid, and the filter does not work so well. In this paper, we use the idea of the ensemble Kalman filter together with the more theoretically valid particle filter. We outline a Gaussian mixture approach based on shrinking the predicted samples to overcome sample degeneracy, while maintaining non-Gaussian nature. A tuning parameter determines the degree of shrinkage. The computational cost is similar to the ensemble Kalman filter. We compare several filtering methods on three different cases: a target tracking model, the Lorenz 40 model, and a reservoir simulation example conditional on seismic and electromagnetic data. 相似文献
14.
Andreas S. Stordal Hans A. Karlsen Geir Nævdal Hans J. Skaug Brice Vallès 《Computational Geosciences》2011,15(2):293-305
The nonlinear filtering problem occurs in many scientific areas. Sequential Monte Carlo solutions with the correct asymptotic
behavior such as particle filters exist, but they are computationally too expensive when working with high-dimensional systems.
The ensemble Kalman filter (EnKF) is a more robust method that has shown promising results with a small sample size, but the
samples are not guaranteed to come from the true posterior distribution. By approximating the model error with a Gaussian
distribution, one may represent the posterior distribution as a sum of Gaussian kernels. The resulting Gaussian mixture filter
has the advantage of both a local Kalman type correction and the weighting/resampling step of a particle filter. The Gaussian
mixture approximation relies on a bandwidth parameter which often has to be kept quite large in order to avoid a weight collapse
in high dimensions. As a result, the Kalman correction is too large to capture highly non-Gaussian posterior distributions.
In this paper, we have extended the Gaussian mixture filter (Hoteit et al., Mon Weather Rev 136:317–334, 2008) and also made the connection to particle filters more transparent. In particular, we introduce a tuning parameter for the
importance weights. In the last part of the paper, we have performed a simulation experiment with the Lorenz40 model where
our method has been compared to the EnKF and a full implementation of a particle filter. The results clearly indicate that
the new method has advantages compared to the standard EnKF. 相似文献
15.
Anca Maria Hanea Maria Gheorghe Remus Hanea Dan Ababei 《Computational Geosciences》2013,17(6):929-949
Reservoir simulation models are used both in the development of new fields and in developed fields where production forecasts are needed for investment decisions. When simulating a reservoir, one must account for the physical and chemical processes taking place in the subsurface. Rock and fluid properties are crucial when describing the flow in porous media. In this paper, the authors are concerned with estimating the permeability field of a reservoir. The problem of estimating model parameters such as permeability is often referred to as a history-matching problem in reservoir engineering. Currently, one of the most widely used methodologies which address the history-matching problem is the ensemble Kalman filter (EnKF). EnKF is a Monte Carlo implementation of the Bayesian update problem. Nevertheless, the EnKF methodology has certain limitations that encourage the search for an alternative method.For this reason, a new approach based on graphical models is proposed and studied. In particular, the graphical model chosen for this purpose is a dynamic non-parametric Bayesian network (NPBN). This is the first attempt to approach a history-matching problem in reservoir simulation using a NPBN-based method. A two-phase, two-dimensional flow model was implemented for a synthetic reservoir simulation exercise, and initial results are shown. The methods’ performances are evaluated and compared. This paper features a completely novel approach to history matching and constitutes only the first part (part I) of a more detailed investigation. For these reasons (novelty and incompleteness), many questions are left open and a number of recommendations are formulated, to be investigated in part II of the same paper. 相似文献
16.
The ensemble Kalman filter (EnKF) has become a popular method for history matching production and seismic data in petroleum
reservoir models. However, it is known that EnKF may fail to give acceptable data matches especially for highly nonlinear
problems. In this paper, we introduce a procedure to improve EnKF data matches based on assimilating the same data multiple
times with the covariance matrix of the measurement errors multiplied by the number of data assimilations. We prove the equivalence
between single and multiple data assimilations for the linear-Gaussian case and present computational evidence that multiple
data assimilations can improve EnKF estimates for the nonlinear case. The proposed procedure was tested by assimilating time-lapse
seismic data in two synthetic reservoir problems, and the results show significant improvements compared to the standard EnKF.
In addition, we review the inversion schemes used in the EnKF analysis and present a rescaling procedure to avoid loss of
information during the truncation of small singular values. 相似文献
17.
为研究观测资料稀少情况下土壤质地及有机质对土壤水分同化的影响,发展了集合卡尔曼平滑(Ensemble Kalman Smooth, EnKS)的土壤水分同化方案。利用黑河上游阿柔冻融观测站2008年6月1日至10月29日的观测数据,使用EnKS算法将表层土壤水分观测数据同化到简单生物圈模型(Simple Biosphere Model 2, SiB2)中,分析不同方案对土壤水分估计的影响,并与集合卡尔曼滤波算法(EnKF)的结果进行比较。研究结果表明,土壤质地和有机质对表层土壤水分模拟结果影响最大而对深层的影响相对较小;利用EnKF和EnKS算法同化表层土壤水分观测数据,均能够显著提高表层和根区土壤水分估计的精度,EnKS算法的精度略高于EnKF且所受土壤质地和有机质的影响小于EnKF;当观测数据稀少时,EnKS算法仍然可以得到较高精度的土壤水分估计。 相似文献
18.
This paper proposes an augmented Lagrangian method for production optimization in which the cost function to be maximized
is defined as an augmented Lagrangian function consisting of the net present value (NPV) and all the equality and inequality
constraints except the bound constraints. The bound constraints are dealt with using a trust-region gradient projection method.
The paper also presents a way to eliminate the need to convert the inequality constraints to equality constraints with slack
variables in the augmented Lagrangian function, which greatly reduces the size of the optimization problem when the number
of inequality constraints is large. The proposed method is tested in the context of closed-loop reservoir management benchmark
problem based on the Brugge reservoir setup by TNO. In the test, we used the ensemble Kalman filter (EnKF) with covariance
localization for data assimilation. Production optimization is done on the updated ensemble mean model from EnKF. The production
optimization resulted in a substantial increase in the NPV for the expected reservoir life compared to the base case with
reactive control. 相似文献
19.
Wiktoria Lawniczak Remus Hanea Arnold Heemink Dennis McLaughlin 《Computational Geosciences》2009,13(2):245-254
Reservoir management requires periodic updates of the simulation models using the production data available over time. Traditionally,
validation of reservoir models with production data is done using a history matching process. Uncertainties in the data, as
well as in the model, lead to a nonunique history matching inverse problem. It has been shown that the ensemble Kalman filter
(EnKF) is an adequate method for predicting the dynamics of the reservoir. The EnKF is a sequential Monte-Carlo approach that
uses an ensemble of reservoir models. For realistic, large-scale applications, the ensemble size needs to be kept small due
to computational inefficiency. Consequently, the error space is not well covered (poor cross-correlation matrix approximations)
and the updated parameter field becomes scattered and loses important geological features (for example, the contact between
high- and low-permeability values). The prior geological knowledge present in the initial time is not found anymore in the
final updated parameter. We propose a new approach to overcome some of the EnKF limitations. This paper shows the specifications
and results of the ensemble multiscale filter (EnMSF) for automatic history matching. EnMSF replaces, at each update time,
the prior sample covariance with a multiscale tree. The global dependence is preserved via the parent–child relation in the
tree (nodes at the adjacent scales). After constructing the tree, the Kalman update is performed. The properties of the EnMSF
are presented here with a 2D, two-phase (oil and water) small twin experiment, and the results are compared to the EnKF. The
advantages of using EnMSF are localization in space and scale, adaptability to prior information, and efficiency in case many
measurements are available. These advantages make the EnMSF a practical tool for many data assimilation problems. 相似文献
20.