共查询到18条相似文献，搜索用时 296 毫秒
1.
顺序数据同化的Bayes滤波框架 总被引：6，自引：2，他引：4
数据同化是在动力学模型的运行过程中不断融合新的观测信息的方法论,Bayes理论是数据同化的基石.从原理、方法和符号系统为Bayes滤波在数据同化中的应用勾勒一个统一的框架.首先对连续数据同化和顺序数据同化的各种方法做了分类,然后给出了非线性系统顺序数据同化的Bayes递推滤波形式,并在此基础上介绍了典型的顺序数据同化方法--粒子滤波和集合Kalman滤波.粒子滤波实质上是一种基于递推Bayes估计和Monte Carlo模拟的滤波方法,而集合Kalman滤波相当于一种权值相等的粒子滤波.Bayes滤波理论为顺序数据同化提供了更广义的理论框架,从基础的数学理论上揭示了数据同化的基本原理. 相似文献
2.
非线性滤波方法与陆面数据同化 总被引：8，自引：4，他引：4
陆面数据同化研究近几年成为地球科学研究的新兴领域,其中以非线性滤波为代表的数据同化方法发展迅速并得到了广泛应用。在贝叶斯理论框架内,从递推贝叶斯估计理论的角度系统地分析了扩展卡尔曼滤波、无迹卡尔曼滤波、集合卡尔曼滤波、SIR粒子滤波等非线性滤波方法的异同;针对应用比较广泛的集合卡尔曼滤波和SIR粒子滤波应用中存在的问题,论述了几种提高滤波性能的实用方法,如协方差矩阵的Localization方法、协方差矩阵的Inflation方法、双集合卡尔曼滤波方法、扰动集合、扰动大气驱动和模型参数、平方根集合卡尔曼滤波以及粒子滤波算法的改进等。最后总结讨论了各种非线性滤波方法应用中的特点、难点以及各种算法在陆面数据同化中的应用前景和发展方向。 相似文献
3.
4.
不同滤波算法在土壤湿度同化中的应用 总被引：1，自引：0，他引：1
为研究不同滤波算法在土壤湿度同化中的有效性,以及土壤湿度模拟结果对模型参数的敏感性,结合简单生物圈模型SiB2,设置敏感性实验,探求土壤饱和水力传导度对土壤湿度模拟结果的影响;并在此基础上,采用集合卡尔曼滤波（EnKF）、无迹卡尔曼滤波（UKF）和无迹粒子滤波（UPF）开展土壤湿度实时同化实验。结果表明：土壤饱和水力传导度能显著影响土壤湿度模拟精度;利用EnKF、UKF、UPF同化站点观测数据,均能改善土壤湿度模拟结果;3种同化方法在不同土壤层的同化效果不同,在土壤表层,EnKF的有效性优于UKF和UPF,在根域层和土壤深层,3种滤波方法有效性在降雨前后相差较大。因此,针对性地选择同化方法,是提高土壤湿度模拟精度的有效手段。 相似文献
5.
模式时间关联误差对集合平方根滤波估算土壤湿度的影响 总被引：2，自引：1，他引：1
为了定量评估模式时间关联误差对NOAH陆面模式同化表层土壤湿度观测估算土壤湿度廓线的影响,采用集合平方根滤波(En SRF)与状态增广相结合的技术,开展同时更新状态变量和订正模式偏差的观测系统模拟试验,结果表明:同化时若不对存在较大系统性偏差的模式时间关联误差进行处理,En SRF就不能有效估算土壤湿度廓线,而采用状态增广和En SRF相结合的技术,可以在更新土壤湿度时同步订正模式偏差,土壤湿度估算精度明显提高。敏感性试验进一步表明:模式偏差大小、同化时间间隔和观测误差会以不同方式对同化结果造成影响。 相似文献
6.
通过数据同化方法合理地将实时水文观测数据融入到洪水预报模型中,可提高洪水预报模型的实时性和精确度。选取沿程断面流量、水位和糙率系数作为代表水流状态的基本粒子,以监测断面实测水位数据作为观测信息,建立了基于粒子滤波数据同化算法的河道洪水实时概率预报模型。模型应用于黄河中下游河道洪水预报计算的结果表明,采用粒子滤波方法同化观测水位后,不仅可以直接校正水位,同时也可以有效地校正流量和糙率,为未来时刻模型预报计算提供更准确的水流初始条件和糙率取值区间,进而有效地提高模型预报结果的精度,给出合理的概率预报区间。不同预报期的预报结果表明,随着预报期的增长,同化效果减弱,模型预报结果的精度会有所降低,水位概率预报结果受粒子间糙率不同的影响不确定性增加,而流量概率预报结果受给定模型边界条件的影响不确定性降低。所提出模型可以有效同化真实水位观测数据,适合应用于实际的洪水预报工作中。 相似文献
7.
地球化学异常提取的自适应衬值滤波法 总被引：1，自引：0，他引：1
地球化学数据分布,特别是微量元素的频率分布,是正向偏斜的,其异常下限的计算也是基于对数正态分布。当定性分析元素的分布状态时,大部分样本分布状态的结果是可以接受的,但是当定量计算异常下限时,一定要结合具体样本的实际分布情况。引入g-h分布函数对样本数据进行拟合,根据所得到的参数确定异常下限,将异常下限与标准离差分割开来。结合子区中位数衬值滤波法,对个旧及其外围成矿区的化探数据进行处理,同时自适应的调节滤波窗口大小,提取化探异常,从而对预测靶区进行圈定。 相似文献
8.
9.
在集合数据同化过程中,由于远距离的观测与同化状态之间存在着虚假相关,局地化方法受到广泛关注.此外,由于集合数的限制,容易引起欠采样和协方差被低估等现象,使得滤波效果欠佳.因此,提出模糊控制算法,模糊控制算法主要用于判断观测点与状态更新点之间的距离来匹配相应的观测权重,进而调整局地化系数来更新背景误差协方差和观测误差协方差矩阵,从而得到有效的状态估计.基于背景误差协方差局地化方法和观测误差协方差局地化方法,耦合模糊控制,形成了新的算法—模糊控制的背景误差协方差局地化方法和模糊控制的观测误差协方差局地化方法.利用Lorenz-96模型,在小集合数和局地化半径下,得出模糊控制的背景误差协方差局地化方法和模糊控制的观测误差协方差局地化方法有较好的同化性能.通过分析泰勒图谱甄别出新算法与观测点具有高度的相关性以及较小的空间变异性.最后,在不同维数的模糊控制器下,新算法的有效性进一步得到验证.为今后数据同化误差处理方面提供了良好的研究平台. 相似文献
10.
把一根线绕到圆柱体上就形成了一条螺旋线。利用螺旋坐标系统可以把笛卡尔坐标空间中的多维滤波和递归滤波转化到一维空间中来处理。相应地多维递归滤波的稳定性问题就可以在一维空间中进行刻画。本文证明了二维滤波和一维螺旋滤波的等价性,并且给出一维螺旋滤波和递归滤波的算法,实际的算例说明了算法的可靠性和有效性 相似文献
11.
利用1998年7月6日至8月9日青藏高原GAME-Tibet试验区MS3608站点的4cm、20cm和100cm的土壤水分观测数据同化SiB2模型输出的表层、根区和深层土壤水分,探讨了一个基于集合卡尔曼滤波和简单生物圈模型的单点土壤水分同化方案。分析和评价了集合大小、同化周期、模型误差、背景场误差以及观测误差对同化系统性能的影响。结果表明:①增加集合数目可以减小土壤水分同化系统的误差,但同时又降低了运行效率;②对于集合卡尔曼滤波,初始场的估计是否准确对同化系统性能影响不大;③模型误差和观测误差的准确估计可以提高土壤水分的估计精度;④利用数据同化的方法对土壤水分的估计有显著提高。 相似文献
12.
为研究观测资料稀少情况下土壤质地及有机质对土壤水分同化的影响,发展了集合卡尔曼平滑(Ensemble Kalman Smooth, EnKS)的土壤水分同化方案。利用黑河上游阿柔冻融观测站2008年6月1日至10月29日的观测数据,使用EnKS算法将表层土壤水分观测数据同化到简单生物圈模型(Simple Biosphere Model 2, SiB2)中,分析不同方案对土壤水分估计的影响,并与集合卡尔曼滤波算法(EnKF)的结果进行比较。研究结果表明,土壤质地和有机质对表层土壤水分模拟结果影响最大而对深层的影响相对较小;利用EnKF和EnKS算法同化表层土壤水分观测数据,均能够显著提高表层和根区土壤水分估计的精度,EnKS算法的精度略高于EnKF且所受土壤质地和有机质的影响小于EnKF;当观测数据稀少时,EnKS算法仍然可以得到较高精度的土壤水分估计。 相似文献
13.
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. 相似文献
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.
集合卡尔曼滤波(Ensemble Kalman Filter,EnKF)作为一种有效的数据同化方法,在众多数值实验中体现优势的同时,也暴露了它使用小集合估计协方差情况下精度较低的缺陷。为了降低取样噪声对协方差估计的干扰并提高滤波精度,应用局域化函数对小集合估计的协方差进行修正,即在协方差矩阵中以舒尔积的形式增加空间距离权重以限制远距离相关。在一个二维理想孔隙承压含水层模型中的运行结果表明,局域化对集合卡尔曼滤波估计地下水参数的修正十分有效,局域化可以很好地过滤小集合估计中噪声的影响,节省计算量的同时又可以防止滤波发散。相关长度较小的水文地质参数(如对数渗透系数)更容易受到噪声的干扰,更有必要进行局域化修正。 相似文献
16.
The ensemble Kalman filter (EnKF), an efficient data assimilation method showing advantages in many numerical experiments, is deficient when used in approximating covariance from an ensemble of small size. Implicit localization is used to add distance-related weight to covariance and filter spurious correlations which weaken the EnKF??s capability to estimate uncertainty correctly. The effect of this kind of localization is studied in two-dimensional (2D) and three-dimensional (3D) synthetic cases. It is found that EnKF with localization can capture reliably both the mean and variance of the hydraulic conductivity field with higher efficiency; it can also greatly stabilize the assimilation process as a small-size ensemble is used. Sensitivity experiments are conducted to explore the effect of localization function format and filter lengths. It is suggested that too long or too short filter lengths will prevent implicit localization from modifying the covariance appropriately. Steep localization functions will greatly disturb local dynamics like the 0-1 function even if the function is continuous; four relatively gentle localization functions succeed in avoiding obvious disturbance to the system and improve estimation. As the degree of localization of the L function increases, the parameter sensitivity becomes weak, making parameter selection easier, but more information may be lost in the assimilation process. 相似文献
17.
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. 相似文献
18.
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. 相似文献