1 / 29

Reservoir Modeling with EnKF: A sensitivity study

Reservoir Modeling with EnKF: A sensitivity study. L. Heidari (IFP), V. Gervais (IFP), M. Le Ravalec (IFP), H. Wackernagel (ENSMP) Ecole Des Mines de Paris (ENSMP) Institut Français du Pétrole (IFP). 4 th Ensemble Kalman Filter Workshop, 22-24 June 2009 - Bergen , Norway. Overview.

angelo
Download Presentation

Reservoir Modeling with EnKF: A sensitivity study

An Image/Link below is provided (as is) to download presentation Download Policy: Content on the Website is provided to you AS IS for your information and personal use and may not be sold / licensed / shared on other websites without getting consent from its author. Content is provided to you AS IS for your information and personal use only. Download presentation by click this link. While downloading, if for some reason you are not able to download a presentation, the publisher may have deleted the file from their server. During download, if you can't get a presentation, the file might be deleted by the publisher.

E N D

Presentation Transcript


  1. Reservoir Modeling with EnKF: A sensitivity study L. Heidari (IFP), V. Gervais (IFP), M. Le Ravalec (IFP), H. Wackernagel (ENSMP) Ecole Des Mines de Paris (ENSMP) Institut Français du Pétrole (IFP) 4th Ensemble Kalman Filter Workshop, 22-24 June 2009 - Bergen , Norway

  2. Overview • Background • Reservoir Study and History Matching • Kalman Filters • Ensemble Kalman Filter – EnKF • Application • 2D synthetic case study • Sensitivity Tests • Conclusions & perspectives

  3. Overview Background Reservoir Study and History Matching Kalman Filters Ensemble Kalman Filter – EnKF Application 2D synthetic case study Sensitivity Tests Conclusions & perspectives 3

  4. Fluid Flow Simulation Upscaling (Phi, K) Reservoir model Geological model (facies, porosity, permeability) Observed Production data Simulated production data Reservoir Study ? = • Aim • Constrain reservoir models to production history by adjusting parameters: porosity, permeability, etc.

  5. Sequential Methods Time t0 t1 t2 Measurements Measurements Assimilation method Assimilation method Initial state Corrected state (t1) Corrected state (t2) … Kalman Filter Forecasted state (t2) Forecasted state (t1) Dynamical model Dynamical model t0 t2 t1 Time

  6. Overview Background Reservoir Study and History Matching Kalman Filters Ensemble Kalman Filter – EnKF Application 2D synthetic case study Sensitivity Tests Conclusions & perspectives 6

  7. Dynamic Model Operator TrueState Measurement error Measurement Operator Kalman filter • Assumptions • Linear dynamic model • Gaussian probability distribution function • First Order Markov processes • Independent and uncorrelated model and measurement errors • Kalman Filter estimate is based on • Forecast Estimate • Measurements

  8. Distance between  and f Distance between  and d Kalman gain Simulated data Kalman filter analysis • Kalman filter estimate, a , minimizes • Kalman filter estimate, a :

  9. Kalman filter forecast • System Dynamics • Linear • Non linear • Slightly Non Linear Models • Extended Kalman Filter • Taylor expansion of nonlinear model operator F • Highly Non Linear Models • Ensemble Kalman Filter • Ensemble representation of error covariance matrix Model operator

  10. Overview Background Reservoir Study and History Matching Kalman Filters Ensemble Kalman Filter – EnKF Application 2D synthetic case study Sensitivity Tests Conclusions & perspectives 10

  11. Best-guess Ensemble Kalman filters (EnKF) • Ensemble representation of error covariances for • Model State • Measurements • Forecast • The time evolution of uncertainty is obtained by integrating each model state forward in time • Analysis • Kalman linear analysis scheme applied to each model state with the sample covariance

  12. Advantages of EnKF for History Matching Only one fluid-flow simulation over the whole production history per model of the ensemble Real-time assimilation of data, for any data frequency Uncertainty quantification Practical implementation Can be plugged with any fluid-flow simulator Suitable for parallel computation Relatively low computational cost 12

  13. Overview Background Reservoir Study and History Matching Kalman Filters Ensemble Kalman Filter – EnKF Application 2D synthetic case study Sensitivity Tests Conclusions & perspectives 13

  14. Application : 2D synthetic case P P I P P • Reservoir • Size : 500 × 500 × 1 m • Grid : 50× 50 × 1 • PetrophysicalProperties • porosity = cte , Kv/Kh = cte • Reference ln(Kh): spherical variogram • FluidFlow • oil-water, IFP simulator: PumaFlow • Water Injection • Well plan • 5 wells: 4 producers and 1 injector • Production data • Bottom hole Pressure • (BHP) • Surface Oil Rate • (QOS) • Watercut • (WCT) 14

  15. Reference data P2 P3 Inj P1 P4 Ref ln(Kh) 15

  16. Assimilation with EnKF Static parameter Log of permeability One value per grid block Pressure Dynamical state Water saturation Surface oil rate One value per well Bottom hole pressure Simulated measurements Water cut • Ensemble of 50 models • State-parameter vector • 17 assimilations, Once every month 16

  17. Evolution of ln(Kh) mean and variance ln(Kh) reference model Mean Initial Model 1st assimilation 4th assimilation 10th assimilation 17th assimilation Variance 17

  18. Simulated answer and prediction – Prod3 Ensemble members Ensemble mean Reference Initial ensemble Prediction Final ensemble 18

  19. Root Mean Square Error - RMS Reference value at each grid Ensemble Mean value for each grid Total number of grids Simulated data for ensemble member Number of data per time step Number of time steps Reference model Measurements 19

  20. Overview Background Reservoir Study and History Matching Kalman Filters Ensemble Kalman Filter – EnKF Application 2D synthetic case study Sensitivity Tests Conclusions & perspectives 20

  21. Effect of Initial Ensemble True Oil Cumulative Oil Production • 10 Ensembles of size 50 • Prediction of cumulative oil production after 5.5 years 21

  22. Effect of Assimilation Time Interval 2 Months Long Assimilation Interval Base Case Short Assimilation Interval Time 22

  23. Effect of Ensemble Size - 1 Initial Mean Final Ensemble Size Reference ln(Kh) Initial Variance Final 50 100 200 500 23

  24. Effect of Ensemble Size - 2 Prediction Ensemble members Ensemble mean Reference 50Members 100Members 200Members 500Members BHP QOS WCT 24

  25. Measurement Uncertainty Mean Initial Model High Case : 10% Low case : 1% Base Case : 5% Variance 25

  26. Overview Background Reservoir Study and History Matching Kalman Filters Ensemble Kalman Filter Theory Implementation Application 2D synthetic case study Sensitivity Tests Conclusions & perspectives 26

  27. Main challenges Finite number of models used to represent uncertainty Quality of initial sampling Variability Preservation Linear analysis for a non linear model Physical inconsistency between parameters and state after analysis Non Gaussian distributions Multimodal distributions are not well preserved (e.g. channel and barrier permeability distribution, facies distributions) Large amount of data 4D seismic 27

  28. Conclusions and perspectives EnKF application to a 2D model and Sensitivity studies Several ensembles of the same size for a sound uncertainty quantification Smaller assimilation time interval : Better Match Lower Variance Larger ensemble size : Better match for production data Variance preservation Larger simulation and computational cost Larger Measurement Uncertainty Variance Preservation Less efficient for matching production data Perspectives Facies History Matching Real field case study 28

  29. Thank you for your attention ! Any Questions ?

More Related