共查询到20条相似文献，搜索用时 656 毫秒
1.
Relation between two common localisation methods for the EnKF 总被引：3，自引：0，他引：3
This study investigates the relation between two common localisation methods in ensemble Kalman filter (EnKF) systems: covariance
localisation and local analysis. Both methods are popular in large-scale applications with the EnKF. The case of local observations
with non-correlated errors is considered. Both methods are formulated in terms of tapering of ensemble anomalies, which provides
a framework for their comparison. Based on analytical considerations and experimental evidence, we conclude that in practice
the two methods should yield very similar results, so that the choice between them should be based on other criteria, such
as numerical effectiveness and scalability. 相似文献
2.
The ensemble Kalman filter (EnKF) has been shown repeatedly to be an effective method for data assimilation in large-scale
problems, including those in petroleum engineering. Data assimilation for multiphase flow in porous media is particularly
difficult, however, because the relationships between model variables (e.g., permeability and porosity) and observations (e.g.,
water cut and gas–oil ratio) are highly nonlinear. Because of the linear approximation in the update step and the use of a
limited number of realizations in an ensemble, the EnKF has a tendency to systematically underestimate the variance of the
model variables. Various approaches have been suggested to reduce the magnitude of this problem, including the application
of ensemble filter methods that do not require perturbations to the observed data. On the other hand, iterative least-squares
data assimilation methods with perturbations of the observations have been shown to be fairly robust to nonlinearity in the
data relationship. In this paper, we present EnKF with perturbed observations as a square root filter in an enlarged state
space. By imposing second-order-exact sampling of the observation errors and independence constraints to eliminate the cross-covariance
with predicted observation perturbations, we show that it is possible in linear problems to obtain results from EnKF with
observation perturbations that are equivalent to ensemble square-root filter results. Results from a standard EnKF, EnKF with
second-order-exact sampling of measurement errors that satisfy independence constraints (EnKF (SIC)), and an ensemble square-root
filter (ETKF) are compared on various test problems with varying degrees of nonlinearity and dimensions. The first test problem
is a simple one-variable quadratic model in which the nonlinearity of the observation operator is varied over a wide range
by adjusting the magnitude of the coefficient of the quadratic term. The second problem has increased observation and model
dimensions to test the EnKF (SIC) algorithm. The third test problem is a two-dimensional, two-phase reservoir flow problem
in which permeability and porosity of every grid cell (5,000 model parameters) are unknown. The EnKF (SIC) and the mean-preserving
ETKF (SRF) give similar results when applied to linear problems, and both are better than the standard EnKF. Although the
ensemble methods are expected to handle the forecast step well in nonlinear problems, the estimates of the mean and the variance
from the analysis step for all variants of ensemble filters are also surprisingly good, with little difference between ensemble
methods when applied to nonlinear problems. 相似文献
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.
Wavenumber-frequency spectral analysis of different atmospheric variables has been carried out using 25 years of data. The
area considered is the tropical belt 25°S–25°N. A combined FFT-wavelet analysis method has been used for this purpose. Variables
considered are outgoing long-wave radiation (OLR), 850 hPa divergence, zonal and meridional winds at 850, 500 and 200 hPa
levels, sea level pressure and 850 hPa geopotential height. It is shown that the spectra of different variables have some
common properties, but each variable also has few features different from the rest. While Kelvin mode is prominent in OLR
and zonal winds, it is not clearly observed in pressure and geopotential height fields; the latter two have a dominant wavenumber
zero mode not seen in other variables except in meridional wind at 200 hPa and 850 hPa divergences. Different dominant modes
in the tropics show significant variations on sub-seasonal time scales. 相似文献
5.
局域化改进集合卡尔曼滤波（EnKF）可以克服EnKF方法在使用小集合时,对参数识别精度较低的缺陷,其能同化
地下水位观测数据有效识别渗透系数场。实际工作中,溶质运移数据也较容易获得。崔凯鹏（2013）尝试增加溶质运移
数据以改进只同化水流数据对渗透系数的估计结果,但是精度提高有限。本文在其基础上修改模型,进一步增加溶质注
入井,探究同时同化水头和溶质运移数据,对渗透系数场识别效果,之后对比了局域化EnKF与非局域化EnKF参数识别结
果,并分析了溶质影响范围与参数识别的关系。结果表明：同时同化溶质运移和水头资料,比同化单一种类观测数据识别
的渗透系数精度更高;相同实现数目下,局域化EnKF比EnKF对渗透系数场的估计结果与真实场更为接近;仅考虑溶质影
响范围内的渗透系数,同化水头数据在最后时刻参数识别结果好于同化溶质运移数据参数识别结果,但差别不大。 相似文献
6.
The study has been focused on examining the usage and the applicability of ensemble Kalman filtering techniques to the history
matching procedures. The ensemble Kalman filter (EnKF) is often applied nowadays to solving such a problem. Meanwhile, traditional
EnKF requires assumption of the distribution’s normality. Besides, it is based on the linear update of the analysis equations.
These facts may cause problems when filter is used in reservoir applications and result in sampling error. The situation becomes
more problematic if the a priori information on the reservoir structure is poor and initial guess about the, e.g., permeability
field is far from the actual one. The above circumstance explains a reason to perform some further research concerned with
analyzing specific modification of the EnKF-based approach, namely, the iterative EnKF (IEnKF) scheme, which allows restarting
the procedure with a new initial guess that is closer to the actual solution and, hence, requires less improvement by the
algorithm while providing better estimation of the parameters. The paper presents some examples for which the IEnKF algorithm
works better than traditional EnKF. The algorithms are compared while estimating the permeability field in relation to the
two-phase, two-dimensional fluid flow model. 相似文献
7.
Mohammadali Tarrahi Siavash Hakim Elahi Behnam Jafarpour 《Computational Geosciences》2016,20(5):929-952
Ensemble methods present a practical framework for parameter estimation, performance prediction, and uncertainty quantification in subsurface flow and transport modeling. In particular, the ensemble Kalman filter (EnKF) has received significant attention for its promising performance in calibrating heterogeneous subsurface flow models. Since an ensemble of model realizations is used to compute the statistical moments needed to perform the EnKF updates, large ensemble sizes are needed to provide accurate updates and uncertainty assessment. However, for realistic problems that involve large-scale models with computationally demanding flow simulation runs, the EnKF implementation is limited to small-sized ensembles. As a result, spurious numerical correlations can develop and lead to inaccurate EnKF updates, which tend to underestimate or even eliminate the ensemble spread. Ad hoc practical remedies, such as localization, local analysis, and covariance inflation schemes, have been developed and applied to reduce the effect of sampling errors due to small ensemble sizes. In this paper, a fast linear approximate forecast method is proposed as an alternative approach to enable the use of large ensemble sizes in operational settings to obtain more improved sample statistics and EnKF updates. The proposed method first clusters a large number of initial geologic model realizations into a small number of groups. A representative member from each group is used to run a full forward flow simulation. The flow predictions for the remaining realizations in each group are approximated by a linearization around the full simulation results of the representative model (centroid) of the respective cluster. The linearization can be performed using either adjoint-based or ensemble-based gradients. Results from several numerical experiments with two-phase and three-phase flow systems in this paper suggest that the proposed method can be applied to improve the EnKF performance in large-scale problems where the number of full simulation is constrained. 相似文献
8.
岩体结构的变形破坏是一个渐近过程,常规的强度设计与极限分析均无法描述这一过程。在回顾总结岩体结构破坏控制、变形控制的基础上,提出岩体结构变形破坏控制最终要落实到开裂控制上来。总结了岩体结构开裂计算的特点和难点,阐明了用不平衡力描述岩体结构变形破坏过程的理论依据。分别从单轴压缩试验、边坡开挖卸荷、拱坝坝踵坝趾开裂等3个具体的数值案例论证不平衡力与开裂破坏的内在相关性。提出了蓄水导致不可逆的谷幅变形也是不平衡力作用结果:裂隙水压力使屈服面收缩、原先处于屈服或临界屈服状态的岩体应力状态超出屈服面,产生不平衡力和不可逆的塑性变形。指出了不平衡力本质是岩体结构非平衡态到平衡态的距离,岩体结构非平衡演化的总体趋势服从最小塑性余能原理,进一步指出了无法消除的不平衡力是岩体结构变形破坏的内在驱动力。不平衡力不仅可以作为岩体结构变形破坏的判据,还能给出相应的加固措施,具有重要的工程意义。 相似文献
9.
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. 相似文献
10.
11.
The ensemble Kalman filter (EnKF) has been successfully applied to data assimilation in steam-assisted gravity drainage (SAGD) process, but applications of localization for the EnKF in the SAGD process have not been studied. Distance-based localization has been reported to be very efficient for assimilation of large amounts of independent data with a small ensemble in water flooding process, but it is not applicable to the SAGD process, since in the SAGD process, oil is produced mainly from the transition zone steam chamber to cold oil instead of the regions around the producer. As the oil production rate is mainly affected by the temperature distribution in the transition zone, temperature-based localization was proposed for automatic history matching of the SAGD process. The regions of the localization function were determined through sensitivity analysis by using a large ensemble with 1000 members. The sensitivity analysis indicated that the regions of cross-correlations between oil production and state variables are much wider than the correlations between production data and model variables. To choose localization regions that are large enough to include the true regions of non-zero cross-covariance, the localization function is defined based on the regions of non-zero covariances of production data to state variables. The non-zero covariances between production data and state variables are distributed in accordance with the steam chamber. This makes the definition of a universal localization function for different state variables easier. Based on the cross-correlation analysis, the temperature range in which oil production is contributed is determined, and beyond or below this range, the localization function reduces from one, and at the critical temperature or steam temperature, the localization function reduces to zero. The temperature-based localization function was obtained through modifying the distance-based localization function. Localization is applied to covariance of data with permeability, saturation, and temperature, as well as the covariance of data with data. A small ensemble (10 ensemble members) was employed in several case studies. Without localization, the variability in the ensemble collapsed very quickly and lost the ability to assimilate later data. The mean variance of model variables dropped dramatically by 95 %, and there was almost no variability in ensemble forecasts, while the prediction was far from the reference with data mismatch keeping up at a high level. At least 50 ensemble members are needed to keep the qualities of matches and forecasts, which significantly increases the computation time. The EnKF with temperature-based localization is able to avoid the collapse of ensemble variability with a small ensemble (10 members), which saves the computation time and gives better history match and prediction results. 相似文献
12.
青藏高原地区多套位势高度和风场再分析资料的对比分析 总被引：5，自引：1，他引：4
针对20CR、CFSR、NCEP1、NCEP2、ERA-Interim、ERA-20CM和JRA-55再分析位势高度和风场资料,基于探空资料,采用计算均方根误差、相关分析等方法从气候均值、长期变化趋势和年际变率三个方面评估再分析资料在青藏高原地区的适用性.结果表明:再分析资料的适用性存在季节、空间和垂直层次上的差异.气候均值方面,NCEP1和ERA-Interim位势高度资料与观测资料最接近,适用性最佳,而ERA-20CM资料偏差最大.春夏季,NCEP2和NCEP1的风速资料质量较优,20CR和ERA-20CM资料质量相对较差;而秋冬季,ERA-20CM风速资料与探空资料最接近,质量最好.气候变化趋势方面,JRA-55、ERA-Interim和NCEP1资料质量存在时空的差异,但均能很好地反映出位势高度的变化趋势.年际变率方面,除了ERA-20CM,各再分析资料与探空资料相关性高,对年际变率的刻画基本一致,其中JRA-55和ERA-Interim位势高度资料与探空资料相关性最好.就季节而言,冬季再分析资料质量最高,适用性好,其次为春季,夏季资料质量最低. 相似文献
13.
Haiyan Zhou Liangping Li Harrie-Jan Hendricks Franssen J. Jaime Gómez-Hernández 《Mathematical Geosciences》2012,44(2):169-185
The ensemble Kalman filter (EnKF) is now widely used in diverse disciplines to estimate model parameters and update model
states by integrating observed data. The EnKF is known to perform optimally only for multi-Gaussian distributed states and
parameters. A new approach, the normal-score EnKF (NS-EnKF), has been recently proposed to handle complex aquifers with non-Gaussian
distributed parameters. In this work, we aim at investigating the capacity of the NS-EnKF to identify patterns in the spatial
distribution of the model parameters (hydraulic conductivities) by assimilating dynamic observations in the absence of direct
measurements of the parameters themselves. In some situations, hydraulic conductivity measurements (hard data) may not be
available, which requires the estimation of conductivities from indirect observations, such as piezometric heads. We show
how the NS-EnKF is capable of retrieving the bimodal nature of a synthetic aquifer solely from piezometric head data. By comparison
with a more standard implementation of the EnKF, the NS-EnKF gives better results with regard to histogram preservation, uncertainty
assessment, and transport predictions. 相似文献
14.
In recent years, data assimilation techniques have been applied to an increasingly wider specter of problems. Monte Carlo
variants of the Kalman filter, in particular, the ensemble Kalman filter (EnKF), have gained significant popularity. EnKF
is used for a wide variety of applications, among them for updating reservoir simulation models. EnKF is a Monte Carlo method,
and its reliability depends on the actual size of the sample. In applications, a moderately sized sample (40–100 members)
is used for computational convenience. Problems due to the resulting Monte Carlo effects require a more thorough analysis
of the EnKF. Earlier we presented a method for the assessment of the error emerging at the EnKF update step (Kovalenko et
al., SIAM J Matrix Anal Appl, in press). A particular energy norm of the EnKF error after a single update step was studied.
The energy norm used to assess the error is hard to interpret. In this paper, we derive the distribution of the Euclidean
norm of the sampling error under the same assumptions as before, namely normality of the forecast distribution and negligibility
of the observation error. The distribution depends on the ensemble size, the number and spatial arrangement of the observations,
and the prior covariance. The distribution is used to study the error propagation in a single update step on several synthetic
examples. The examples illustrate the changes in reliability of the EnKF, when the parameters governing the error distribution
vary. 相似文献
15.
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. 相似文献
16.
重质非水相有机污染物(DNAPL)泄漏到地下后,其运移与分布特征受渗透率非均质性影响显著。为刻画DNAPL污染源区结构特征,需进行参数估计以描述水文地质参数的非均质性。本研究构建了基于集合卡尔曼滤波方法(EnKF)与多相流运移模型的同化方案,通过融合DNAPL饱和度观测数据推估非均质介质渗透率空间分布。通过二维砂箱实际与理想算例,验证了同化方法的推估效果,并探讨了不同因素对同化的影响。研究结果表明:基于EnKF方法同化饱和度观测资料可有效地推估非均质渗透率场;参数推估精度随观测时空密度的增大而提高;观测点位置分布对同化效果有所影响,布置在污染集中区域的观测数据对于参数估计具有较高的数据价值。 相似文献
17.
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. 相似文献
18.
The ensemble Kalman filter has been successfully applied for data assimilation in very large models, including those in reservoir
simulation and weather. Two problems become critical in a standard implementation of the ensemble Kalman filter, however,
when the ensemble size is small. The first is that the ensemble approximation to cross-covariances of model and state variables
to data can indicate the presence of correlations that are not real. These spurious correlations give rise to model or state
variable updates in regions that should not be updated. The second problem is that the number of degrees of freedom in the
ensemble is only as large as the size of the ensemble, so the assimilation of large amounts of precise, independent data is
impossible. Localization of the Kalman gain is almost universal in the weather community, but applications of localization
for the ensemble Kalman filter in porous media flow have been somewhat rare. It has been shown, however, that localization
of updates to regions of non-zero sensitivity or regions of non-zero cross-covariance improves the performance of the EnKF
when the ensemble size is small. Localization is necessary for assimilation of large amounts of independent data. The problem
is to define appropriate localization functions for different types of data and different types of variables. We show that
the knowledge of sensitivity alone is not sufficient for determination of the region of localization. The region depends also
on the prior covariance for model variables and on the past history of data assimilation. Although the goal is to choose localization
functions that are large enough to include the true region of non-zero cross-covariance, for EnKF applications, the choice
of localization function needs to balance the harm done by spurious covariance resulting from small ensembles and the harm
done by excluding real correlations. In this paper, we focus on the distance-based localization and provide insights for choosing
suitable localization functions for data assimilation in multiphase flow problems. In practice, we conclude that it is reasonable
to choose localization functions based on well patterns, that localization function should be larger than regions of non-zero
sensitivity and should extend beyond a single well pattern. 相似文献
19.
Elisabeth Peters 《Computational Geosciences》2011,15(2):345-358
The ensemble Kalman filter (EnKF) appears to give good results for matching production data at existing wells. However, the
predictive power of these models outside of the existing wells is much more uncertain. In this paper, for a channelized reservoir
for five different cases with different levels of information the production history is matched using the EnKF. The predictive
power of the resulting model is tested for the existing wells and for new wells. The results show a consistent improvement
for the predictions at the existing wells after assimilation of the production data, but not for prediction of production
at new well locations. The latter depended on the settings of the problem and prior information used. The results also showed
that the fit during the history match was not always a good predictor for predictive capabilities of the history match model.
This suggests that some form of validation outside of observed wells is essential. 相似文献
20.
为减少非恒定水流计算中的不确定性,在水流随机运动系统状态空间模型基础上,应用集合卡尔曼滤波(EnKF)技术建立了非恒定水流分析的实时更新(校正)方法。该方法适用于非线性的随机微分方程,过程和观测噪声可以是非正态分布。同时,为充分利用水位、流量等误差量级相差巨大的观测中所蕴含的有效信息,导出了EnKF多变量分析格式。以明渠单峰洪水过程的合成数据为例,考察了运用建立的实时更新方法分析预报一维洪水演进的性能。重点对比了采用不同精度等级下的水位和流量观测资料进行滤波的效果。在中国现行标准规定的允许观测误差范围内,以水位观测进行一维洪水动力学模型的滤波分析可有效地控制误差、估计流量、识别水流运动系统状态。长江干流清溪场至万县江段实际洪水计算还证实:该方法通过插入即时观测,可实时更新模型状态,给出与实际更为接近的计算。 相似文献