首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 46 毫秒
1.
This paper presents the application of a multimodel method using a wavelet‐based Kalman filter (WKF) bank to simultaneously estimate decomposed state variables and unknown parameters for real‐time flood forecasting. Applying the Haar wavelet transform alters the state vector and input vector of the state space. In this way, an overall detail plus approximation describes each new state vector and input vector, which allows the WKF to simultaneously estimate and decompose state variables. The wavelet‐based multimodel Kalman filter (WMKF) is a multimodel Kalman filter (MKF), in which the Kalman filter has been substituted for a WKF. The WMKF then obtains M estimated state vectors. Next, the M state‐estimates, each of which is weighted by its possibility that is also determined on‐line, are combined to form an optimal estimate. Validations conducted for the Wu‐Tu watershed, a small watershed in Taiwan, have demonstrated that the method is effective because of the decomposition of wavelet transform, the adaptation of the time‐varying Kalman filter and the characteristics of the multimodel method. Validation results also reveal that the resulting method enhances the accuracy of the runoff prediction of the rainfall–runoff process in the Wu‐Tu watershed. Copyright © 2004 John Wiley & Sons, Ltd.  相似文献   

2.
The least squares estimation procedures used in different disciplines can be classified in four categories:
  • a. Wiener filtering,
  • b. b. Autoregressive estimation,
  • c. c. Kalman filtering,
  • d. d. Recursive least squares estimation.
The recursive least squares estimator is the time average form of the Kalman filter. Likewise, the autoregressive estimator is the time average form of the Wiener filter. Both the Kalman and the Wiener filters use ensemble averages and can basically be constructed without having a particular measurement realisation available. It follows that seismic deconvolution should be based either on autoregression theory or on recursive least squares estimation theory rather than on the normally used Wiener or Kalman theory. A consequence of this change is the need to apply significance tests on the filter coefficients. The recursive least squares estimation theory is particularly suitable for solving the time variant deconvolution problem.  相似文献   

3.
Environmental indices (EI) constitute a common communication tool that is often used to describe the overall status of environmental systems (air, water and soil). EI development entails the use of mathematical operators to aggregate various non-commensurate input parameters in a logical manner. The ordered weighted averaging (OWA) operator is a general mean type operator that provides flexibility in the aggregation process such that the aggregated value is bounded between minimum and maximum values of the input parameters. This flexibility of the OWA operator is realized through the concept of orness, which is a surrogate for decision maker’s attitude. The type of input parameters also affects the choice of aggregation operators. If the input parameters are linguistic or fuzzy, the aggregation through OWA operators is not possible, and the use of fuzzy arithmetic is warranted. The concept of fuzzy number OWA (FN-OWA) operators is explored to handle situations in which one or more input parameter has fuzzy (or linguistic) values. The proposed approach is demonstrated using data provided in an earlier study by Swamee and Tyagi (ASCE J Environ Eng 126(5):451–455, 2000) for establishing water quality indices. Multiple hypothetical scenarios are also generated to highlight the utility and sensitivity of the proposed approach.  相似文献   

4.
The properties of linear spatial interpolators of single realizations and trend components of regionalized variables are examined in this work. In the case of the single realization estimator explicit and exact expressions for the weighting vector and the variances of estimator and estimation error were obtained from a closed-form expression for the inverse of the Lagrangian matrix. The properties of the trend estimator followed directly from the Gauss-Markoff theorem. It was shown that the single realization estimator can be decomposed into two mutually orthogonal random functions of the data, one of which is the trend estimator. The implementation of liear spatial estimation was illustrated with three different methods, i.e., full information maximum likelihood (FIML), restricted maximum likelihood (RML), and Rao's minimum norm invariant quadratic unbiased estimation (MINQUE) for the single realization case and via generalized least squares (GLS) for the trend. The case study involved large correlation length-scale in the covariance of specific yield producing a nested covariance structure that was nearly positive semidefinite. The sensitivity of model parameters, i.e., drift and variance components (local and structured) to the correlation length-scale, choice of covariance model (i.e., exponential and spherical), and estimation method was examined. the same type of sensitivity analysis was conducted for the spatial interpolators. It is interesting that for this case study, characterized by a large correlation length-scale of about 50 mi (80 km), both parameter estimates and linear spatial interpolators were rather insensitive to the choice of covariance model and estimation method within the range of credible values obtained for the correlation length-scale, i.e., 40–60 mi (64–96 km), with alternative estimates falling within ±5% of each other.  相似文献   

