Utilizing the Ensemble Kalman Filter and Ensemble Kalman Smoother for Combined State and Parameter Estimation of a Three-Dimensional Towed Underwater Cable Model

Author(s):  
Jan Vidar Grindheim ◽  
Inge Revhaug ◽  
Egil Pedersen

A finite difference method (FDM) solving the coupled partial differential equations governing three-dimensional (3D) motions of a towed underwater cable has been implemented in a combined ensemble Kalman filter (EnKF) and ensemble Kalman smoother (EnKS), as a new approach to combined state and parameter estimation for towed underwater cables. A simulation study of the method applied to a seismic streamer has been performed. Cable state variables as well as model parameters are estimated. Parameters estimated are crossline ocean current varying with time as well as cable tangential drag coefficient. The presented results indicate that the method is able to estimate state as well as parameters for seismic streamers.

Author(s):  
Kamalanand Krishnamurthy

Parameter estimation is a central issue in mathematical modelling of biomedical systems and for the development of patient specific models. The technique of estimating parameters helps in obtaining diagnostic information from computational models of biological systems. However, in most of the biomedical systems, the estimation of model parameters is a challenging task due to the nonlinearity of mathematical models. In this chapter, the method of estimation of nonlinear model parameters from measurements of state variables, using the extended Kalman filter, is extensively explained using an example of the three-dimensional model of the HIV/AIDS system.


2018 ◽  
pp. 690-713
Author(s):  
Kamalanand Krishnamurthy

Parameter estimation is a central issue in mathematical modelling of biomedical systems and for the development of patient specific models. The technique of estimating parameters helps in obtaining diagnostic information from computational models of biological systems. However, in most of the biomedical systems, the estimation of model parameters is a challenging task due to the nonlinearity of mathematical models. In this chapter, the method of estimation of nonlinear model parameters from measurements of state variables, using the extended Kalman filter, is extensively explained using an example of the three-dimensional model of the HIV/AIDS system.


2011 ◽  
Vol 15 (8) ◽  
pp. 2437-2457 ◽  
Author(s):  
S. Nie ◽  
J. Zhu ◽  
Y. Luo

Abstract. The performance of the ensemble Kalman filter (EnKF) in soil moisture assimilation applications is investigated in the context of simultaneous state-parameter estimation in the presence of uncertainties from model parameters, soil moisture initial condition and atmospheric forcing. A physically based land surface model is used for this purpose. Using a series of identical twin experiments in two kinds of initial parameter distribution (IPD) scenarios, the narrow IPD (NIPD) scenario and the wide IPD (WIPD) scenario, model-generated near surface soil moisture observations are assimilated to estimate soil moisture state and three hydraulic parameters (the saturated hydraulic conductivity, the saturated soil moisture suction and a soil texture empirical parameter) in the model. The estimation of single imperfect parameter is successful with the ensemble mean value of all three estimated parameters converging to their true values respectively in both NIPD and WIPD scenarios. Increasing the number of imperfect parameters leads to a decline in the estimation performance. A wide initial distribution of estimated parameters can produce improved simultaneous multi-parameter estimation performances compared to that of the NIPD scenario. However, when the number of estimated parameters increased to three, not all parameters were estimated successfully for both NIPD and WIPD scenarios. By introducing constraints between estimated hydraulic parameters, the performance of the constrained three-parameter estimation was successful, even if temporally sparse observations were available for assimilation. The constrained estimation method can reduce RMSE much more in soil moisture forecasting compared to the non-constrained estimation method and traditional non-parameter-estimation assimilation method. The benefit of this method in estimating all imperfect parameters simultaneously can be fully demonstrated when the corresponding non-constrained estimation method displays a relatively poor parameter estimation performance. Because all these constraints between parameters were obtained in a statistical sense, this constrained state-parameter estimation scheme is likely suitable for other land surface models even with more imperfect parameters estimated in soil moisture assimilation applications.


2021 ◽  
Author(s):  
Mengtian Lu ◽  
Sicheng Lu ◽  
Weihong Liao ◽  
Xiaohui Lei ◽  
Zhaokai Yin ◽  
...  

Abstract Although field measurements and using long hydrological datasets provide a reliable method for parameters' calibration, changes in the underlying basin surface and lack of hydrometeorological data may affect parameter accuracy in streamflow simulation. The ensemble Kalman filter (EnKF) can be used as a real-time parameter correction method to solve this problem. In this study, five representative Xin'anjiang model parameters are selected to study the effects of the initial parameter ensemble distribution and the specific function form of the parameter on EnKF parameter estimation process for both single and multiple parameters. Results indicate: (1) the method of parameter calibration to determine the initial distribution mean can improve the assimilation efficiency; (2) there is mutual interference among the parameters during multiple parameters' estimation which invalidates some conclusions of single-parameter estimation. We applied and evaluated the EnKF method in Jinjiang River Basin, China. Compared to traditional approaches, our method showed a better performance in both basins with long hydrometeorological dataset (an increase of Kling–Gupta efficiency (KGE) from 0.810 to 0.887 and a decrease of bias from −1.08% to −0.74%); and in basins with a lack of hydrometeorological data (an increase of KGE from 0.536 to 0.849 and a decrease of bias from −15.55% to −11.42%).


