An ensemble Kalman filter approach to parameter estimation for patient-specific cardiovascular flow modeling

Many previous studies have shown that the fidelity of three-dimensional cardiovascular flow simulations depends strongly on inflow and outflow boundary conditions that accurately describe the characteristics of the larger vascular network. These boundary conditions are generally based on lower-dimensional models that represent the upstream or downstream flow behavior in some aggregated fashion. However, the parameters of these models are patient-specific, and no clear technique exists for determining them. In this work, an ensemble Kalman filter (EnKF) is implemented for the purpose of estimating parameters in cardiovascular models through the assimilation of specific patients’ clinical measurements. Two types of models are studied: a fully zero-dimensional model of the right heart and pulmonary circulation, and a coupled 0D–1D model of the lower leg. Model parameters are estimated using measurements from both healthy and hypertensive patients, and demonstrate that the EnKF is able to generate distinct parameter sets whose model predictions produce features unique to each measurement set. Attention is also given toward the quality of model predictions made in the absence of direct clinical counterparts, as well as techniques to improve filter robustness against shrinking ensemble covariance.

[1]  Pablo J. Blanco,et al.  Partitioned Analysis for Dimensionally-Heterogeneous Hydraulic Networks , 2011, Multiscale Model. Simul..

[2]  M Danielsen,et al.  Modeling ventricular contraction with heart rate changes. , 2003, Journal of theoretical biology.

[3]  R. Srinivasan,et al.  A global sensitivity analysis tool for the parameters of multi-variable catchment models , 2006 .

[4]  Berend E. Westerhof,et al.  The arterial Windkessel , 2009, Medical & Biological Engineering & Computing.

[5]  Geir Evensen,et al.  The ensemble Kalman filter for combined state and parameter estimation: MONTE CARLO TECHNIQUES FOR DATA ASSIMILATION IN LARGE SYSTEMS , 2009 .

[6]  R. Himeno,et al.  Biomechanical characterization of ventricular-arterial coupling during aging: a multi-scale model study. , 2009, Journal of biomechanics.

[7]  D. Boyers,et al.  Simulation of the human cardiovascular system: a model with normal responses to change of posture, blood loss, transfusion, and autonomic blockade , 1972 .

[8]  T Kööbi,et al.  Pulse wave velocity reference values in healthy adults aged 26–75 years , 2007, Clinical physiology and functional imaging.

[9]  Horst Olschewski,et al.  Magnetic Resonance–Derived 3-Dimensional Blood Flow Patterns in the Main Pulmonary Artery as a Marker of Pulmonary Hypertension and a Measure of Elevated Mean Pulmonary Arterial Pressure , 2008, Circulation. Cardiovascular imaging.

[10]  Michael Danielsen,et al.  Modeling of feedback mechanisms which control the heart function in a view to an implementation in cardiovascular models , 1998 .

[11]  Lucas O Müller,et al.  A global multiscale mathematical model for the human circulation with emphasis on the venous system , 2014, International journal for numerical methods in biomedical engineering.

[12]  D. F. Young,et al.  Computer simulation of arterial flow with applications to arterial and aortic stenoses. , 1992, Journal of biomechanics.

[13]  W. Yeh Review of Parameter Identification Procedures in Groundwater Hydrology: The Inverse Problem , 1986 .

[14]  S. Osher,et al.  Efficient implementation of essentially non-oscillatory shock-capturing schemes,II , 1989 .

[15]  Jeffrey R. Gohean A closed-loop multi-scale model of the cardiovascular system for evaluation of ventricular assist devices , 2007 .

[16]  G. Evensen The ensemble Kalman filter for combined state and parameter estimation , 2009, IEEE Control Systems.

[17]  P J Blanco,et al.  On the integration of the baroreflex control mechanism in a heterogeneous model of the cardiovascular system , 2012, International journal for numerical methods in biomedical engineering.

[18]  S. Sherwin,et al.  Lumped parameter outflow models for 1-D blood flow simulations: Effect on pulse waves and parameter estimation , 2008 .

[19]  S. Maka,et al.  Sensitivity analysis and parameter estimation of cardiovascular model , 2016, 2016 International Conference on Systems in Medicine and Biology (ICSMB).

[20]  V. Fuster,et al.  Characterization of right ventricular remodeling and failure in a chronic pulmonary hypertension model. , 2014, American journal of physiology. Heart and circulatory physiology.

[21]  Yuguo Chen,et al.  State and parameter estimation of hydrologic models using the constrained ensemble Kalman filter , 2009 .

[22]  O. Frank,et al.  Die grundform des arteriellen pulses , 1899 .

[23]  Charles A. Taylor,et al.  Outflow boundary conditions for one-dimensional finite element modeling of blood flow and pressure waves in arteries , 2004 .

[24]  Marwan A. Simaan,et al.  Estimation of systemic vascular bed parameters for artificial heart control , 1998 .

[25]  Michel Verhaegen,et al.  Estimation of Three- and Four-Element Windkessel Parameters Using Subspace Model Identification , 2010, IEEE Transactions on Biomedical Engineering.