5.
 The efficiency of a sequential data assimilation scheme relies on the capability to describe the error covariance. This aspect is all the more relevant if one needs accurate statistics on the estimation error. Frequently an ad hoc function depending on a few parameters is proposed, and these parameters are tuned, estimated or updated. This usually requires that the covariance is second-order stationary (i.e. depends only on the distance between two points). In this paper, we discuss this feature and show that even in simple applications (such as one-dimensional hydrodynamics), this assumption does not hold and may lead to poorly described estimation errors. We propose a method relying on the analysis of the error term and the use of the hydrodynamical model to generate one part of the covariance function, the other part being modeled using a second-order stationary approach. This method is discussed using a twin experiment in the case where a physical parameter is erroneous, and improves significantly the results: the model bias is strongly reduced and the estimation error is well described. Moreover, it enables a better adaptation of the Kalman gain to the actual estimation error.  相似文献   

6.
Sequential kriging and cokriging: Two powerful geostatistical approaches   总被引:1,自引:0,他引:1  
A sequential linear estimator is developed in this study to progressively incorporate new or different spatial data sets into the estimation. It begins with a classical linear estimator (i.e., kriging or cokriging) to estimate means conditioned to a given observed data set. When an additional data set becomes available, the sequential estimator improves the previous estimate by using linearly weighted sums of differences between the new data set and previous estimates at sample locations. Like the classical linear estimator, the weights used in the sequential linear estimator are derived from a system of equations that contains covariances and cross-covariances between sample locations and the location where the estimate is to be made. However, the covariances and cross-covariances are conditioned upon the previous data sets. The sequential estimator is shown to produce the best, unbiased linear estimate, and to provide the same estimates and variances as classic simple kriging or cokriging with the simultaneous use of the entire data set. However, by using data sets sequentially, this new algorithm alleviates numerical difficulties associated with the classical kriging or cokriging techniques when a large amount of data are used. It also provides a new way to incorporate additional information into a previous estimation.  相似文献   

7.
Selection of noise parameters for Kalman filter   总被引:1,自引:1,他引:0  
The Bayesian probabilistic approach is proposed to estimate the process noise and measurement noise parameters for a Kalman filter. With state vectors and covariance matrices estimated by the Kalman filter, the likehood of the measurements can be constructed as a function of the process noise and measurement noise parameters. By maximizing the likelihood function with respect to these noise parameters, the optimal values can be obtained. Furthermore, the Bayesian probabilistic approach allows the associated uncertainty to be quantified. Examples using a single-degree-of-freedom system and a ten-story building illustrate the proposed method. The effect on the performance of the Kalman filter due to the selection of the process noise and measurement noise parameters was demonstrated. The optimal values of the noise parameters were found to be close to the actual values in the sense that the actual parameters were in the region with significant probability density. Through these examples, the Bayesian approach was shown to have the capability to provide accurate estimates of the noise parameters of the Kalman filter, and hence for state estimation.  相似文献   

8.
This paper introduces an extension of the traditional stationary linear coregionalization model to handle the lack of stationarity. Under the proposed model, coregionalization matrices are spatially dependent, and basic univariate spatial dependence structures are non-stationary. A parameter estimation procedure of the proposed non-stationary linear coregionalization model is developed under the local stationarity framework. The proposed estimation procedure is based on the method of moments and involves a matrix-valued local stationary variogram kernel estimator, a weighted local least squares method in combination with a kernel smoothing technique. Local parameter estimates are knitted together for prediction and simulation purposes. The proposed non-stationary multivariate spatial modeling approach is illustrated using two real bivariate data examples. Prediction performance comparison is carried out with the classical stationary multivariate spatial modeling approach. According to several criteria, the prediction performance of the proposed non-stationary multivariate spatial modeling approach appears to be significantly better.  相似文献   

9.
Local seismic event slopes contain subsurface velocity information and can be used to estimate seismic stacking velocity. In this paper, we propose a novel approach to estimate the stacking velocity automatically from seismic reflection data using similarity‐weighted k‐means clustering, in which the weights are local similarity between each trace in common midpoint gather and a reference trace. Local similarity reflects the local signal‐to‐noise ratio in common midpoint gather. We select the data points with high signal‐to‐noise ratio to be used in the velocity estimation with large weights in mapped traveltime and velocity domain by similarity‐weighted k‐means clustering with thresholding. By using weighted k‐means clustering, we make clustering centroids closer to those data points with large weights, which are more reliable and have higher signal‐to‐noise ratio. The interpolation is used to obtain the whole velocity volume after we have got velocity points calculated by weighted k‐means clustering. Using the proposed method, one obtains a more accurate estimate of the stacking velocity because the similarity‐based weighting in clustering takes into account the signal‐to‐noise ratio and reliability of different data points in mapped traveltime and velocity domain. In order to demonstrate that, we apply the proposed method to synthetic and field data examples, and the resulting images are of higher quality when compared with the ones obtained using existing methods.  相似文献   