2016 ◽  
Vol 144 (7) ◽  
pp. 2605-2621 ◽  
Author(s):  
Lili Lei ◽  
Jeffrey S. Whitaker

Abstract The analysis produced by the ensemble Kalman filter (EnKF) may be dynamically inconsistent and contain unbalanced gravity waves that are absent in the real atmosphere. These imbalances can be exacerbated by covariance localization and inflation. One strategy to combat the imbalance in the analyses is the incremental analysis update (IAU), which uses the dynamic model to distribute the analyses increments over a time window. The IAU has been widely used in atmospheric and oceanic applications. However, the analysis increment that is gradually introduced during a model integration is often computed once and assumed to be constant for an assimilation window, which can be seen as a three-dimensional IAU (3DIAU). Thus, the propagation of the analysis increment in the assimilation window is neglected, yet this propagation may be important, especially for moving weather systems. To take into account the propagation of the analysis increment during an assimilation window, a four-dimensional IAU (4DIAU) used with the EnKF is presented. It constructs time-varying analysis increments by applying all observations in an assimilation window to state variables at different times during the assimilation window. It then gradually applies these time-varying analysis increments through the assimilation window. Results from a dry two-layer primitive equation model and the NCEP GFS show that EnKF with 4DIAU (EnKF-4DIAU) and 3DIAU (EnKF-3DIAU) reduce imbalances in the analysis compared to EnKF without initialization (EnKF-RAW). EnKF-4DIAU retains the time-varying information in the analysis increments better than EnKF-3DIAU, and produces better analysis and forecast than either EnKF-RAW or EnKF-3DIAU.


SPE Journal ◽  
2010 ◽  
Vol 16 (02) ◽  
pp. 331-342 ◽  
Author(s):  
Hemant A. Phale ◽  
Dean S. Oliver

Summary When the ensemble Kalman filter (EnKF) is used for history matching, the resulting updates to reservoir properties sometimes exceed physical bounds, especially when the problem is highly nonlinear. Problems of this type are often encountered during history matching compositional models using the EnKF. In this paper, we illustrate the problem using an example in which the updated molar density of CO2 in some regions is observed to take negative values while molar densities of the remaining components are increased. Standard truncation schemes avoid negative values of molar densities but do not address the problem of increased molar densities of other components. The results can include a spurious increase in reservoir pressure with a subsequent inability to maintain injection. In this paper, we present a method for constrained EnKF (CEnKF), which takes into account the physical constraints on the plausible values of state variables during data assimilation. In the proposed method, inequality constraints are converted to a small number of equality constraints, which are used as virtual observations for calibrating the model parameters within plausible ranges. The CEnKF method is tested on a 2D compositional model and on a highly heterogeneous three-phase-flow reservoir model. The effect of the constraints on mass conservation is illustrated using a 1D Buckley-Leverett flow example. Results show that the CEnKF technique is able to enforce the nonnegativity constraints on molar densities and the bound constraints on saturations (all phase saturations must be between 0 and 1) and achieve a better estimation of reservoir properties than is obtained using only truncation with the EnKF.


2011 ◽  
Vol 8 (1) ◽  
pp. 1433-1468
Author(s):  
S. Nie ◽  
J. Zhu ◽  
Y. Luo

Abstract. The performance of the ensemble Kalman filter (EnKF) in soil moisture assimilation applications is investigated in the context of simultaneous state-parameter estimation in the presence of uncertainties from model parameters, initial soil moisture condition and atmospheric forcing. A physically-based land surface model is used for this purpose. Using a series of idealized twin experiments, model generated near-surface soil moisture observations are assimilated to estimate soil moisture state and three hydraulic parameters (the saturated hydraulic conductivity, the saturated soil moisture suction and a soil texture empirical parameter) in the model. The single imperfect parameter can be successfully estimated using the EnKF. Results show that all the three estimated parameters converge toward their respective true values, while the root mean squared errors (RMSE) of soil moisture associated with these parameters is on average reduced by 54% and 53% comparing with the non-parameter-estimation benchmark RMSE for near-surface layer and root zone layer, respectively. The performance of simultaneous multi-parameter estimation is significant degraded, mainly because the inherent balance relationship of these parameters is broken and the degree of freedom increases in assimilation processes. By introducing constraints between estimated parameters, the performance of the constraint-based simultaneous multi-parameter estimations are as good as that of single-parameter cases even assimilating temporal-sparse observations. In terms of the relative root mean squared error (RRE), the constraint-based estimation cases can achieve 36% to 53% in near-surface layer and 25% to 50% in root zone layer for different assimilation intervals ranging from 1-day to 40-days. This result suggests that the greatest advantage of this method can be displayed with a proper temporal-sparse assimilation interval of 10-days as actual measurement interval of conventional in situ soil moisture observations. As these obtained constraints are mostly in statistical sense, this constraint-based simultaneous state-parameter estimation scheme is supposed to be suitable for other land surface models in soil moisture assimilation applications.


