In order to detect the deformation in real-time of the GPS time series and improve its reliability, the multiple Kalman filters model with shaping filter was proposed. Two problems were solved: firstly, because the GP...In order to detect the deformation in real-time of the GPS time series and improve its reliability, the multiple Kalman filters model with shaping filter was proposed. Two problems were solved: firstly, because the GPS real-time deformation series with a high sampling rate contain coloured noise, the multiple Kalman filter model requires the white noise, and the multiple Kalman filters model is augmented by a shaping filter in order to reduce the colored noise; secondly, the multiple Kalman filters model with shaping filter can detect the deformation epoch in real-time and improve the quality of GPS measurements for the real-time deformation applications. Based on the comparisons of the applications in different GPS time series with different models, the advantages of the proposed model were illustrated. The proposed model can reduce the colored noise, detect the smaller changes, and improve the precision of the detected deformation epoch.展开更多
During geodetic monitoring with GNSS technology one of important steps is the correct processing and analysis of the measured displacements. We used the processing method of Kalman filter smoothing algorithm, which al...During geodetic monitoring with GNSS technology one of important steps is the correct processing and analysis of the measured displacements. We used the processing method of Kalman filter smoothing algorithm, which allows to evaluate not only displacements, but also the speed, acceleration, and other characteristics of the deformation model. One of the important issues is the calculation of the obser- vations weight matrix in the Kalman filter. Recurrence algorithm of Kalman filtering can calculate and specify the weights during processing. However, the weights obtained in such way do not always exactly correspond to the actual observation accuracy. We established the observations weights based on the accuracy of baseline measurements. In the presented study, we offered and investigated different models of establishing the accuracy of the baselines. The offered models and the processing of the measured displacements were tested on an experimentally geodetic GNSS network. The research results show that despite of different weight models, changing weights up to 2 times do not change Kalman filtering ac- curacy extremely. The significant improvements for Kalman filtering accuracy for baselines shorter than 10 km were not got. Therefore, for typical GNSS monitoring networks with baseline range 10-15 km, we recommend to use any kind of models. The compulsory condition for getting correct and reliable results is checking results on blunders. For baselines, which are longer than 15 km we propose to use weight model which include baseline standard deviation from network adjustment and corrections for baseline length and its accuracy.展开更多
The design problem of the state filter for the generalized stochastic 2-D Roesser models, which appears when both the state and measurement are simultaneously subjected to the interference from white noise, is discuss...The design problem of the state filter for the generalized stochastic 2-D Roesser models, which appears when both the state and measurement are simultaneously subjected to the interference from white noise, is discussed. The well-known Kalman filter design is extended to the generalized 2-D Roesser models. Based on the method of “scanning line by line”,the filtering problem of generalized 2-D Roesser models with mode-energy reconstruction is solved. The formula of the optimal filtering, which minimizes the variance of the estimation error of the state vectors, is derived. The validity of the designed filter is verified by the calculation steps and the examples are introduced.展开更多
Data assimilation in agricultural remote sensing research is of great significance to integrate with remote sensing observations and model simulations for parameters estimation. The present investigation not only desi...Data assimilation in agricultural remote sensing research is of great significance to integrate with remote sensing observations and model simulations for parameters estimation. The present investigation not only designed and realized the Ensemble Kalman Filtering algorithm (EnKF) assimilation by combing the crop growth model (CERES-Wheat) with remote sensing data, but also optimized and updated the key parameters (LAI) of winter wheat by using remote sensing data. Results showed that the assimilation LAI and the observation ones agreed with each other, and the R2 reached 0.8315. So assimilation remote sensing and crop model could provide reference data for the agricultural production.展开更多
Kalman filter data assimilation technique is incorporated into a standard two-dimensional linear storm surge model. Imperfect model equation and imperfect meteorological forcimg are accounted for by adding noise terms...Kalman filter data assimilation technique is incorporated into a standard two-dimensional linear storm surge model. Imperfect model equation and imperfect meteorological forcimg are accounted for by adding noise terms to the momentum equations. The deterministic model output is corrected by using the available tidal gauge station data. The stationary Kalman filter algorithm for the model domain is calculated by an iterative procedure using specified information on the inaccuracies in the momentum e- quations and specified error information for the observations. An application to a real storm surge that occurred in the summer of 1956 in the East China Sea is performed by means of this data assimilation technique. The result shows that Kalman filter is useful for storm surge forecast and hindcast.展开更多
The computational cost required by the Ensemble Kalman Filter (EnKF) is much larger than that of some simpler assimilation schemes, such as Optimal Interpolation (OI) or three-dimension variational as- similation ...The computational cost required by the Ensemble Kalman Filter (EnKF) is much larger than that of some simpler assimilation schemes, such as Optimal Interpolation (OI) or three-dimension variational as- similation (3DVAR). Ensemble optimal interpolation (EnOI), a crudely simplified implementation of EnKF, is sometimes used as a substitute in some oceanic applications and requires much less computational time than EnKF. In this paper, to compromise between computational cost and dynamic covariance, we use the idea of "dressing" a small size dynamical ensemble with a larger number of static ensembles in order to form an approximate dynamic covariance. The term "dressing" means that a dynamical ensemble seed from model runs is perturbed by adding the anomalies of some static ensembles. This dressing EnKF (DrEnKF for short) scheme is tested in assimilation of real altimetry data in the Pacific using the HYbrid Coordinate Ocean Model (HYCOM) over a four-year period. Ten dynamical ensemble seeds are each dressed by 10 static ensemble members selected from a 100-member static ensemble. Results are compared to two EnKF assimilation runs that use 10 and 100 dynamical ensemble members. Both temperature and salinity fields from the DrEnKF and the EnKF are compared to observations from Argo floats and an OI SST dataset. The results show that the DrEnKF and the 100-member EnKF yield similar root mean square errors (RMSE) at every model level. Error covariance matrices from the DrEnKF and the 100-member EnKF are also compared and show good agreement.展开更多
To accurately estimate winter wheat yields and analyze the uncertainty in crop model data assimilations, winter wheat yield estimates were obtained by assimilating measured or remotely sensed leaf area index (LAI) v...To accurately estimate winter wheat yields and analyze the uncertainty in crop model data assimilations, winter wheat yield estimates were obtained by assimilating measured or remotely sensed leaf area index (LAI) values. The performances of the calibrated crop environment resource synthesis for wheat (CERES-Wheat) model for two different assimilation scenarios were compared by employing ensemble Kalman filter (EnKF)-based strategies. The uncertainty factors of the crop model data assimilation was analyzed by considering the observation errors, assimilation stages and temporal-spatial scales. Overalll the results indicated a better yield estimate performance when the EnKF-based strategy was used to comprehen- sively consider several factors in the initial conditions and observations. When using this strategy, an adjusted coefficients of determination (R2) of 0.84, a root mean square error (RMSE) of 323 kg ha-1, and a relative errors (RE) of 4.15% were obtained at the field plot scale and an R2 of 0.81, an RMSE of 362 kg ha-1, and an RE of 4.52% were obtained at the pixel scale of 30 mx30 m. With increasing observation errors, the accuracy of the yield estimates obviously decreased, but an acceptable estimate was observed when the observation errors were within 20%. Winter wheat yield estimates could be improved significantly by assimilating observations from the middle to the end of the crop growing seasons. With decreasing assimilation frequency and pixel resolution, the accuracy of the crop yield estimates decreased; however, the computation time decreased. It is important to consider reasonable temporal-spatial scales and assimilation stages to obtain tradeoffs between accuracy and computation time, especially in operational systems used for regional crop yield estimates.展开更多
Aiming at improving the estimation accuracy and real-time of nonlinear system with linear Gaussian sub-structure,a novel marginalized cubature Kalman filter is proposed in Bayesian estimation framework. Firstly,the ma...Aiming at improving the estimation accuracy and real-time of nonlinear system with linear Gaussian sub-structure,a novel marginalized cubature Kalman filter is proposed in Bayesian estimation framework. Firstly,the marginalized technique is adopted to model the target system dynamics with nonlinear state and linear state separately,and the two parts are estimated by cubature Kalman filter and standard Kalman filter respectively. Therefore,the linear part avoids the generation and propagation process of cubature points. Accordingly,the computational complexity is reduced.Meanwhile,the accuracy of state estimation is improved by taking the difference of nonlinear state estimation as the measurement of linear state. Furthermore,the computational complexity of marginalized cubature Kalman filter is discussed by calculating the number of floating-point operation. Finally,simulation experiments and analysis show that the proposed algorithm can improve the performance of filtering precision and real-time effectively in target tracking system.展开更多
This study investigated the feasibility of updating prior uncertain geologic models using Ensemble Kalman filter for controlling water coning problems in horizontal wells. Current downhole data acquisition technol-ogy...This study investigated the feasibility of updating prior uncertain geologic models using Ensemble Kalman filter for controlling water coning problems in horizontal wells. Current downhole data acquisition technol-ogy allows continuous updating of the reservoir models and real-time control of well operations. Ensemble Kalman Filter is a model updating algorithm that permits rapid assimilation of production response for res-ervoir model updating and uncertainty assessment. The effect of the type and amount of production data on the updated geologic models was investigated first through a synthetic reservoir model, and then imple-mented on a laboratory experiment that simulated the production of a horizontal well affected by water con-ing. The worth of periodic model updating for optimized production and oil recovery is demonstrated.展开更多
A disease transmission model of SEIR type is discussed in a stochastic point of view. We start by formulating the SEIR epidemic model in form of a system of nonlinear differential equations and then change it to a sys...A disease transmission model of SEIR type is discussed in a stochastic point of view. We start by formulating the SEIR epidemic model in form of a system of nonlinear differential equations and then change it to a system of nonlinear stochastic differential equations (SDEs). The numerical simulation of the resulting SDEs is done by Euler-Maruyama scheme and the parameters are estimated by adaptive Markov chain Monte Carlo and extended Kalman filter methods. The stochastic results are discussed and it is observed that with the SDE type of modeling, the parameters are also identifiable.展开更多
This paper introduces a stochastic hemodynamic system to describe the brain neural activity based on the balloon model. A continuous-discrete extended Kalman filter is used to estimate the nonlinear model states. The ...This paper introduces a stochastic hemodynamic system to describe the brain neural activity based on the balloon model. A continuous-discrete extended Kalman filter is used to estimate the nonlinear model states. The stability, controllability and observability of the proposed model are described based on the simulation and measurement data analysis. The observability and controllability characteristics are in- troduced as significant factors to validate the preference of different hemodynamic factors to be considered for diagnosis and monitoring in clinical applications. This model also can be efficiently applied in any monitoring and control platform include brain and for study of hemodynamics in brain imaging modalities such as pulse oximetry and functional near infrared spectroscopy. The work is on progress to extend the proposed model to cover more hemodynamic and neural brain signals for real-time in-vivo application.展开更多
基金Project(20120022120011)supported by the Specialized Research Fund for the Doctoral Program of Higher Education of ChinaProject(2652012062)supported by the Fundamental Research Funds for the Central Universities,China
文摘In order to detect the deformation in real-time of the GPS time series and improve its reliability, the multiple Kalman filters model with shaping filter was proposed. Two problems were solved: firstly, because the GPS real-time deformation series with a high sampling rate contain coloured noise, the multiple Kalman filter model requires the white noise, and the multiple Kalman filters model is augmented by a shaping filter in order to reduce the colored noise; secondly, the multiple Kalman filters model with shaping filter can detect the deformation epoch in real-time and improve the quality of GPS measurements for the real-time deformation applications. Based on the comparisons of the applications in different GPS time series with different models, the advantages of the proposed model were illustrated. The proposed model can reduce the colored noise, detect the smaller changes, and improve the precision of the detected deformation epoch.
文摘During geodetic monitoring with GNSS technology one of important steps is the correct processing and analysis of the measured displacements. We used the processing method of Kalman filter smoothing algorithm, which allows to evaluate not only displacements, but also the speed, acceleration, and other characteristics of the deformation model. One of the important issues is the calculation of the obser- vations weight matrix in the Kalman filter. Recurrence algorithm of Kalman filtering can calculate and specify the weights during processing. However, the weights obtained in such way do not always exactly correspond to the actual observation accuracy. We established the observations weights based on the accuracy of baseline measurements. In the presented study, we offered and investigated different models of establishing the accuracy of the baselines. The offered models and the processing of the measured displacements were tested on an experimentally geodetic GNSS network. The research results show that despite of different weight models, changing weights up to 2 times do not change Kalman filtering ac- curacy extremely. The significant improvements for Kalman filtering accuracy for baselines shorter than 10 km were not got. Therefore, for typical GNSS monitoring networks with baseline range 10-15 km, we recommend to use any kind of models. The compulsory condition for getting correct and reliable results is checking results on blunders. For baselines, which are longer than 15 km we propose to use weight model which include baseline standard deviation from network adjustment and corrections for baseline length and its accuracy.
基金National Nature Science Foundation of China Under Grants (60474078 ,60304001 ,60574015)
文摘The design problem of the state filter for the generalized stochastic 2-D Roesser models, which appears when both the state and measurement are simultaneously subjected to the interference from white noise, is discussed. The well-known Kalman filter design is extended to the generalized 2-D Roesser models. Based on the method of “scanning line by line”,the filtering problem of generalized 2-D Roesser models with mode-energy reconstruction is solved. The formula of the optimal filtering, which minimizes the variance of the estimation error of the state vectors, is derived. The validity of the designed filter is verified by the calculation steps and the examples are introduced.
基金supported by the National Natural Science Foundation of China (40701120)the Beijing Natural Science Foundation, China (4092016)the Beijing Nova, China (2008B33)
文摘Data assimilation in agricultural remote sensing research is of great significance to integrate with remote sensing observations and model simulations for parameters estimation. The present investigation not only designed and realized the Ensemble Kalman Filtering algorithm (EnKF) assimilation by combing the crop growth model (CERES-Wheat) with remote sensing data, but also optimized and updated the key parameters (LAI) of winter wheat by using remote sensing data. Results showed that the assimilation LAI and the observation ones agreed with each other, and the R2 reached 0.8315. So assimilation remote sensing and crop model could provide reference data for the agricultural production.
基金This project was supported by the National Natural Science Foundation of China under contract No. 40176001 and theHi-tech Rese
文摘Kalman filter data assimilation technique is incorporated into a standard two-dimensional linear storm surge model. Imperfect model equation and imperfect meteorological forcimg are accounted for by adding noise terms to the momentum equations. The deterministic model output is corrected by using the available tidal gauge station data. The stationary Kalman filter algorithm for the model domain is calculated by an iterative procedure using specified information on the inaccuracies in the momentum e- quations and specified error information for the observations. An application to a real storm surge that occurred in the summer of 1956 in the East China Sea is performed by means of this data assimilation technique. The result shows that Kalman filter is useful for storm surge forecast and hindcast.
基金supported by the Knowledge Innovation Program of Chinese Academy of Sciences (Grant No. KZCX1-YW-12-03)National Basic Research Program of China (2006CB403600)+3 种基金Project of Young Scientists Fund by National Natural Sciences Foundation of China (Grant No. 40606008)National Science and Technology Infrastructure Program(2006BAC03B04)supported by National Natural Sciences Foundation of China (Grant No.40531006)supported by a private donation from Trond Mohn c/o Frank Mohn AS, Bergenand the MERSEA project from the European Commission (Grant No. SIP3-CT-2003-502885)
文摘The computational cost required by the Ensemble Kalman Filter (EnKF) is much larger than that of some simpler assimilation schemes, such as Optimal Interpolation (OI) or three-dimension variational as- similation (3DVAR). Ensemble optimal interpolation (EnOI), a crudely simplified implementation of EnKF, is sometimes used as a substitute in some oceanic applications and requires much less computational time than EnKF. In this paper, to compromise between computational cost and dynamic covariance, we use the idea of "dressing" a small size dynamical ensemble with a larger number of static ensembles in order to form an approximate dynamic covariance. The term "dressing" means that a dynamical ensemble seed from model runs is perturbed by adding the anomalies of some static ensembles. This dressing EnKF (DrEnKF for short) scheme is tested in assimilation of real altimetry data in the Pacific using the HYbrid Coordinate Ocean Model (HYCOM) over a four-year period. Ten dynamical ensemble seeds are each dressed by 10 static ensemble members selected from a 100-member static ensemble. Results are compared to two EnKF assimilation runs that use 10 and 100 dynamical ensemble members. Both temperature and salinity fields from the DrEnKF and the EnKF are compared to observations from Argo floats and an OI SST dataset. The results show that the DrEnKF and the 100-member EnKF yield similar root mean square errors (RMSE) at every model level. Error covariance matrices from the DrEnKF and the 100-member EnKF are also compared and show good agreement.
基金supported by the National Natural Science Foundation of China (41401491,41371396,41301457,41471364)the Introduction of International Advanced Agricultural Science and Technology,Ministry of Agriculture,China (948 Program,2016-X38)+1 种基金the Agricultural Scientific Research Fund of Outstanding Talentsthe Open Fund for the Key Laboratory of Agri-informatics,Ministry of Agriculture,China (2013009)
文摘To accurately estimate winter wheat yields and analyze the uncertainty in crop model data assimilations, winter wheat yield estimates were obtained by assimilating measured or remotely sensed leaf area index (LAI) values. The performances of the calibrated crop environment resource synthesis for wheat (CERES-Wheat) model for two different assimilation scenarios were compared by employing ensemble Kalman filter (EnKF)-based strategies. The uncertainty factors of the crop model data assimilation was analyzed by considering the observation errors, assimilation stages and temporal-spatial scales. Overalll the results indicated a better yield estimate performance when the EnKF-based strategy was used to comprehen- sively consider several factors in the initial conditions and observations. When using this strategy, an adjusted coefficients of determination (R2) of 0.84, a root mean square error (RMSE) of 323 kg ha-1, and a relative errors (RE) of 4.15% were obtained at the field plot scale and an R2 of 0.81, an RMSE of 362 kg ha-1, and an RE of 4.52% were obtained at the pixel scale of 30 mx30 m. With increasing observation errors, the accuracy of the yield estimates obviously decreased, but an acceptable estimate was observed when the observation errors were within 20%. Winter wheat yield estimates could be improved significantly by assimilating observations from the middle to the end of the crop growing seasons. With decreasing assimilation frequency and pixel resolution, the accuracy of the crop yield estimates decreased; however, the computation time decreased. It is important to consider reasonable temporal-spatial scales and assimilation stages to obtain tradeoffs between accuracy and computation time, especially in operational systems used for regional crop yield estimates.
基金Supported by the National Natural Science Foundation of China(No.61771006)the Open Foundation of Key Laboratory of Spectral Imaging Technology of the Chinese Academy of Sciences(No.LSIT201711D)+1 种基金the Outstanding Young Cultivation Foundation of Henan University(No.0000A40366)the Excellent Chinese and Foreign Youth Exchange Programme of China Science and Technology Association(2017CASTQNJL046)
文摘Aiming at improving the estimation accuracy and real-time of nonlinear system with linear Gaussian sub-structure,a novel marginalized cubature Kalman filter is proposed in Bayesian estimation framework. Firstly,the marginalized technique is adopted to model the target system dynamics with nonlinear state and linear state separately,and the two parts are estimated by cubature Kalman filter and standard Kalman filter respectively. Therefore,the linear part avoids the generation and propagation process of cubature points. Accordingly,the computational complexity is reduced.Meanwhile,the accuracy of state estimation is improved by taking the difference of nonlinear state estimation as the measurement of linear state. Furthermore,the computational complexity of marginalized cubature Kalman filter is discussed by calculating the number of floating-point operation. Finally,simulation experiments and analysis show that the proposed algorithm can improve the performance of filtering precision and real-time effectively in target tracking system.
文摘This study investigated the feasibility of updating prior uncertain geologic models using Ensemble Kalman filter for controlling water coning problems in horizontal wells. Current downhole data acquisition technol-ogy allows continuous updating of the reservoir models and real-time control of well operations. Ensemble Kalman Filter is a model updating algorithm that permits rapid assimilation of production response for res-ervoir model updating and uncertainty assessment. The effect of the type and amount of production data on the updated geologic models was investigated first through a synthetic reservoir model, and then imple-mented on a laboratory experiment that simulated the production of a horizontal well affected by water con-ing. The worth of periodic model updating for optimized production and oil recovery is demonstrated.
文摘A disease transmission model of SEIR type is discussed in a stochastic point of view. We start by formulating the SEIR epidemic model in form of a system of nonlinear differential equations and then change it to a system of nonlinear stochastic differential equations (SDEs). The numerical simulation of the resulting SDEs is done by Euler-Maruyama scheme and the parameters are estimated by adaptive Markov chain Monte Carlo and extended Kalman filter methods. The stochastic results are discussed and it is observed that with the SDE type of modeling, the parameters are also identifiable.
文摘This paper introduces a stochastic hemodynamic system to describe the brain neural activity based on the balloon model. A continuous-discrete extended Kalman filter is used to estimate the nonlinear model states. The stability, controllability and observability of the proposed model are described based on the simulation and measurement data analysis. The observability and controllability characteristics are in- troduced as significant factors to validate the preference of different hemodynamic factors to be considered for diagnosis and monitoring in clinical applications. This model also can be efficiently applied in any monitoring and control platform include brain and for study of hemodynamics in brain imaging modalities such as pulse oximetry and functional near infrared spectroscopy. The work is on progress to extend the proposed model to cover more hemodynamic and neural brain signals for real-time in-vivo application.