10.
The key problem in nonparametric frequency analysis of flood and droughts is the estimation of the bandwidth parameter which defines the degree of smoothing. Most of the proposed bandwidth estimators have been based on the density function rather than the cumulative distribution function or the quantile that are the primary interest in frequency analysis. We propose a new bandwidth estimator derived from properties of quantile estimators. The estimator builds on work by Altman and Léger (1995). The estimator is compared to the well-known method of least squares cross-validation (LSCV) using synthetic data generated from various parametric distributions used in hydrologic frequency analysis. Simulations suggest that our estimator performs at least as well as, and in many cases better than, the method of LSCV. In particular, the use of the proposed plug-in estimator reduces bias in the estimation as compared to LSCV. When applied to data sets containing observations with identical values, typically the result of rounding or truncation, the LSCV and most other techniques generally underestimates the bandwidth. The proposed technique performs very well in such situations.  相似文献   

11.
Due to the complicated nature of environmental processes, consideration of uncertainty is an important part of environmental modelling. In this paper, a new variant of the machine learning-based method for residual estimation and parametric model uncertainty is presented. This method is based on the UNEEC-P (UNcertainty Estimation based on local Errors and Clustering – Parameter) method, but instead of multilayer perceptron uses a “fuzzified” version of the general regression neural network (GRNN). Two hydrological models are chosen and the proposed method is used to evaluate their parametric uncertainty. The approach can be classified as a hybrid uncertainty estimation method, and is compared to the group method of data handling (GMDH) and ordinary kriging with linear external drift (OKLED) methods. It is shown that, in terms of inherent complexity, measured by Akaike information criterion (AIC), the proposed fuzzy GRNN method has advantages over other techniques, while its accuracy is comparable. Statistical metrics on verification datasets demonstrate the capability and appropriate efficiency of the proposed method to estimate the uncertainty of environmental models.  相似文献   

12.
为提高弹性波CT层析成像中反演成像的质量,同时考虑到CT走时方程中系数矩阵的病态性和BPT-SIRT算法重建图像质量较低的问题,提出基于主元加权预处理和对BPT-SIRT算法加权的方法。该方法对BPT得到的初始迭代值进行加权预处理,在CT走时方程中的系数矩阵主元上叠加一个权值,再通过非线性的指数加权函数对射线路径矩阵中的元素进行加权。通过数值模拟和预制含缺陷的钢管混凝土试件对该方法进行验证,结果表明本文的方法具有更好的收敛速度和精度,对缺陷的识别效果更好。   相似文献   

13.
This paper analyses the skills of fuzzy computing based rainfall–runoff model in real time flood forecasting. The potential of fuzzy computing has been demonstrated by developing a model for forecasting the river flow of Narmada basin in India. This work has demonstrated that fuzzy models can take advantage of their capability to simulate the unknown relationships between a set of relevant hydrological data such as rainfall and river flow. Many combinations of input variables were presented to the model with varying structures as a sensitivity study to verify the conclusions about the coherence between precipitation, upstream runoff and total watershed runoff. The most appropriate set of input variables was determined, and the study suggests that the river flow of Narmada behaves more like an autoregressive process. As the precipitation is weighted only a little by the model, the last time‐steps of measured runoff are dominating the forecast. Thus a forecast based on expected rainfall becomes very inaccurate. Although good results for one‐step‐ahead forecasts are received, the accuracy deteriorates as the lead time increases. Using the one‐step‐ahead forecast model recursively to predict flows at higher lead time, however, produces better results as opposed to different independent fuzzy models to forecast flows at various lead times. Copyright © 2004 John Wiley & Sons, Ltd.  相似文献   

14.
In this paper, a new state-parameter estimation approach is presented based on the dual ensemble Kalman smoother(DEn KS) and simple biosphere model(Si B2) to sequentially estimate both the soil properties and soil moisture profile by assimilating surface soil moisture observations. The Arou observation station, located in the upper reaches of the Heihe River in northwestern China, was selected to test the proposed method. Three numeric experiments were designed and performed to analyze the influence of uncertainties in model parameters, atmospheric forcing, and the model's physical mechanics on soil moisture estimates. Several assimilation schemes based on the ensemble Kalman filter(En KF), ensemble Kalman smoother(En KS), and dual En KF(DEn KF) were also compared in this study. The results demonstrate that soil moisture and soil properties can be simultaneously estimated by state-parameter estimation methods, which can provide more accurate estimation of soil moisture than traditional filter methods such as En KF and En KS. The estimation accuracy of the model parameters decreased with increasing error sources. DEn KS outperformed DEn KF in estimating soil moisture in most cases, especially where few observations were available. This study demonstrates that the DEn KS approach is a useful and practical way to improve soil moisture estimation.  相似文献   

15.
为减小地震数据缺失给地震后续处理工作带来的影响,需要对地震数据进行插值重建.针对反假频插值重建这个难点问题,进行了相关研究,并由此提出了一种改进的R-P(半径-斜率)域加权反假频地震数据插值重建方法.该方法将F-K(频率-波数)谱变换到R-P域,在R-P域设计一个权函数并将其作用于每次的迭代插值过程.通过模型数据和实际数据的测试,证明了该方法具有较好的反假频插值重建能力.  相似文献   