2017 ◽  
Vol 65 (4) ◽  
pp. 479-488 ◽  
Author(s):  
A. Boboń ◽  
A. Nocoń ◽  
S. Paszek ◽  
P. Pruski

AbstractThe paper presents a method for determining electromagnetic parameters of different synchronous generator models based on dynamic waveforms measured at power rejection. Such a test can be performed safely under normal operating conditions of a generator working in a power plant. A generator model was investigated, expressed by reactances and time constants of steady, transient, and subtransient state in the d and q axes, as well as the circuit models (type (3,3) and (2,2)) expressed by resistances and inductances of stator, excitation, and equivalent rotor damping circuits windings. All these models approximately take into account the influence of magnetic core saturation. The least squares method was used for parameter estimation. There was minimized the objective function defined as the mean square error between the measured waveforms and the waveforms calculated based on the mathematical models. A method of determining the initial values of those state variables which also depend on the searched parameters is presented. To minimize the objective function, a gradient optimization algorithm finding local minima for a selected starting point was used. To get closer to the global minimum, calculations were repeated many times, taking into account the inequality constraints for the searched parameters. The paper presents the parameter estimation results and a comparison of the waveforms measured and calculated based on the final parameters for 200 MW and 50 MW turbogenerators.


2021 ◽  
Vol 11 (7) ◽  
pp. 2898
Author(s):  
Humberto C. Godinez ◽  
Esteban Rougier

Simulation of fracture initiation, propagation, and arrest is a problem of interest for many applications in the scientific community. There are a number of numerical methods used for this purpose, and among the most widely accepted is the combined finite-discrete element method (FDEM). To model fracture with FDEM, material behavior is described by specifying a combination of elastic properties, strengths (in the normal and tangential directions), and energy dissipated in failure modes I and II, which are modeled by incorporating a parameterized softening curve defining a post-peak stress-displacement relationship unique to each material. In this work, we implement a data assimilation method to estimate key model parameter values with the objective of improving the calibration processes for FDEM fracture simulations. Specifically, we implement the ensemble Kalman filter assimilation method to the Hybrid Optimization Software Suite (HOSS), a FDEM-based code which was developed for the simulation of fracture and fragmentation behavior. We present a set of assimilation experiments to match the numerical results obtained for a Split Hopkinson Pressure Bar (SHPB) model with experimental observations for granite. We achieved this by calibrating a subset of model parameters. The results show a steady convergence of the assimilated parameter values towards observed time/stress curves from the SHPB observations. In particular, both tensile and shear strengths seem to be converging faster than the other parameters considered.


Water ◽  
2019 ◽  
Vol 11 (7) ◽  
pp. 1520
Author(s):  
Zheng Jiang ◽  
Quanzhong Huang ◽  
Gendong Li ◽  
Guangyong Li

The parameters of water movement and solute transport models are essential for the accurate simulation of soil moisture and salinity, particularly for layered soils in field conditions. Parameter estimation can be achieved using the inverse modeling method. However, this type of method cannot fully consider the uncertainties of measurements, boundary conditions, and parameters, resulting in inaccurate estimations of parameters and predictions of state variables. The ensemble Kalman filter (EnKF) is well-suited to data assimilation and parameter prediction in Situations with large numbers of variables and uncertainties. Thus, in this study, the EnKF was used to estimate the parameters of water movement and solute transport in layered, variably saturated soils. Our results indicate that when used in conjunction with the HYDRUS-1D software (University of California Riverside, California, CA, USA) the EnKF effectively estimates parameters and predicts state variables for layered, variably saturated soils. The assimilation of factors such as the initial perturbation and ensemble size significantly affected in the simulated results. A proposed ensemble size range of 50–100 was used when applying the EnKF to the highly nonlinear hydrological models of the present study. Although the simulation results for moisture did not exhibit substantial improvement with the assimilation, the simulation of the salinity was significantly improved through the assimilation of the salinity and relative solutetransport parameters. Reducing the uncertainties in measured data can improve the goodness-of-fit in the application of the EnKF method. Sparse field condition observation data also benefited from the accurate measurement of state variables in the case of EnKF assimilation. However, the application of the EnKF algorithm for layered, variably saturated soils with hydrological models requires further study, because it is a challenging and highly nonlinear problem.


Sign in / Sign up

Export Citation Format

Share Document