[26]  T D Nauser,et al.  Diagnosis and treatment of pulmonary hypertension. , 2001, American family physician.

[27]  L. Formaggia,et al.  Numerical modeling of 1D arterial networks coupled with a lumped parameters description of the heart , 2006, Computer methods in biomechanics and biomedical engineering.

[28]  Soroosh Sorooshian,et al.  Dual state-parameter estimation of hydrological models using ensemble Kalman filter , 2005 .

[29]  J K Raines,et al.  A computer simulation of arterial dynamics in the human leg. , 1974, Journal of biomechanics.

[30]  Matteo Semplice,et al.  On a third order CWENO boundary treatment with application to networks of hyperbolic conservation laws , 2018, Appl. Math. Comput..

[31]  Alison L. Marsden,et al.  A modular numerical method for implicit 0D/3D coupling in cardiovascular finite element simulations , 2013, J. Comput. Phys..

[32]  S. Sherwin,et al.  One-dimensional modelling of a vascular network in space-time variables , 2003 .

[33]  Horst Olschewski,et al.  Blood flow vortices along the main pulmonary artery measured with MR imaging for diagnosis of pulmonary hypertension. , 2015, Radiology.

[34]  J. Whitaker,et al.  Evaluating Methods to Account for System Errors in Ensemble Data Assimilation , 2012 .

[35]  Mohammadali Tarrahi,et al.  Assessing the performance of the ensemble Kalman filter for subsurface flow data integration under variogram uncertainty , 2011 .

[36]  J P Mynard,et al.  A simple, versatile valve model for use in lumped parameter and one‐dimensional cardiovascular models , 2012, International journal for numerical methods in biomedical engineering.

[37]  P. Nardinocchi,et al.  A Heterogeneous Approach for Modelling Blood Flow in an Arterial Segment , 2003 .

[38]  Geir Evensen,et al.  The Ensemble Kalman Filter: theoretical formulation and practical implementation , 2003 .

[39]  Justin S Tran,et al.  Automated Tuning for Parameter Identification and Uncertainty Quantification in Multi-scale Coronary Simulations. , 2017, Computers & fluids.

[40]  Charles A. Taylor,et al.  On Coupling a Lumped Parameter Heart Model and a Three-Dimensional Finite Element Aorta Model , 2009, Annals of Biomedical Engineering.

[41]  G. Evensen,et al.  Analysis Scheme in the Ensemble Kalman Filter , 1998 .

[42]  J AlastrueyArimon Numerical modelling of pulse wave propagation in the cardiovascular system : development, validation and clinical applications. , 2006 .

[43]  J.R. Boston,et al.  Control of heart assist devices , 2003, 42nd IEEE International Conference on Decision and Control (IEEE Cat. No.03CH37475).

[44]  R J Cohen,et al.  System identification of closed-loop cardiovascular control: effects of posture and autonomic blockade. , 1997, The American journal of physiology.

[45]  N. Stergiopulos,et al.  Validation of a one-dimensional model of the systemic arterial tree. , 2009, American journal of physiology. Heart and circulatory physiology.

[46]  H. Tran,et al.  Blood pressure and blood flow variation during postural change from sitting to standing: model development and validation. , 2005, Journal of applied physiology.

[47]  Alastair Gillies,et al.  Better Management of Cardiovascular Diseases by Pulse Wave Velocity: Combining Clinical Practice with Clinical Research using Evidence-Based Medicine , 2007, Clinical Medicine & Research.

[48]  Charles A. Taylor,et al.  Patient-specific modeling of cardiovascular mechanics. , 2009, Annual review of biomedical engineering.

[49]  Ryutaro Himeno,et al.  Multi-scale modeling of the human cardiovascular system with applications to aortic valvular and arterial stenoses , 2009, Medical & Biological Engineering & Computing.

[50]  A Noordergraaf,et al.  Analog studies of the human systemic arterial tree. , 1969, Journal of biomechanics.

[51]  A. Shoukas,et al.  Load Independence of the Instantaneous Pressure‐Volume Ratio of the Canine Left Ventricle and Effects of Epinephrine and Heart Rate on the Ratio , 1973, Circulation research.

[52]  Juanzhen Sun,et al.  Impacts of Initial Estimate and Observation Availability on Convective-Scale Data Assimilation with an Ensemble Kalman Filter , 2004 .

[53]  Alfio Quarteroni,et al.  Analysis of a Geometrical Multiscale Model Based on the Coupling of ODE and PDE for Blood Flow Simulations , 2003, Multiscale Model. Simul..

[54]  N. Stergiopulos,et al.  Determinants of stroke volume and systolic and diastolic aortic pressure. , 1996, The American journal of physiology.

[55]  Francesco Vicario,et al.  Cardiovascular system identification: Simulation study using arterial and central venous pressures , 2015, 2015 37th Annual International Conference of the IEEE Engineering in Medicine and Biology Society (EMBC).

[56]  J. Eldredge,et al.  A regulated multiscale closed‐loop cardiovascular model, with applications to hemorrhage and hypertension , 2018, International journal for numerical methods in biomedical engineering.