16.
Under a climate change, the physical factors that influence the rainfall regime are diverse and difficult to predict. The selection of skilful inputs for rainfall forecasting models is, therefore, more challenging. This paper combines wavelet transform and Frank copula function in a mutual information‐based input variable selection (IVS) for non‐linear rainfall forecasting models. The marginal probability density functions (PDFs) of a set of potential rainfall predictors and the rainfall series (predictand) were computed using a wavelet density estimator. The Frank copula function was applied to compute the joint PDF of the predictors and the predictand from their marginal PDFs. The relationship between the rainfall series and the potential predictors was assessed based on the mutual information computed from their marginal and joint PDFs. Finally, the minimum redundancy maximum relevance was used as an IVS stopping criterion to determine the number of skilful input variables. The proposed approach was applied to four stations of the Nigerien Sahel with rainfall series spanning the period 1950–2016 by considering 24 climate indices as potential predictors. Adaptive neuro‐fuzzy inference system, artificial neural networks, and random forest‐based forecast models were used to assess the skill of the proposed IVS method. The three forecasting models yielded satisfactory results, exhibiting a coefficient of determination between 0.52 and 0.69 and a mean absolute percentage error varying from 13.6% to 21%. The adaptive neuro‐fuzzy inference system performed better than the other models at all the stations. A comparison made with KDE‐based mutual information showed the advantage of the proposed wavelet–copula approach.  相似文献   

17.
18.
This paper presents a finite element (FE) model updating procedure applied to complex structures using an eigenvalue sensitivity‐based updating approach. The objective of the model updating is to reduce the difference between the calculated and the measured frequencies. The method is based on the first‐order Taylor‐series expansion of the eigenvalues with respect to some structural parameters selected to be adjusted. These parameters are assumed to be bounded by some prescribed regions which are determined according to the degrees of uncertainty that exist in the parameters. The changes of these parameters are found iteratively by solving a constrained optimization problem. The improvement of the current study is in the use of an objective function that is the sum of a weighted frequency error norm and a weighted perturbation norm of the parameters. Two weighting matrices are introduced to provide flexibility for individual tuning of frequency errors and parameters' perturbations. The proposed method is applied to a 1/150 scaled suspension bridge model. Using 11 measured frequencies as reference, the FE model is updated by adjusting ten selected structural parameters. The final updated FE model for the suspension bridge model is able to produce natural frequencies in close agreement with the measured ones. Copyright © 2000 John Wiley & Sons, Ltd.  相似文献   

19.
基于稀疏反演的地震插值方法是一种重要的插值方法,然而大多数这类方法只针对无噪声数据或者高信噪比数据插值.实际上,地震数据含有各种噪声,使得插值问题变得更加困难.凸集投影方法是一种高效的插值算法,但是对于含噪声数据的插值效果不理想,针对含噪声数据提出的加权凸集投影方法能够实现同时插值和去噪,但是除了最小阈值需要认真选取外,增加一个权重因子来实现去噪功能.本文由迭代阈值算法推导出加权凸集投影方法,证明其是解无约束优化问题的一种方法,加权因子可以看作拟合误差项的系数.本文还提出了一种改进的凸集投影方法,与原始凸集投影方法相比该方法不需要增加任何计算量,只要通过阈值的选择来进行插值和去噪.数值模拟证明了该算法的计算效率,并且对含噪声数据能够实现较好的插值效果;先插值后去噪的结果证明了同时去噪和插值算法的可靠性和稳定性.  相似文献   

20.
Studies have illustrated the performance of at-site and regional flood quantile estimators. For realistic generalized extreme value (GEV) distributions and short records, a simple index-flood quantile estimator performs better than two-parameter (2P) GEV quantile estimators with probability weighted moment (PWM) estimation using a regional shape parameter and at-site mean and L-coefficient of variation (L-CV), and full three-parameter at-site GEV/PWM quantile estimators. However, as regional heterogeneity or record lengths increase, the 2P-estimator quickly dominates. This paper generalizes the index flood procedure by employing regression with physiographic information to refine a normalized T-year flood estimator. A linear empirical Bayes estimator uses the normalized quantile regression estimator to define a prior distribution which is employed with the normalized 2P-quantile estimator. Monte Carlo simulations indicate that this empirical Bayes estimator does essentially as well as or better than the simpler normalized quantile regression estimator at sites with short records, and performs as well as or better than the 2P-estimator at sites with longer records or smaller L-CV.  相似文献   

设为首页 | 免责声明 | 关于勤云 | 加入收藏

Copyright©北京勤云科技发展有限公司  京ICP备09084417号